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ABSTRACT 

A micro computer program to solve the complex task of 
aii»_borne navigation v/as developed to demonstrate the 
practicality of replacing costly general purpose digital 
computers with relatively inexpensive dedicated micro- 
computers on-board naval aircraft . The microcomputer 
program showed that microcomputers have sufficient speed 
and accuracy to solve the navigation problem. In order 
to overcome the microcomputer’s major deficiencies 3 speed 
and accuracy 3 special arithmetic subprograms based on 
table look— up were developed to trade inexpensive memory 
for more speed. An application of graph theory in the form 
of process graphs v/as made to facilitate the development 
and documentation of the navigation program. To aid in 
the testing of the microcomputer program, a Fortran simu- 
lation program was developed to confirm the results of an 
error bound analysis of the navigation program. 
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I. INTRODUCTION 



Current naval aircraft are depending more and more on 
airborne digital computers. The digital computer has proven 
to be a great aid, providing the precision and speed required 
to perform many calculations. The digital computer is very 
versatile, because the same computer can be programmed to 
perform new or different tasks. The only limitations to 
the use of digital computers are those of cost and 
maintenance requirements. 

The digital computers currently used by such naval 
aircarft as the P3C, A6E, A7E, and S3A are all large general 
purpose computers. These computers are very fast, flexible, 
but expensive. The large size, large power requirements, 
and great cost has limited each aircraft to one such com- 
puter. The complexity of these computers has made main- 
tenance difficult. The large cost of each unit makes spares 
prohibitive, therefore one computer going down results in 
the operational loss of one aircraft. 

It is advantageous to utilize several small distributed 
computers to meet all of the systems requirements. In 
order to minimize the probability of the entire system 
failing due to the failure of a single component in a cri- 
tical computer circuit, completely separate computers could 
be used for the various system requirements with back up 
computers ready to fill in when needed. In addition to 



reliability improvements, a distributed approach offers 
the possibility of using less complex equipment, flexibility 
in matching equipment to system requirements, and increased 
standardization. The limitation of this approach is the 
increased cost, weight, size and the complexity of system 
interconnections . 

The creation of the microcomputer , using new developments 
in the Large Scale Integration (LSI) technology, has made 
the distributed computer system possible. The microcomputer 
is a complete general purpose computer on a set of four 
standard LSI chips. The LSI chip measures 200 mils by 200 
mils, requires less than one watt of power and costs about 
$30. The processors Central Processing Unit (CPU), constructed 
on a single chip, is designed to be used in a multiprocessor 
environment . The limitations of a microcomputer are a 
limited instruction set and slow speed. 

This thesis is a design study of the possibility of 
using the MCS-4 microcomputer as the Avionics Navigation 
Computer in a complex navigation system. Section II dis- 
cusses what a microcomputer is, what microcomputers are 
currently available, and what is the system makeup of the 
MCS-J| microcomputer used in this report. Section III des- 
cribes the navlgatin systems currently available, what their 
advantages and limitations are, and the navigation system 
chosen for this report . The navigation equations used to 
compute the current position of the system's carrier are 
discussed in Section IV. The actual programming of the 
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microcomputer is discussed in detail in Section V. This 
includes the programming aids developed in this report, a 
program analysis of the requirements of this system, a 
detailed discussion of the executive routine and the sub- 
routines used in this program, and an error bound analysis 
of the final program. Section VI describes the FORTRAN 
simulation program written to aid in the writing, debugging, 
and testing of the final navigation program. The conclusions 
of this report are summed up in Section VII. The listing 
of both the MCS-^l microcomputer navigation program and the 
FORTRAN simulation program - are included in the Appendices. 
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II . MICROCOMPUTER 



A microcomputer is a general-purpose digital computer 
constructed from a set of LSI chips. It has a complete 
instruction set and is capable of addressing sizeable memo- 
ries. It can interface with a full complement of input and 
output devices. 

The main component of a microcomputer is the micro- 
processor. The microprocessor is a CPU on a chip that 
interprets and executes instructions in a bit-parallel 
fashion. Included on the CPU chip are the Index registers. 
Arithmetic unit and Input/Output control logic. 

The main feature that distinguishes a microprocessor 
from a general purpose or minicomputer CPU is that the 
entire CPU is on one chip. This has been made possibly by 
Large Scale Integration, where over 14,000 Metal-Oxide Semi- 
conductor transistors can be put on one chip. A minicomputer 
CPU with the capability of one microprocessor would require 
over 100 TTL packages. The major advantages of micropro- 
cessors are low cost, low power requirement, and less com- 
plexity in system design. Figure 1 is a photomicrograph 
of the 4004 CPU chip used in the MCS-4 Microcomputer. 

The microprocessor becomes a microcomputer when a con- 
trol program, memory, and input/output circuits are added 
to the system. The control program is usually Metal-Masked 
on a Read-Only Memory (ROM) chip, however', a Programmable 
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Figure 1. Photomicrograph of 4 00 *4 CPU 
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Read-Only Memory (PROM) may be used for initial system 
development. The storing of data or variables in a program 
is handled by the Random Access Memory (RAM) chip which is 
a read/write type memory. The input/output ports are incor- 
porated on the RAM, ROM, or CPU chip or any combination 
of these chips. The total chip area of a typical micro- 
computer is 200 mils by 200 mils and- costs less than $200. 

Microcomputers have been used in specialized, single 
user data processing systems and as components in digital 
products. Microcomputers are currently being used in cal- 
culators, special purpose terminals, measurement systems, 
intelligent traffic controls, small business computers, and 
digital cash registers. Research is in progress to use 
microcomputers in intelligent peripherals, multiplexors 
and communications controllers, automotive control systems, 
and educational systems. Microcomputers can also be used 
in large computer systems to relieve the large central 
processor of much of the overhead associated with scheduling, 
text editing, and file management. 

Microcomputers major application will be in the area of 
dedicated computations. They provide all the data processing 
power needed for these applications. These areas include 
computational tasks required by NTDS, TSC , and most current 
naval aircraft. 

The following is a list of advantages of microcomputers: 

1. Reduced costs due to reduction in number of logic 
card types. 
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2. Self test capability. 

3. Equipment modularity. 

4. Equipment commonality. 

5. Ease in design changes through microprogram. 

6. Reduced logistic support. 

7. Standardization of peripheral interfaces. 

8. Multiprocessor capability. 

The total savings in cost, power required, and size, 
together with the flexibility inherent in the system, make 
the microcomputer a powerful tool for system design of 
dedicated computational tasks. 

A. COMPARISON OF MICROCOMPUTER WITH GENERAL PURPOSE COMPUTER 

The general purpose computer and microcomputer represent 
the two extremes on the computer scale. The general purpose 
computer is large in size, very expensive, requires much 
power, and is very fast. The microcomputer is small, inex- 
pensive, requires little power, and is relatively slow in 
comparison. The general purpose computer is built using 
Integrated Circuits (I.C. Technology), while the microcom- 
puter is built from Large Scale Integration Techniques (LSI 
Technology). The great cost savings of LSI Technology in 
computer design can best be shown by comparing the cost of 
designing and manufacturing a 3000 gate logic unit using 
I.C. techniques versus the LSI techniques. 

The cost of manufacturing a computer system can be 

f 

expressed by the following formula: 
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CS = S [ (CM(i) + CT( i) )V( i) + CG(i)] 
i=l 



where 

CS E Cost of System 

CM(i) E Cost of Manufacturing the i-th Replaceable 
Module 

CT(i) = Cost of Testing the i-th Replaceable 
Module 

V(i) = Number of Modules used in the system of 
the i-th type 

CG(i) = Cost of generating the i-th Replaceable 
Module, including design, layout, test 
condition generation, etc. 

The cost of manufacturing the computer system from I.C. 
circuits was estimated using figures from people knowledgeable 
in this area, however, documentation of these figures was 
not available. 

The cost of manufacturing and the cost of testing an 
individual I.C. module was assumed to be independent of the 
type of I.C. module because of standard means of manufac- 
turing and testing these units. This assumption simplifies 
the "cost of system" equation as follows : 



CM(i) s CM 
CT(i) = CT 
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n 

CS = (CM + CT) V + l CG (i) , 

i=l 



where 



n 

V = l V(i) . 
i=l 

I.C. technology does not permit highly complex functional 
circuitry to be placed on one module. The typical module 
consisting of sixteen pins contains eight gates. A 3000 
gate logic pattern would therefore require V = 360 modules. 

The technology typically used by computer manufacturers 
permits 60 modules to be assembled on one panel which pro- 
vides all the interconnections between the modules. The 
system would then require six panels at an estimated cost 
of $100 per panel. 

In determining the number of distinct modules used, it 
was assumed that about one-half of the modules used in the 
system are duplicates of previously designed modules, n - 180. 
The cost of designing one new module was estimated to be 
$230. The total cost of one system was computed as follows: 

CM = $20 

CT = $ 5 
v = 360 

CG ( i ) = $250 
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n = 180 



CS = (CM + CT) + CG(1) + 600 

CS = 9,000 + 45,000 + 600 = 54,600 

To produce this unit In large numbers, the cost of manu- 
facturing and testing the modules and panels is multiplied 
by the number of units produced while the cost of designing 
remains unchanged. The cost of 1000 units could be predicted 
as follows : 



CS(1000) = 9000(1000) + 45,000 + 600(1000) 

Total Cost of 1000 units =$9,645,000 

Next the cost of developing and manufacturing a 3000 
gate logic unit using LSI technology was investigated. The 
cost estimates for this analysis were obtained from a trip 
to Intel Corp., Santa Clara, California. The LSI design 
cost of the 8008 CPU chip used in the MCS-8 Microcomputer 
was used as an indication of the cost associated with manu- 
facturing the required logic unit by LSI technology. 

The entire design and layout of the 8008 CPU was done 
exclusively by hand and took a total of four man-years. 

The approximate cost of developing this chip is estimated 
at $100,000. The total manufacturing and testing costs of 
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the chip is $15. Since all the logic gates are on one chip, 
there are no other costs associated with the unit. The 
cost of one system is therefore: 

CS = (15) (1) + 100,000 + 0 
CS = $100,015 

Again, the cost of producing the system in large quan- 
tities, the design cost remains fixed while the manufacturing 
and testing costs are multiplied by the number of units. 

The cost of 1000 units using LSI technology is as follows: 

CS(1000) = (15) (1000) + 100,000 
Total cost of 1000 units = $115,000 

It was concluded from these calculations that LSI 
technology could produce the system in quantity about 80 
times cheaper than by I.C. technology. 

It was concluded from this analysis that the microcom- 
puter has a tremendous cost advantage over general purpose 
computers in systems produced in large quantities with a 
need for a dedicated computer. The object of this thesis 
is to demonstrate that the microcomputer also has the speed 
and computational power to handle complex tasks. 



B. COMPARISON OF MICROCOMPUTER WITH MINICOMPUTER 

The minicomputer was developed to meet the need for 
dedicated systems to handle data processing requirements 
such as communications control, data acquisition, and small 
business accounting. The microcomputer is proving itself 
capable of handling these tasks with improved price/perfor- 
man c e , c o mpac t ne ss, and relia b i 1 i t y . 

The greatest advantage the minicomputer has over the 
microcomputer is speed. The LSI chips used for microcom- 
puters are made using Metal Oxide Semiconductor, MOS, tech- 
nology. MOS technology allows smaller size of individual 
transistors, and logical circuits together with low power 
consumption. The electrical properties of MOS circuitry, 
however, make it slower than the Bipolar circuitry used in 
mini computers . In order to increase the computational spe ed 
of microcomputers, many microcomputer makers are switching 
from the slower P-channel MOS devices to the much fas ter 
N- channel MOS devices. By using microprogramming techniques 
such as pipelining, the computati onal_ speed of microcomputers 
can become as fast or faster than minicomputer >-> « 

The next advantage the minicomputer' has over the micro- 
computer is in the instruction set size. The first genera- 
tion microcomputers were limited to between *10 to 60 instruc- 
tions while the minicomputers had instruction sets in the 
100-120 range. The gap here is also closing with the second 
generation microcomputers having instruction sets of 50-100 
instructions . 
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The last major advantage of the minicomputer is the 
existing software developed over the past several years to 
assist program development. With the development of high 
level programming languages (PL/M for the Intel 8008 and 
8080), the microcomputers are quickly eliminating this 
advantage of the minicomputers. 

The microcomputer offers better price/performance, 
lower power consumption, smaller size and higher reliabil- 
ity than minicomputers. Although a single microcomputer 
can not match the power of a minicomputer, several micro- 
computers can be combined to share the workload at a cost 
still less than a minicomputer. An added advantage of 
multi-microcomputers is processor reliability which can be 
increased through the use of back-up processors, providing 
a self-test and repair capability. 

Table I is a summary of the comparison of microcomputers 
and minicomputers. 

C. CURRENTLY AVAILABLE MICROCOMPUTERS 

The number of microcomputers being designed and currently 
available is increasing rapidly. This section covers eight 
microcomputers that are currently available or will be in 
the near future. Table II is a summation of the capabilities 
and support available for each microcomputer. 
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MINICOMPUTER 


MICROCOMPUTER 


CPU Cost 




$1500 


$30 


Instruction Speed 




2-5 msec 


10-20 msec 


Execution Speed 
(Memory to Memory 


Add) 


5-20 msec 


15-60 msec 


Instruction Set 




100-150 


50-100 


Price/Instruction 


Ratio 


$80-$300 


$40 


Registers 




1-30 


16-100 


Price/Register Ratio 


$300 


$15 


Memory Size 




64K 


64K 


Performance/Price 
.Word Length 
^Add Time x Price 




200 (Bits 
Sec-Dollar ) 


100 (Bits 
Sec-Dollar) 


Power Consumption 




4 Watts 


1 Watt 


Reliability (CPU) 




Less (Due to 
100 TTL Packages) 


Greater (Due to 
One LSI Chip) 



TABLE I. MICROCOMPUTER VERSUS MINICOMPUTER 
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TABLE II. SUMMARY OF AVAILABLE MICROCOMPUTERS 



1. Intel MCS-4 



The MCS-4 is a four-bit parallel processor with a 
fixed instruction set. The 4004 CPU is a P-Channel MOS LSI 
chip mounted on a l6-pln package. A four-bit data bus 
connects the CPU with up to 16 ROMs and l6 RAMS. The 
instruction cycle is 10.8 microseconds. 

The instruction set used by the MCS-4 consits of 
45 instructions grouped into three categories: Machine, 

Input/output, and accumulator. System development aids 
include a cross assembler available on the IBM 360/67 located 
at the Naval Postgraduate School. The SIMQ4-02 Hardware 
Prototyping System with Assembling, Programmable -ROM 
Programming, and Debugging capability is also available 
at the NPS microcomputer lab. 

The MCS-4 Microcomputer was chosen as the navigation 
microcomputer in this thesis. The MCS-4 is covered in more 
depth in Section II. D. 

2. Intel MCS-8 

The MCS-8 uses the 8008 single CPU chip. The 8008 
is an eight-bit fixed instruction set parallel processor 
mounted in an 18-pin package. The 8008 executes a single 
instruction in 20 microseconds. The 8008 is capable of 
addressing 16K bytes of memory. An eight-bit data bus 
interfaces the processor with memory. A total of 48 instruc- 
tions are available broken into four groups : Instruction 
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Register, Accumulator, Program Counter and stack control, 
and Input/Output. 

System Development aids Include a cross assembler 
available on the IBM 360/67 at the Naval Postgraduate School. 
The SIM 8-01 Prototyping System with programmable and erasable 
ROMs is available in the NPS microcomputer lab for develop- 
ment and check-out of microprograms . The development of 
the PL/M Higher Level Language patterned after PL/I for 
the MCS -8 greatly facilitates the programming task. The 
PL/M compiler is also available on the IBM 360/67 at NPS. 
Lastly the INTELLEC -8 system, available in the NPS micro- 
computer lab, makes available to the programmer a resident 
software monitor, assembler, PROM programmer, and text 
editing capability to aid in the development of micropro- 
grams for the MCS- 8 . 

3. Intel 8080 — 

The Intel 8080 is the first of the second generation 
microcomputers. The 8080 derives more speed and capability 
by using the more efficient N-Channel MOS . The 8080 CPU 
is a single-chip eight-bit parallel processor in a 40 pin 
package. The 8080 contains six eight-bit data registers, 
an eight-bit accumulator, three eight-bit temporary regis- 
ters, four testable flags, and an eight-bit arithmetic/logic 
unit. The execution time for one instruction is 2 micro- 
seconds. 

The 8080 can directly access 6HK bytes of memory. 

A separate 16-bit address bus is provided as well as ten j 
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control lines that indicate CPU and I/O bus status. Up to 
256 I/O devices can be directly addressed. Multiprocessor 
capability is designed into the MCS-8080. 

The 8080 is software compatible with the 8008 micro- 
processor. The 8080 instruction set contains the 48 instruc 
tions of the 8008 plus 26 new instructions for a total of 74 
The 8080 is capable of unlimited subroutine nesting. 

System development aids include a cross assembler, 
INTELLEC 8080 simulator, and the PL/M higher level program- 
ming language . 

4. AMI 7300 ~ ' 

The AMI 7300 is an eight-bit fully parallel, bus- 
oriented processor. The processor consists of two chips, 
the Micro instruction ROM chip (MIR) and the registers- 
adder logic unit chip (RALU). Both chips are P-Channel MOS 
and each is packaged in a 40-pin dual-inline package (DIP). 

The MIR contains a mask-programmable 512 word x 
22 bit microinstruction ROM and a programmable instruction 
mapping array, allowing up to 50 microprogram locations to 
be predefined for macroinstruction decoding. This allows 
the instruction set to be tailored to suit a particular 
application. A hardware address stack and loop counter 
allows subroutine nesting to seven levels. 

The RALU contains 48 registers which may be utilized 
as one or two first-in/last-out stacks or as one or two 
files of general registers. The eight-bit adder/subtractor 
performs over 30 arithmetic and logical operations. The 
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instruction set consists of three basic formats: Register 

control, literal and branching. The processor can address 
up to 64 K of memory. 

Software development aids include a cross assembler 
and instruction simulator. 

5 . Fairchild PPS-25 

The PPS-25 is a BCD oriented 25-digit .serial/parallel 
processor. The system is best suited to decimal applications 
such as calculators, keyboard/printer interface, and vending 
machines. Four level subroutine nesting and three way condi- 
tional branching are provided. The system can be micro- 
programmed with a custom instruction set to best meet the 
needs of the user. 

Program storage consists of up to 26 ROMs, each 
capable of storing 256 twelve-bit words. Seven general 
purpose 25-digit registers are provided and an external 
interrupt capability is included. 

The 3805 Arithmetic Chip includes the Adder/Subtractor 
plus a 25-digit register. The instructions are located in 
the 3810 ROM. A total of 30 arithmetic/logic and l6 I/O 
instructions are available with the standard set. The 1/0 
format permits expansion to 63 I/O instructions. 

Software development aids include a cross assembler 
and instruction simulator. 

6 . National Semiconductor 0PCP/IKP -I6 

The GPCP sixteen-bit microcomputer processor consists 
of five MOS LSI chips, each mounted in a 40-pin DIP. The 
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five chips consist of four RALU chips and one CROM (Control 
Read Only Memory). Each RALU chip is a four-bit slice of 
CPU with its own registers, ALU logic, and I/O data lines. 
Multiple CROMs may be used to increase the size of the 
instruction set. 

The RALU consists of a four-bit program counter, 
four-bit memory data register, four-bit memory address 
register, four-bit accumulator, a pushdown stack, data 
multiplexer, and four-bit arithmetic and logic unit. The 
system may be expanded to 32 bits by adding four more RALU 
chips . 

The CROM contains the control instructions for the 
RALU chips. The CROM is broken into two parts, an instruc- 
tion ROM containing 100 twenty-three-bit words and an address 
control ROM consisting of 12 programmable ten-bit words. 

The IMP-16 is a sixteen-bit microcomputer developed 
to use the GPCP microprocessor. The CROM provided with the 
system contain a *13 word instruction set. The instruction 
set may be expanded to meet the system designer's specific 
needs. Communication between the RALU chips and the CROM 
chip is over a sixteen-bit data bus and a four-bit control 
bus. This requires sixteen extra TTL packages excluding 
memory and timing. 

Software development aids include a cross assembler 
and prototyping system with resident monitor, assemblers, 
and linking loader. 
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7. North American Rockwell PPS-4 



The PPS-4 microcomputer consists of a set of six 
MOS LSI chips. The CPU is a four-bit single chip processor 
mounted in a 42-pin package. The other five chips support 
the CPU and consist of: A 256 x 4-bit RAM, 1024 x 8-bit 

RAM, an I/O buffer, and a two-phase clock generator. 

The CPU can drive up to 4K bytes of ROM and 4K 
bytes of RAM over its 12-bit parallel address bus. The 
basic instruction set contains 50 instructions with an 
execution time of five microseconds. 

Twenty-one multiplexed lines interconnect the CPU 
with ROM, RAM, and I/O circuits. These lines are func- 
tionally grouped into twelve parallel address lines, eight 
parallel data lines, and one write command and I/O enable 
line. The ROM has two chip-select inputs and the RAM has 
one chip-select input, which may be directly controlled by 
discrete outputs from the CPU or I/O circuits to expand on 
memory without the need for auxiliary circuits. Each I/O 
chip can handle up to 12 inputs and 12 outputs with a total 
system capability of 16 I/O chips. 

Software development aids include a cross assembler 
and simulator available on a national time-sharing network. 

8. Signetics PIP 

The programmable integrated processor (PIP) is a 



single-chip eight-bit CPU made with N-channel MOS technology. 
This second generation CPU is packaged in a 40-pin DIP. The 



address logic, control memory, and ALU are organized around 
an eight-bit bidirectional data bus. There are fifteen 
address lines for handling external memory and I/O circuitry 
The PIP instruction set contains over 64 instructions with 
an execute time of less than ten microseconds for the most 
complex instruction. 

The PIP chip can be broken into four parts, the 
address logic, the RALU, and the control section. The 
address logic section handles all instructions and includes 
a return address stack that allows subroutine nesting to 
eight levels. The RALU section contains four 8-bit general 
purpose registers and executes all arithmetic, boolean, 
compare, and rotate operations. The control section manages 
operation of all external control lines, decodes all instruc 
tions, and coordinates the activities of all other internal 
circuitry. 

Software development aids include a cross assembler 
and instruction simulator. 

9. Summary 

Microcomputers have the capability of replacing 
both special function logic modules and large computational 
machines. Certain microcomputers are more suited for one 
application than for others. In selecting a microcomputer, 
such parameters as data word length and type, instruction 
power, and interface structure must be considered. Table II 
is a summary of the microcomputers discussed in the section. 
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The rapid developments in the LSI technology should 
bring about many new improvements in microcomputers such as 
wider word lengths, larger memory capacity, and more flexi- 
ble and convenient I/O interfacing. Improvements in LSI 
structure will make possible single chip microcomputers with 
the CPU, I/O, and memory all on one chip. Work on a Bipolar 
LSI processor could develop into a microcomputer with more 
speed and memory than most present minicomputers. 

D. THE MCS-4 MICROCOMPUTER 

The MCS-JJ was the first microcomputer made. Compared 
to other microcomputers, it is slow with a limited instruc- 
tion set. The most complex function available on the MCS-JJ 
is a four-bit add. The MCS-4 is also the least expensive 
microcomputer. The structure of the MCS-4 is similar to a 
general purpose computer, making it compatible with the 
requirements of a navigation computer. 

The MCS-JJ was chosen as the microcomputer in this design 
study for two major reasons. The first reason is that the 
MCS-JJ is the least powerful and hence serves as a lower 
bound of the microcomputers. To prove that the MCS-^ is 
capable of handling the required navigational computations, 
would in itself prove microcomputers capable of handling 
complex tasks. The second major reason for choosing the 
MCS-*4 microcomputer is that it is available, has been tested, 
and has the required software aids to complete a design study. 
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1 . MCS-^ System Description 



The MCS-4 microcomputer is built up from a standard 
set of off-the-shelf chips. The only custom part is the 
ROM chip which stores the specific program defined by the 
user and requires a metal mask option for each new program. 

The MCS-*J consists of four chip types, each packaged 
in a conventional l6-pin DIP: 

(1.) A Central Processor Unit Chip-CPU- 
(2.) A Read Only Memory Chip-ROM-^OOl 
(3.) A Random Access Memory Chip-RAM-^1002 
(i|.) A Shift Register Chip-SR-4003 

The CPU contains the control unit and the arithmetic 
unit. The ROM stores the program and data tables, the RAM 
stores input data and variables, and the Shift Register is 
used in conjunction with I/O devices to effectively increase 
the number of I/O lines. 

A qomplete microcomputer can be built using only 
a single CPU chip and a single ROM, and the only external 
circuitry required is a two phase clock. The CPU is capable 
of driving a system up to 16 ROMs (^K bytes), 16 RAMs (6^10 

bytes) and 128 I/O lines, with no additional interfacing 

circuitry. The CPU communicates w ith the RAM’s and ROM’s 
by means of a four-line data bus. This single data bus is 
used for all information flow between the chips except for 
control signals which are sent to RAM and ROM over five 

additional lines. Figure 2 shows the MCS-'I System Inter- 

connection . 
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Figure 2 



MCS-JJ System Interconnection 



The MCS-4 uses a 10.8 microsecond instruction cycle. 
The cycle is broken into eight steps. In the first three 
steps, the CPU sends the memory address to the ROM in three 
4-bit nibbles. The ROM then sends back 8 bits of instruction 
in two 4-bit nibbles during steps 4 and 5. The instruction 
is then interpreted and executed by the CPU during the last 
three' steps . 

a. 4004 CPU 

The heart of the MCS-4 microcomputer is the 4004 
CPU. The 4004 CPU contains the following functional blocks: 
(1.) Address Register and Address Incrementer 
(2.) Index Register 
( 3. ) 4-bit adder 

(4.) Instruction Register, Decoder and Control 
(5.) Peripheral Circuitry 

The Address Register is a RAM array of 4 x 12 
bits. One level is used to store the current instruction 
address, leaving three levels to store the addresses of 
nested subroutines. As each byte of address is sent onto 
the data bus, the address is incremented by a 4-bit carry 
look-ahead circuit. The incremented address is then trans- 
ferred back into the address register. 

The Index Register is a RAM array of 16 x 4 
bits and has two modes of operation. In one mode of opera- 
tion the index register provides sixteen directly addressable 
storage locations. In the second mode, the index registers 
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provide eight pairs of addressable storage locations for 
addressing RAM and ROM as well as for storing data fetched 
from ROM. The index registers can thus provide 64 bits of 
RAM to a minimum MCS-4 system of one CPU and one ROM. 

The 4-bit Adder is of the ripple-through carry 
type. The output of the adder is transferred to the accumu- 
lator and carry flip-flop. The accumulator is provided with 
a shifter to implement rotate right and rotate left instruc- 
tions. The accumulator also communicates with the command 
control register, the condition flip-flop and the 4-bit 
internal data bus. The condition logic allows the execution 
of conditional instructions based on the contents of the 
accumulator, index registers, or the status of the control 
lines . 

The Instruction Register is an 8-bit register 
which is loaded with the two 4-bit nibbles of instruction 
read from the ROM. The instructions are decoded in the 
instruction decoder and appropriately gated with timing 
signals to provide the control signals for the various 
functional blocks. 

The peripheral circuitry consists of the 4-bit 
internal data bus, the timing and SYNC generator, one ROM 
command control and the four RAM command control output 
buffers, and the reset flip-flop. 
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b. 4001 ROM 

The 4001 is a 204 8-bit metal mask programmable 
ROM. The 4001 performs two functions. As a ROM, its first 
function is to store 256 x 8-bit words of program or data 
tables. The second function of the 4001 is to act as a 
vehicle of communication between the data bus lines and 
peripheral devices through the 4-bit input-output port 
located on each chip . 

c. 4002 RAM 

The 4002 also performs two functions. As a 
RAM it stores 320 bits arranged as four registers of twenty 
4-bit characters each. As a vehicle of communication with 
peripheral devices, it is provided with 4 output lines and 
associated control logic to perform output operations. 

d. 4003-SR 

The 4003 is a 10-bit serial-in/parallel-out 
serial-out shift register. Its function is to increase the 
number of output lines to interface with I/O devices. 

2 . Instruction Set 

The MCS-4 Instruction Set consists of a total of 
45 instructions grouped into three sets: Machine Instruc- 

tions, Input/Output and RAM Instructions, and Accumulator 
Group Instructions. 

The Machine Instructions are the Housekeeping 
instructions of the MCS-4. They consist of two types of 
instructions: one-word instructions which are 8-bits wide 
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and require one instruction cycle, and two-word instructions 
whieh are 1 6 -bits wide and require two instruction cycles. 
Table III is a list of the Machine Instructions. 

The Input/Output Instructions are used to transfer 
data between the Accumulator and RAM. These instructions 
are also used to transfer data between the Accumulator and 
the I/O ports located on the ROMs and RAMs . Table IV is 
a list of the Input/Output Instructions. 

The Accumulator Instructions are used to perform 
bit-by-bit manipulation of the data in the Accumulator. 

Table V is a list of the Accumulator Instructions. 

3 . System Development Aids 

The program written for the MCS-4 is metal-masked 
into the ROMs. Before doing this, the program needs to be 
tested to insure that it functions correctly in all situa- 
tions. To aid in testing programs written for the MCS-* 1 , 
there are available three system development aids. These 
aids include an assembler and interpreter, a complete hardwar 
prototyping system, and a resident software monitor . 

The MCS- Assembler and Interpreter is an ALGOL-W 
program which can be used to test and debug programs for 
the HCS-lJ microcomputer. The system consists of an assem- 
bler which allows symbolic programming of the MCS-*I ROM 
and an interpreter which simulates the actions of the MCS-*J. 
The program provides extensive diagnostic facilities for 
monitoring the actions of the MCS-'i program. 



(Those instructions preceded by an asterisk (*) are 2 word instructions that occupy 2 successive locations in ROM] 
MACHINE INSTRUCTIONS (Logic 1 - Low Voltage - Negative Voltage; Logic 0 - High Voltage - Ground I 



MNEMONIC 


OPR 

D 3 D 2 D 1 °0 


OPA 

° 3 ° 2°1 °0 


DESCRIPTION OF OPERATION 


NOP 


0 0 0 0 


0 0 0 0 


No operation. 


•JCN 


0 0 0 1 
A j A2 Aj A2 


C 1 c 2 C 3 C 4 
A, A 1 A, A| 


Jump to ROM address A 2 A 2 A 2 A 2 . A, A, A, A j (within the same 
ROM that contains this JCN instruction) if condition Cj 020304^) 
is true, otherwise skip (go to the next Instruction In sequence) 


•FIM 


0010 
°2 °2 °2 °2 


R R R 0 
°1 D 1 D 1 D 1 


Fetch Immediate (direct) from ROM Data 0 2 . Di to Index register pair 
location RRR.'?) 


SRC 


0010 


R R R 1 


Send register control. Send the address (contents of index register pair RRR) 
to ROM and RAM at X 2 and X3 time in the Instruction Cycle. 


FIN 


0 0 11 


R R R 0 


Fetch indirect from ROM. Send contents of Index register pair location 0 
out as an address. Data fetched is placed Into register pair location RRR. 


JIN 


0 0 11 


R R R 1 


Jump indirect. Send contents of reg ( ster pair R R R out as an address 
at At and A2 time In the Instruction Cycle. 


•JUN 


0 10 0 
^2 ^2 ^2 ^2 


a 3 A .1 a 3 A 3 
A 1 A, A 1 A y 


Ju.i.h uiK.undiuor.di 10 Kum address A3, A 2, Aj. 


•JMS 


0101 

A2 A2 A2 A2 


A 3 A 3 A 3 *3 
A 1 A 1 A 1 A 1 


Jump to subroutine ROM address A3, A 2 . A 1 , save old address. (Up 1 level 
In stack.) 


INC 


0110 


R R R R 


Increment contents of register RRRR. <3) 


• IS 2 


0111 

A 2 Aj A 2 A? 


R R R R 
A 1 A 1 A 1 A 1 


Increment contents of register RRRR. Go to ROM address A 2 . A 1 
(within the same ROM that contains this ISZ Instruction) if result # 0 , 
otherwise skip (go to the next Instruction In sequence) 


AOD 


10 0 0 


R R R R 


Add contents of register RRRR to accumulator with carry. 


SUB 


10 0 1 


R R R R 


Subtract contents of register RRRR to accumulator with borrow. 


LO 


10 10 


R R R R 


Load contents of register RRRR to accumulator. 


XCH 


10 11 


R R R R 


Exchange contents of index register RRRR and accumulator. 


661 


1100 


DODD 


6ranch back (down 1 level in stack) and load data OODO to accumulator. 


IDM 


110 1 


0 0 0 0 


Load data 0000 to accumulator. 



Table III 



MCS- 1 ! Machine Instructions 



INPUT/OUTPUT AND RAM INSTRUCTIONS 

(The RAM'iand ROM’s operated on in the I/O and RAM instructions have been previously selected by the last SRC instruction executed.) 



MNEMONIC 


OPR 

D3O2O1D0 


°3 


OPA 

°2°i 


°0 


DESCRIPTION OF OPERATION 


WRM 


1110 


0 


0 


0 


0 


Write the contents ol the accumulator into the previously selected 
RAM mam memory character. 


WMP 


1110 


0 


0 


0 


1 


Write the contents of the accumulator into the previously selected 
RAM output port. (Output Lines) 


WRR 


1110 


0 


0 


1 


0 


Write the contents of the accumulator into the previously selected 
ROM output port. (I/O Lines) 


WPM 


1110 


0 


0 


1 


1 


Write the contents of the accumulator into the previously selected 
half byte of read/write program memory (for use with 4008/4009 only) 


VYR£ ,4) 


1110 


0 


1 


0 


0 


Write the contents of the accumulator into the previously selected 
RAM status character 0 


WR1 <41 


1110 


0 


1 


0 


1 


Write the contents of the accumulator into the previously selected 
RAM status character 1. 


WR2 ,4> 


1110 


0 


1 


1 


0 


Write the contents of the accumulator into the previously selected 
RAM status character 2 


WR3 ,4) 


1110 


0 111 


Write the contents of the accumulator into the previously selected 
RAM :•.*!•;' 3 


SBM 


1110 


1 


0 


0 


0 


Subtract the previously selected RAM mam memory character from 
accumulator with borrow 


ROM 


1110 


1 


0 


0 


1 


Read the previously selected RAM main memory character 
into the accumulator. 


ROR 


1110 


1 


0 


1 


0 


Read the contents ol the previously selected ROM input port 
into the accumulator. (I/O. Lines) 


AOM 


1110 


10 11 


Add the previously selected RAM main memory character to 
accumulator with carry. 


rd*< 4 > 


1110 


1 


1 


0 


0 


Read the previously selected RAM status character 0 into accumulator. 


roi ,4) 


1110 


110 1 


Read the previously selected RAM status character 1 into accumulator. 


RD2 ,4) 


1110 


1110 


Read the previously selected RAM status character 2 into accumulator. 


R03 ,4 > 


1110 


1111 


Read the previously selecied RAM status character 3 into accumulator. 



Table IV 



MCS-4 Input/Output Instructions 



accumulator group instructions 



iac 

CMC 



1 1 1 1 



1 1 1 1 



111 ' 



0 0 0 0 
0 0 0 1 
0 



0 0 



0 0 1 



0 10 0 



Clear both. (Accumulator and carry) 



Clear carry. 



Increment accumulator. 



Complement carry. 



Complement accumulator. 



0 10 1 



1111 



1111 



0 110 



10 0 0 



Rotate left. (Accumulator and carry) 



(Accumulator and carry) 



Rotate right 



Transmit carry to accumulator and clear carry. 



Decrement accumulator 



* tubuact and dear carry. 



10 10 



1 0 0 



Decimal adjutl accumulator 



Keyboard procett 

one out of lour c ode to a binary codc_ 



Convert! the content ol the accumulator from a 



Designate command line. 



notes 



(llfhe condition code It attigned at follow* 

Cj - 1 Invert jump condition 
C y . 0 Not Invert lump cond.tion 
WRRR it l*e of 1 of 8 inde* register 



Cj * 1 Jump H accumulator It *ero 
- 1 Jump If carry/MnV it a 1 

palrt In the CPU. 



C 4 - 1 Jump If tett tlgnal It a 0 



RR „R „ ,h, .<**... o. . o. 16 in*. "* CPU ch.. ...... .nd 4 ...lu. 

t *<n ram cn,P "«* ro. ...p «'<««. - "**"• h0 — 

ct.pnomb... RAM ^ 

locnon. ... «l.c..d by Ih. In.lrucllon cod. IOPAI. 



Table V. MCS-1) Accumulator Instructions 
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There are two hardware prototyping systems avail- 
able. Each system uses the electrically programmable and 
erasable ROM (PROM). The 1701 or 1702 PROM is programmed on 
the prototyping system and simulates the exact action of 
a programmed ROM. The SIM^— 02 prototyping board consists 
of a CPU and is capable of programming and testing up to 
sixteen 1701 or 1702 PROMs. The HOOH CPU also controls 
sixteen RAMs together with TTL simulations of eight ROM 
output ports and eight ROM input ports. The SIMH-01 is 
designed as a prototype for small systems. This board 
contains provisions for up to four 1701s or four 1702s. 

It -also provides up to four RAM output ports, four ROM 
ourput ports and four ROM input ports. For small quantity 
systems where the cost of designing a metal-masked ROM can 
not be justified, the 1701 or 1702 PROM can be used instead. 

The recently announced INTELLECT system makes avail- 
able to the programmer a resident software monitor, assem- 
bler, PROM programming, and text editing capability. 

4 . Costs 

The cost of developing a system that uses the MCS-*I 
Microcomputer can be divided into two areas: Hardware and 

software. The hardware costs have been shown to be small 
for the microcomputer. This is a list of the MCS-4 hardware 
costs : 
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(1) 


(over 100) 


(over 1000) 


CPU 


$60 


$30 each 


$15 each 


ROM 


$60 


$15 each 


$ 5 each 


RAM 


$30 


$15 each 


$ 5 each 



The software costs of programming the microcomputer 
are not really known. The assembler type language used to 
"'program the -MCS-4 requires the programmer to keep track of 
the contents of all registers as the program steps through 
its set of instructions. Even one bit, such as the carry 
bit, can cause large errors to be produced if not accounted 
for. All indexing and transferring of data between the 
processor and the assigned locations in memory must be 
written into the program. 

The task of writing a small program on the MCS-iJ 
is fairly straight forward, however, as the computational 
complexity grows, the task of writing the program increases 
rapidly. 

The sections following this discuss the problems 
encountered in programming the MCS-^I to handle the computa- 
tional tasks of a large navigation system. Aids in software 
development were investigated and are also covered. 



III. NAVIGATION SYSTEM 



Air Navigation is the process of directing the movement 
of an aircraft from an initial point to a desired final 
point. A coordinate reference frame must be established 
and the initial point and the final point must be located 
in the reference coordinates. A Navigation System must 
provide timely coordinate measurement and computation of 
desired aircraft position. 

The prime navigation problem years ago was to reach the 
desired destination. Today many aircraft wish to use the 
same airspace at the same time. This problem has required 
separation standards to be put into effect. Aircraft arc 
assigned a block of space, measured by lateral, longitudinal, 
and vertical dimensions, which moves at the speed indicated 
in the pilot's flight plan. It is the pilot's responsibility 
to remain within this block of space. As more aircraft 
wish to use the same airspace, the separation standards will 
have to be reduced. The desired navigation system is one 
that gives a continuous, real time indication of where the 
aircraft is located. Navigation systems must be required 
to have greater accuracy, greater automation, and simplicity 
in operation and display. 

The need for a navigation system to be reliable is a 
necessity. Equipment failure must not endanger the aircraft. 



The reliability of a navigation system can not be measured 
solely by the fault-free operation of one piece of equipment. 
Reliability is a function of the operation of the total navi- 
gation system. In case of partial equipment failure, the 
navigation system must automatically switch to an alternate 
source of information. Therefore, decision making circuitry 
must also be a part of the navigation system to insure a 
continuous flow of accurate navigation information. 

The following is a review of current navigation systems 
available. Prom this information, the possibility of inte- 
grating the MCS-4 microcomputer into a total navigation 
system will be investigated. 

A. SELF-CONTAINED NAVIGATION SYSTEM . 

Only automated navigation systems can satisfy the require- 
ment for continuous and accurate navigation information. 

The need is for a pictorial type of display which will give 
the pilot an accurate and immediate indication of the air- 
craft’s present position. Economic constraints have en- 
couraged that this navigation system be standard throughout 
the military and have a close commonality with civil aviation 
needs. Many of the military and civil aircraft in use today 
have space provided for this type of system. 

The navigation systems in greatest use today are radio 
navigation systems. In order to keep the cost and complexity 
of airborne radio navigation equipment down, a comprehensive 
system of large and costly radio stations must be kept on 
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the ground. Some of the major systems in use today are 
listed in Table VI. When such a system is created, it is 
almost impossible to abolish it. The high initial cost of 
up a new system and the large number of aircraft 
equipped for and dependent on the system in use, make it 
hard for new systems based on new technology to become 
established. 

Radio navigation systems have many limitations. Some 
of these limitations are due to the short distance of ground 
wave propagation, sky wave contamination, atmospheric noise, 
multipath effects, and site error which can result in ambi- 
guous position fixes. Errors in radio navigation signals 
can not be predicted because they are a function of daily 
and seasonal changes in environment, temperature, ionosphere 
location, and local weather. Even the placement and type 
of radio used in the system is affected by political factors. 
These limitations, together with frequency interference 
between radio stations and the susceptability of radio 
navigation systems to jamming, all support the desirability 
of aircraft to have a self-contained navigation system. 

Military aircraft navigation systems require world-wide 
flexibility with the ability to navigate without reliance 
on ground-based aids or the use of equipment susceptible 
to jamming. Civil aircraft that fly over oceanic or desert 
routes require a self-contained navigation system as its 
primary source of position fixing. The development of 
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hovercraft that must cross large bodies of water without 
land based navigation aids and the need for cruise missiles 
that can independently maneuver over hostile territor all 
point to the need for self-contained navigation systems. 

B. AIR DATA SYSTEM 

An air-data system consists of aerodynamic and thermo- 
dynamic sensors inputting to a central air-data computer. 

The sensors measure the characteristics of the air surrounding 
the vehicle and input them to the computer. The computer 
calculates flight parameters such as true airspeed, free- 
stream outside-air temperature, and Mach number. The sensors 
required in this system are angle-of-attack vane, static 
pressure source, pitot tube, and total temperature probe. 

In order to use the air-data computations for dead- 
reckoning, the attitude and heading of the vehicle must be 
supplied from some external source onboard. This can be in 
the form of a simple directional and vertical gyroscope or 
a stable-platform configuration such as an Inertial Navigation 
System. This information together with information on the 
velocity and direction of the wind allows the navigator to 
determine the position of the vehicle by extrapolating from 
a previously known fixed position. 

C. DOPPLER NAVIGATION SYSTEM 

A Doppler Navigation system radiates a pattern of beams 
to the surface of the earth and receives the reflection of 
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this energy back. The difference between the frequency of 
the signals transmitted and the frequency of the signals 
received is called the Doppler Effect and can be used to 
compute the vehicle's velocity along each beam. The 
frequency shift in each of the beams is detected and used 
by the system computer to calculate the distance traveled 
along and across the vehicle's true heading. 

The advantages of the Doppler Navigation System, as 
listed by Kayton and Fried [Ref. 21], are as follows: 

1. ) It provides continuous velocity and position with 

respect to the ground. 

2. ) It is completely self contained. 

3. ) Its average-velocity information is extremely 

accurate . 

h.) Doppler-Radar information is obtainable anywhere 

on earth, including over oceans. 

5. ) It is an all-weather system. 

6. ) Doppler Radars are amenable to high-reliability 

all-solid-state design because of their low 

radiated power . 

7. ) The Doppler Navigator does not require preflight 

alignment or warmup . 

8. ) It radiates at microwave frequencies. 

The disadvantages of the Doppler Navigation System are: 

1. ) It is dependent on an external direction sensor 

for azimuth information. 

2. ) The position Information derived from a Doppler 

Navigator degrades as the distance traveled 

increases . 

3. ) The short-term or instantaneous velocity information 

is not as accurate as the average velocity. 
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D. INERTIAL NAVIGATION SYSTEM 

Airborne inertial navigation systems, hereafter re- 
ferred to as INS, operate on the principle that every time 
the vehicle changes speed or direction it is said to experi- 
ence an acceleration. By measuring this acceleration, the 
velocity and distance traveled by the aircraft can be found 
by integrations. 

The INS is a gyro-stabilized platform with accelerometers 
mounted so that accelerations of the vehicle are measured 
in the north/south, east/west and vertical directions. The 
total system is mounted on gimbals to allow the vehicle to 
rotate without disturbing the attitude of the stable platform. 

The INS must be kep tangent to the earth’s surface at 
all times or the accelerometers w ill experience an accelera- 
tion error due to gravity. This problem is solved by con- 
structing a theoretical pendulum in the system with its 
bob at the center of the earth and its other end on the 
surface of the earth. The aligning of the INS before flight 
will bring this pendulum to rest making it impossible to 
be set in motion by any force at the surface of the earth. 

Any small error in the alignment causes a sinusoidal error, 
called the Schuler Pendulum effect to be introduced in all 
inertials. This effect is periodic with a period of 8 J J 
minutes . 

The Inertial Navigation System has the following advan- 
tages as listed by Kayton and Fried [Ref. 21]: 
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1. ) Its indications of position and velocity are 

instantaneous and continuous. 

2. ) It is completely self-contained. 

3. ) It is nonradiating and nonj ammable . 

4. ) Navigation information is obtainable anywhere on 

the earth. 

5. ) It Is an all-weather system. 

6. ) •-Navigation - information is -independent of the 

vehicle’s maneuvers. 

7. ) It directly provides outputs of position, ground- 

speed, and vertical position. 

8. ) It is the most accurate source of roll, pitch, 

and attitude of the vehicle. 



The Inertial has the following disadvantages: 



1 .) 

2 .) 



3. ) 



The' position and velocity information degrades 
with time. 

The equipment is expensive and relatively difficult 
to maintain and service. 

The INS must be initially aligned. 



E. POSITION FIXING SYSTEMS 

Position fixing is the determination of the vehicle's 
position without reference to any former position. Short 
range radio systems such as ADF, TACAM, DECCA, VOR, and 
LORAN A are not considered applicable for long range 
navigation . 

1 . Long Range Radio Navigation Systems 

The two primary long range Radio Navigation Systems 
in use today arc LORAN C and OMEGA. LORAN C and OMEGA are 
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both hyperbolic fixing systems. The vehicle’s fix position 

is obtained by measuring the relative distance from two or 

more stations of know location. 

L 

a. The LORAN C Network 

The LORAN C network is comprised of master ground 
station and at least two associated slave ground stations. 
Pulsed transmissions radiated from the master station are 
received at the slave stations; each slave station then 
transmits similar groups of pulses, synchronized accurately 
with the signals received from the master station. Receivers 
on-board the vehicle measure the time differences between 
the master and slave's transmissions which allows an automatic 
fix computation. 

The following are advantages and capabilities 

of LORAN C: 

1. ) Currently operational with minaturized receivers 

developed . 

2. ) Range up to 1500 nm. 

3. ) Continuous fixing provided. 

. ) Passive at the receivers. 

5. ) All weather capability 

6. ) Accurate position fix to v >200 ft. 

7. ) Total cost of on-board receiving equipment 

approximately $10,000. 

The disadvantages of LORAN C are: 

1. ) Position fix in hyperbolic coordinates. 

2. ) Dependent on fixed ground stations. 



3. ) Susceptible to atmospheric noise and sky-wave 

contamination . 

4. ) Susceptible to jamming. 

5. ) Total world coverage not available. 

b. The OMEGA System 

OMEGA is an earth reference hyperbolic phase- 
matching navigation system operating at 10.2, 11.33, and 
13.6 kHz. The system is designed to provide world-wide 
coverage with accuracies better than one nautical mile. 

The position of the vehicle is determined by measuring the 
relative phase of the signals transmitted from two or more 
stations . 

OMEGA has the following advantages: 

1. ) Range 8000 nm. 

2. ) Coverage is global with just eight stations. 

3. ) Continuous fixing available. 

4 . ) Passive receiver on the vehicle. 

5. ) All weather capability. 

The disadvantages of OMEGA are: 

1. ) Position fix is in hyperbolic coordinates. 

2. ) Dependent on fixed ground stations. 

3. ) Susceptible to Ionospheric disturbances. 

4. ) Susceptible to jamming. 

2 . The Satellite navigation System 

Satellite navigation is a method of fixing the 
vehicle's position from data obtained from an artificial 
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earth satellite. The position of an artificial satellite 
can be predicted by the orbit it is traveling. A relative 
fix can be made from the satellite by measuring the elevation 
angle, the azimuth angle, range to the satellite, or the 
rate of change of altitude, azimuth or range. If a single 
variable is measured, successive measurements are required. 

There are three basic types of satellite navigation: 
Doppler Systems, which measure rate of change of slant range; 
angle measuring systems, which measure elevation angle; and 
ranging systems, which measure the slant range. 

Satellite Navigation Systems have the following 
advantages : 

1. ) All-weather position fixing. 

2. ) Communication service to the vehicle while in flight. 

3. ) Worldwide service. 

*1.) Capable of accuracies up to 200 ft. 

5. ) Can provide vehicle with velocity and heading 

information as well as position information. 

6. ) Passive at the receiver. 

The disadvantages of a Satellite Navigation System 
are as follows: 

1. ) Special receiving equipment required including 

antennas with vertical patterns. 

2. ) The satellite orbit must be predicted and 

continuously up-dated. 

3. ) The satellite's signal may be susceptible to 

jamming. 

4. ) The satellite signal is susceptible to Ionospheric 

refraction . 



57 



5.) Navigation errors are introduced by drifts in the 
frequencies of the satellite transmitters. 

3. The Terrain Mapping System 

A Terrain Mapping System operates on the principle 
that the geographic location of any place on the land surface 
of the earth is uniquely defined by the vertical contours 
or. topography of the surrounding area. The Terrain Mapping 
System measured the vertical contour of the terrain along 
its flight path, using a radar altimeter to measure clearance 
above the terrain and a reference altitude of the vehicle. 

By subtracting instantaneous radar measured altitude from 
the reference altitude, the Terrain Mapping System determines 
the terrain contour. The system then searches its computer 
memory to find a stored terrain contour, whose coordinates 
are known, which closely matches the measured one. This 
serves to fix the vehicle's position. 

The advantages of the Terrain Mapping System are 
as follows: 

l 

1. ) It is an all-weather system. 

2. ) It is completely self-contained. 

3. ) It is inherently resistant to jamming. 

H.) It is accurate up to 50 ft. 

5. ) Accuracy demands imposed on the radar altimeter are 

not excessive. 

6. ) It is unaffected by man-made changes in topography. 
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The Terrain Mapping System has the following 
disadvantages : 

1. ) The terrain contour to be flown over must be 

previously mapped. 

2. ) To increase the accuracy of the system, there must 

be an increase in the computer memory capacity. 

3. ) A dead-reckoning system must also be carried so the 

Terrain Mapping System knows where to search for a 
contour match. 

l\.) It can not be used to obtain position fixes over 
the ocean or large bodies of water. 



F. THE INTEGRATED SYSTEM 

The basic parameters of the navigation systems currently 

available have been reviewed. In order to select the proper 

long range navigation system for military use, a minimum 

set of system requirements was established: 

Range 3000 nm. 

Accuracy 5 nm. 

Velocity 0 - 2000 kts. 

Coverage Worldwide 

Weather Restriction All-weather capability 

Flexibility Capable of self-testing and 

switching mode when needed 
Equipment Non-j ammable , completely 

self-contained 



A review of the previously mentioned navigation systems 
indicates that no single system available is capable or 
flexible enough to meet the system requirements. The Doppler 
Navigation System is self-contained and provides very accu- 
rate average velocity measurements, however, the system can 
not sense short term velocity fluctuations and must depend 
on an external heading source. The Inertial Navigation 
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System is self-contained and provides excellent short- 
term position, velocity, and heading information. The 
Inertial Navigation System, however, has a long-term error 
build-up due to the inherent Schuler Cycle in the system. 

In order to meet the system requirements, an integrated- 
navigation-system was chosen. An integrated-navigation- 
system feeds the outputs of several navigation systems into 
a central computer which then provides a single more accu- 
rate output. The central computer provides decision making 
ability to test the accuracy of the individual navigation 
systems. By inputting fixed positions to the central com- 
puter at regular intervals, the computer can determine 
system drift parameters and provide the proper system bias 
velocity. 

The navigation system developed in this paper is an 
Inertial/Doppler system integrated by an MCS-lJ microcomputer. 
The microcomputer combines the short-term accuracy of the 
Inertial with the long-term accuracy of the Doppler to 
obtain the most probable position of the vehicle. The 
micro-computer is programmed to provide decision making 
flexibility to ensure the outputs remain accurate during 
partial system failure. The wind influencing the vehicle 
is continually computed and updated in the MCS-JJ memory so 
that the computer can automatically switch to an air-data 
mode of operation in case of Inertial and Doppler failure. 
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The navigation equations used to compute the vehicle’s 
position from the inputted data are covered in the next 
section. 
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IV. NAVIGATION EQUATIONS 



The navigation equations necessary to calculate the 
position of the vehicle relative to the earth are presented 
in this section. The mathematical model presented takes 
the dead-reckoning outputs of the Inertial, Doppler, and 
Air-Mass Systems and extrapolates the present position of 
the vehicle from the last known position. Direct position 
data can be inputted into the algorithm and is used to update 
the position of the vehicle . The system drift error is 
computed by comparing the known position with the dead- 
reckoned computed position. 



In describing the equations in detail, a graphical pic- 
ture of each equation's function is included. The graphical 
picture represents each equation as a black box with speci- 
fied inputs which produce the desired outputs. A complete 
description of graph theory as applied to this paper is 
given in Section V-A. 



A. THE EARTH MODEL AND THE COORDINATE REFERENCE FRAME 

The first step in developing the Navigation equations 
is to prepare a mathematical model of the earth. The earth 
can be approximated as a sphere of radius R, the nominal 
equatorial radius. However, this approximation would produce 
an unacceptable error as great as 3 miles for every 1,000 
nm. traveled. This is due to the flattening of the earth 
at the poles and the bulging of the earth at the equator. 
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A more exact model of the earth is the reference 



spheroid where the earth is assumed to be an ellipsoid 
of revolution with the semi-major axis R, the nominal 
equatorial radius, and the semi -minor axis of radius P, 
the polar radius. Accuracy to within 30 ft. for every 
1,000 nm. traveled is possible with this model. The 
determination of angels, however, in this system is very 
difficult. 

The model chosen in this system is a combination of the 
two previously stated models. All angles are determined 
from the spherical model thus allowing the convenient lati- 
tude/longitude coordinate reference frame to be used to 
map each point on the real earth to the earth model. The 
reference spheroid is then used to determine arc lengths 
that more closely represent the actual shape of the earth. 
The arc lengths are applied to the determined latitude and 
longitude in the form of conversion constants. This method 
of computation is very simple and has a maximum error of 
300 ft. for every 1000 nm. traveled. 

B. DERIVATION OP CONVERSION CONSTANTS 

The international nautical mile equals 6076.1033 it. 

One minute of arc length measured at the equator equals 
6087.08 ft. The longitude conversion constant, KLO, is used 
to relate nautical measurement to the actual geographical 
measurement : 
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KLO = 


feet/degree 
feet/nautical mile 


KLO = 


6087.08 feet per minute/60 minutes per degree 
6076.1033 feet/1 nautical mile 


KLO = 


60.108^1 nautical miles/degree 



The average length of a minute of arc length on a 
meridian Is 6076.82 ft. To determine the latitude conver- 
sion constant, KLA, which is also used to relate nautical 
measurement to the actual geographical measurement, the 
following equation is used: 



KLA = 


feet/ degree 
feet/nautical mile 


KLA = 


6076.82 feet per minute/60 minutes per degree 
6076.1033 feet/1 nautical mile 


KLA = 


60.0071 nautical miles/degree 



The distance traveled in nautical miles, east-west and 
north-south, is divided by the conversion constants resulting 
in distance traveled measurements in degrees, longitude and 
latitude . 

C. WIND CALCULATIONS 



The wind effecting the motion of the vehicle can be 
obtained by comparing the north-south and east-west vector 



components of the true airspeed, as computed by the air- 
mass system, with the vector components of ground speed, as 
computed by the inertial system. The calculated wind com- 
ponents must be continually updated so the current wind 
components stored in the computers memory are the most cur- 
rent. In case of failure of the Inertial and Doppler Sys- 
tems, the computer will be capable of computing the dead- 
reckoning position from the outputs of the air-mass system 
and the last stored wind components. 

The wind velocity components are computed from the 
following equations: 




■v 

*■ 



VWXR = VGXI - VAX 
VWYR = VGYI - VAY 

The wind velocity components computed are subject to 
gusts, maneuvers, and other short term disturbances. In 
order to decrease the effects of random fluctuations in the 
wind velocity, the new computed wind velocity is averaged 
with the previously computed wind velocity in a smoothing 
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routine. The smooth wind velocity is derived from the 
following equations: 



-VWXR 



•VWYR 



•VWX. 



i-1 



■VWY. 



i-1 



-v - 1 



SMOOTH 

WIND 

VELOCITY 

COMPUTATION 



•wx.- 



-VWY .- 
1 









VWX = G VWX + ( 1-G) VWXR 
VWY = G VWY + (1-G) VWYR 



D. DEAD-RECKONING CALCULATION 

Dead-reckoning is a means of navigation in the absence 
of position fixes. The vehicle's position can be estimated 
by measurements of the groundspeed components. The ground- 
speed components are integrated over a given time interval 
to give the distance traveled in that time interval. The 
estimated position of the vehicle is the summation of the 
distance increments traveled from the last known fixed 
position . 

1 . Air Mass Mode 

The air mass mode of operation is a back up for the 
two primary dead-reckoning navigation systems, the Inert jal 
and Doppler. The components of airspeed, true heading, and 
stored wind information are combined to output distance 
increments in the north-south and east-west directions. 
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The true air-speed inputs are influenced by random 
fluctuations in the Pitot-Static System. In order not to 
influence the navigation cycle by these fluctuations, the 
inputted true air-speed is averaged with the previously 
stored true air-speed. The smooth true air-speed is derived 
from the following equation: 



TASR 



TAS 



i-1 



SMOOTH TRUE 

AIR-SPEED 

CALCULATION 



TAS£ > 



TAS 



TAS + TASR 
2 



The north-south and east-west components of air-speed 
are calculated using the previously calculated true-air-speed 
and the true heading inputted from the INS. In case of INS 
failure, magnetic heading may be inputted from the flux 
valve together with the magnetic variations applicable set 
in by hand. The air-speed components are calculated as 
follows : 



TH- 

TAS 



COMPONENT 

AIR-SPEED 

CALCULATION 



VAX > 

VAY > 



VAX = TAS SIN (TH) 
VAY = TAS COS (TH) 
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The air-speed components represent the velocity and 
direction of the vehicle through the air— mass. The ground 
track of the vehicle can be found by adding vectorally the 
motion of the air-mass, wind velocity, to the air-speed as 
shown in Figure 3. The equations for the ground velocity 
in the air-mass mode are as follows : 



VGXW v 

VGYW »- 



VGXW = VWY + VAX 
VGYW = VWY + VAY 



VAX 
VAY • 
■ VWX ■ 
VWY- 



GROUND 
VELOCITY 
IN AIR-MASS 
MODE 



Lastly, the ground distance traveled by the vehicle 
during one navigation cycle of the computer can be determined 
by integrating the ground velocity over the time interval of 
the navigation cycle. 

VGXW *- 

VGYW >■ 

T v 



DISTANCE 
CALCULATION 
AIR-MASS MODE 



DXW - 
DYW- 
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2 . Doppler-Mode 



The Doppler System obtains very accurate average 
velocity measurements and converts them to distance measure- 
ments along and across the vehicle’s true heading. These 
measurements along with true heading, from an external 
source, can be used to calculate the distance traveled by 
the vehicle as shown in Figure 3* 

The dead reckoning position increments computed from 
the doppler are derived from the following equations: 



SDD > 

SHD *■ 

TH *- 

DXD = SDD x COS (TH) + SHD x (TH) 

DYD = SHD x COS (TH) - SDD x SIN (TH) 

3. Inertial Mode 

The Inertial Navigation System, (INS), provides 
excellent short-term velocity and heading information. The 
outputs of the INS are the velocity components of the vehicle 
north-south and east-west. The ground track of the vehicle 
can be determined as shown in Figure 3- 

The dead reckoning position increments computed 
from the Inertial are derived from the following equations: 



Doppler 



Inertial 



Doppler 

Distance 

Computation 



■DXD 

•DYD 
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VGXI 

VGYI 



DXI 

DYI 



+- 



Inertial 



■> 

> 



Inertial 

Distance 

Computation 



DXI = T • VGXI 
DYI = T • VGYI 



4 . Integrated System Mode 

To take advantage of the extremely accurate average- 
distance measurement of the Doppler and the precise short- 
term distance measurement of the Inertial, the distances 
measurement by each is averaged in a smoothing routine to 
give the most accurate dead-reckoning distance measurements 
possible. The smoothing factor, G, used in the smoothing 
equation is set depending on the accuracy of each system. 



DXI 

DYI 

DXD 

DYD 



Distance 

Smoothing 

Routine 



DY * 

DY * 



DX = G x DXI + ( 1-G ) • DXD 



DY = G x DYI + (1-G) • DYD 
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E. POSITION CALCULATION 

1. Geographic Fix Position 

Any time the vehicle's position is fixed, the updated 
fixed position will be loaded into the computer's memory 
along with the time the fix was taken. The dead-reckoning 
position coordinates at the time of the fix will be up-dated 
to the new position coordinates and the computer will continue 
to dead-reckon the vehicle's position from the last inputted 
fix. 



LOF 

LAF 



Update 

Nav 

Position 



LON 

LAT 



•v 






LON = LOF 



LAT = LAF 



2 . Up-Date Geographic Position 

The initial position of the vehicle is loaded into 
the computer memory before launch. The geographic position 
is up-dated each time a fix is taken. The geographic posi- 
tion between fixes is updated by dead-reckoning distance 
computations from the Doppler and Inertial. The zero sub- 
script indicates the vehicle's previous geographic position. 
The computed system drift increments are also inputted to 
account for system drift errors in the dead-reckoning 
computation . 
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LON >. 

LAT *- 



LAT 



LAT + 
o 



DY+SDY 

KLA 



LON 



rnM DX + SDX 

^ o KLO x COSC LAT) 



3. System Drift Error 

Each time a new position fix is taken, the fixed 
position is compared with the calculated dead-reckoned 
position. From this information, the amount of drift in 
the dead-reckoning computations between fixes is determined. 



LOF 

LAF 

LON 

LAT 



System 

Drift 



SDXC 

SDYC 



SDXC = (LOF-LON) x KLO x COS (LAT)' 

SDYC = (LAF-LAT) x KLA 

^ . System Bias Velocity 

The system drift error computed from the previous 
fix is used to update the system bias velocity. The system 
bias velocity is used to cancel the system drift error. 
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SDXC 

SDYC 



System 

Bias 

Velocity 



VSDX 

VSDY > 



VSDX = VSDX + (SDXC x TF) 



VSDY = VSDY + (SDYC x TF) 



5 . System Bias Distance 

The system bias distance is computed by integrating 
the system bias velocity over the time of one nav cycle. 

The system bias distance is applied to the up-date geographic 
position to cancel the effects of the system drift. 



VSDX 

VSDY 



- — v 


System 






Bias 






Distance 





SDX 

SDY 



-> 

■*- 



SDX = VSDX x T 
SDY = VSDY x T 



The navigation program combines these equations as 
described in the Navigation Functional Flow Chart, Figure 5 - 
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V. PROGRAMMING THE MICROCOMPUTER FOR NAVIGATION 



The Navigation equations discussed in Section IV des- 
cribed how to process the outputs of the Inertial, Doppler, 
and Air-Mass Systems in order to calculate the current 
position of the vehicle. In order to make a working system, 
the Navigation equations had to be programmed into the 
software of the MCS-4 microcomputer. 

The program was designed to incorporate the following 
principles and requirements: 

1. ) The total program had to be written in the assembly 

language of the MCS-4 microcomputer. 

2. ) The total run time of one navigation cycle could 

not exceed 200 milliseconds. 

3. ) The design of the program should be modular to 

facilitate the addition or deletion of new code. 

4. ) The program should consist of an executive routine 

that calls upon various subroutines as indicated 

on a process graph of the navigation program. 

5. ) The design allows for substitution or addition 

of n ew functions, routines, and/or programs in 

a straight forward manner. 

6. ) For each built-in function, a single subroutine 

is called to perform the calculation. 

The instruction set of the MCS-JJ microcomputer, as des- 
cribed in Section II. C., contains only addition and subtrac- 
tion of two four digit binary numbers as the fundamental 
arithmetic operation. In order to make the detailed pro- 
gramming easier, a so-called "process graph" method was 
used to develop the program. 
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A. GRAPH THEORY 

Programming the navigation equations would be an easy 
task in a higher level language such as FORTRAN. To program 
the same equations in an assembly language requires a detailed 
knowledge of the current contents of all index registers 
and memory locations. The status of even one single bit 
such as the carry bit can net be left unaccounted. The 
requirements of keeping minute details in mind led to the 
development of a graphical means of representing the 
microcomputer program. 

Graph theory provides a simple and powerful tool for 
constructing mathematical models of discrete arrangements 
of objects. The process graph consists of vertices which 
are pairwise connected by a directed line. 

The best means to illustrate the concept of graph theory 
applied to microprogramming is through an example. Consider 
the problem of adding two numbers together which reside at 
given locations in memory and whose sum is to be located 
in a specified location in memory depending on its sign. 

This process can be described by the following flowchart: 




This notation does not give much information about how the 
processor is to accomplish this process, how much memory is 



required by the program, the input and output variables, 
and how much time is needed to carry out the process. 

The same example written as a process graph could be 
as follows: 




The process graph indicates that the variables X and Y are 
inputs to the function add and Z is the output of that func- 
tion. A test is made of the output, Z, to decide where Z 
will be located. The amount of time and memory involved in 
this computation can be easily established by referencing 
the known memory space and time required for each operation 
involved. The depth to which the programmer draws the re- 
quired graph for his computation is a function of the capa- 
bility of the operations he has available to him together 
with his own preferences. The same example discussed could 
be graphed in more depth, as follows : 



X 




20 



Y 



HAM 30 



Another programmer reading this graph only needs to know 
the following in order to write the required program: 



X 





— ^ 


ACC 







Load the variable X 
into the accumulator 



ACC 
Y - 






-ACC 



Add the variable Y to 
the contents of the 
accumulator 



ACC 




Test the sign of the 
contents of the 
accumulator 



+ 


> 


RAM 


ACC 





If the sign of the 
accumulator is positive , 
store the contents of 
the accumulator into 
RAM location (2,0) 







RAM 


ACC 


> 





30 



If the sign of the 
accumulator is negative, 
store the contents of 
the accumulator into 
RAM location (3,0). 



The detailed process graphs used in developing the navi- 
gation program were a great aid in the design, analysis and 
documentation of the program. The following is a list of 
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the problems that could be analysed directly from the 
process graph: 

1. ) What portions of the problem are sequential. 

2. ) What portions of the problem may be processed in 

parallel by independently operating processors. 

3. ) Where and how intercommunication may take place 

when a multiplicity of processors is used. 

4. ) What capacity must the processor possess before it 

can be successfully used to solve the problem. 

5. ) What is the maximum allowable time for the operation 

of each subprogram. 

6. ) What is the memory space requirement for each 

subroutine . 

B. DEVELOPMENT OF NAVIGATION PROCESS GRAPHS 

The MCS-JJ microcomputer was first thought of as a black 
box taking inputs from the Inertial, Doppler, Air-Mass, and 
Position Fixing Systems and outputting the vehicles current 
position in latitude and longitude. The graphical represen- 
tation of this system is shown in Figure 4. 

Figure 4 represents the total navigation system. The 
Inertial, Doppler, Air-Mass, and Position Fixing Systems 
provide input to the system as described in Section IV. 

The microcomputer program was required to complete the sys- 
tem. The programming was accomplished by continually break- 
ing down the functional operations into smaller and smaller 
parts until the operations were simple enough to be directly 
written in the MCS-ll machine language. The functional proces 
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NAVIGATION SYSTEM BLOCK DIAGRAM 



graph. Figure 5 , represents how the microcomputer program 
was initially broken down into the fundamental navigation 
equations as described in Section V. 

The functional process graph was analyzed to define the 
flow of the input variables and to determine the feasibility 
of a multi-processor system in order to shorten the required 
computational time. It was noted that the program could be 
broken into two parts, the calculation of the distance in- 
crements, and the calculation of the latitude and longitude 
from the distance increments. The analysis demonstrated 
that the computation time could be nearly halved by having 
two microcomputers working simultaneously to produce the 
desired outputs. The system was designed to have one micro- 
computer receive the given inputs and compute the distance 
increments traveled, while at the same time, the second 
microcomputer computes the latitude and longitude from the 
previously calculated distance increments. 

The functional process graph was then broken down into 
the operational process graphs for each microcomputer, 

Figure 6, and Figure 7. The operational process graphs 
represented the desired program for each microcomputer. 

The program described by the operational process graph, 
Figure 6, was written to investigate the time and effort 
required to develop the required software. 

The program was written in a modular form as required. 
Each operation represented in the operational process graph. 
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OPERATIONAL PROCESS GRAPH 
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NAVIGATION POSITION PROCESS GRAPH 



was written as a separate subroutine. An executive routine 
was then developed to call each subroutine in the proper 
order and account for each variable as required by the 
operational process graph. 



C. PROGRAM ANALYSIS 

A program analysis was developed to define those pro- 
blem areas that had to be solved before programming the 
microcomputer. The areas that were investigated were com- 
putational speed, memory space available, and accuracy 
required. 

The operational process graph. Figure 6, represented 
the tasks to be accomplished. The type of operations and 
number of operations that wore required are listed as follows: 



Operations 



Times Called 



Multiply 12 
Cosine 1 
Sine 1 
Addition 12 
Subtraction 3 
Division by two 7 

Total 36 



The operational process graph was used to determine the 
speed limitations on each operation programmed. The total 
navigation cycle was limited to 200 milliseconds. Since the 
number of operations required to be performed differed depen- 
ding on what sensors were operational, the critical path of 
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the operational process graph was determined. The critical 
path occurred when the Inertial and Doppler were both opera- 
tional, Figure 8. The operations involved in the critical 
path were ten multiplies. Cosine and Sine calculations, ten 
additions, three subtractions, and seven divisions by two. 
The total computational time for the critical path was 
limited to 200 . milliseconds . 

The critical operations which had to be developed were 
the Multiply, Cosine, and Sine subroutines. Investigating 
previous work on the MCS-4 indicated that previously pro- 
grammed multiply. Cosine, and Sine routines were requiring 
50 msec., 650 msec., and 750 msec, respectively. In order 
to program the microcomputer for navigation, the development 
of the folloiving routines were required with the following 
constraints : 



The time used by each of these subroutines determined the 
number of processors required to meet the 200 millisecond 
time constraint. 

The amount of memory available to program the navigation 
routine was a function of the h 0 0 ^4 CPU and the number of 
microprocessors used. One ^00^1 CPU can directly drive six- 
teen ROMs and sixteen RAMs. The number of instructions in 



Multiply 

Cosine 

Sine 



less than 200 msec, 
less than 200 msec, 
less than 200 msec. 
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Figure 8. OPERATIONAL CRITICAL PATH 



the navigation routine was limited by the space available 
in the ROMs. Since each ROM could hold 256 instructions, 
the program was limited to 4,096 instructions per micro- 
processor. 

The time required to sequentially execute every instruc- 
tion in the sixteen ROMs would be 44 mlliseconds. Since the 
program was required to be written with an executive routine 
and a set of subroutines that would repeat the same set of 
instructions several times, the limiting time constraint 
would be reached before using up the available ROM space 
in one microprocessor. It was determined from this analysis 
that ROM space would not be a limiting factor in writing 
the navigation program. 

The amount of memory space available to store the values 
of each variable was determined by the space available in 
the RAMs. Since one 4004 CPU could drive 5120 bits of RAM, 
the navigation program was limited to 320 variables of 16 
bits for each microprocessor. The number of variables re- 
quired was determined from the operational process graph, 
Figure 6, where each line connecting a pair of vertices 
represents one variable. There were 54 lines indicating 
that a maximum of 54 variables were required plus those 
variables used in any single operation. By overlaying 
variables in the same RAM memory space, the memory space 
requirement was reduced. It was determined from this analysis 
that RAM space available would not be a limiting factor. 
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The required accuracy of the navigation program was a 
function of the accuracy of the input variables. The accu- 
racies of Inertial, Doppler, and Air-Mass systems used on- 
board the P3C aircraft were used in this analysis as repre- 
senting the state-of-the-art systems in naval aircraft today. 
The specifications for these systems are as follows: 



System 

Inertial 
True Heading 
Doppler 

True Air Speed 



Designation 

ASN-8^ 

ASN-8^ 

APN-187 

Pitot-Static 



Accuracy 

± 1.5 knots RMS 
± 9 ARC -MIN RMS 
± 1.0 knots RMS 
± 2.0 knots RMS 



The accuracy of the navigation program was a function 
of the accuracy of input data as well as the bit size 
assigned to each variable. The limited accuracy of the 
input data permitted each variable to be no greater than 
16 bits. This allowed each variable to be represented by 
four hexadecimal-digits with the first bit assigned as the 
sign bit. The hexadecimal point for speed measurements was 
fixed so that there is one hexadecimal digit to the right 
of the decimal point. This allowed the accuracy of the 
speed inputs to be within ± .0625 knots. The range on the 
inputs due to a 16 bit variable limitation was ± 20^17.99 
knots. The accuracy requirement was not considered a major 
limitation in the program analysis. 

The program analysis pointed out the major areas that 
had to be demonstrated feasible before programming of the 
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microcomputer for navigation was initiated. The program 
analysis was used to set the design goals at each step and 
the process graph proved to be a major tool in performing 
the program analysis. 



D. SUBROUTINES 

A large number of the mathematical operations required 
in the navigation program are repeated many times. In order 
to decrease the total programming effort and also decrease 
the memory-capacity requirements, many of the operations 
required were written as subroutines. These routines con- 
sist of a series of instructions dedicated to a specific 
task. The subroutines developed for the navigation program 
were divided into two groups, these involving complex mathe- 
matical operations and those involving more common func- 
tional operations. Each subroutine was written in a general 
form to permit its recurring use by the executive program. 

1 . The 'Multiplication Routine 

The major limitations of the MCS-^ microcomputer to 
be overcome were the limited instruction set and slow speed 
of calculation. The most powerful arithmetic instruction 
available was the single addition two hex-digits. The mul- 
tiplication routine was written to increase the capability 
of the MCS-lJ microcomputer in order to satisfy the require- 
ments of the navigation program. 

A multiplication routine had been written involving 
multiplication by a series of additions. The program required 
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only fifty instructions, however, the computational time 
was 40 milliseconds. Since the navigation program required 
a minimum of ten multiplications, this method was unsatisfactory. 

One of the advantages the MCS-4 microcomputer has 
is its inexpensive memory. It was decided to investigate 
a different w ay of programming the microcomputer that would 
take advantage of available memory. It was discovered that 
memory space could be traded for speed by a table look-up. 

The multiplication subroutine was written using a table 
look-up scheme. 

a. Capability of Subroutine 

The multiplication subroutine was designed to 
take a four hex-digit number, X, and multiply it by a four 
hex-digit number, Y, resulting in a four hex-digit chopped 
number. The input X is loaded in Index Registers R8 -*■ RB 
with the least significant digit, Xo, loaded in R8. The 
Input Y is loaded in Index Registers RC -»■ RF with the least 
significant digit, Yo, loaded in RC. The product will appear 
in Index Registers RC -*• RF with the least significant digit 
located in RC. 

b. Description of the Routine 

The process was done by a table look-up scheme 
similar to longhand multiplication. An example of this 
method follows: 
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X 3 Y 2 


X 2 Y 2 


X 1 Y 2 


X 0 Y 2 






x 2 y 3 


X 1 Y 3 


X 0 Y 3 









The table of values used in the multiplication routine con- 
sits of a 16 x 16 matrix of product values. Each product 
value is an exact value of a multiplication of two single 
hex-digit numbers. The row that would normally contain the 
products of a zero multiplication is used for instructions 
within the ROM containing the table. A test for a zero input 
is made in the body of the subroutine. Each product value 
within the table consists of a two hex-digit number with 
the second value being the least significant digit. During 
the execution of the program, each digit in the two hex- 
digit number is loaded in a separate index register. In 
this write up, an "M" stands for the most significant digit, 
an "L" stands for the least significant digit, and a "CY" 
stands for a carry. The longhand multiplication with the 
portion chopped in this program can be shown as follows: 
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A numerical example of this procedure in decimal is as follows 

9 7 5 3 



7 

54 4 _ 

36 28 2 _ 

18 i4 10 0 - 

6 10 13 

2 4 0 5 



c. Development of the Program 

The method used in this procedure becomes very 
complex because of the large number of separate hex-digits 
involved and the small number of index registers available 
to store each digit. This problem becomes more complex 
since each addition of two hex-digits creates a possible 
carry. A solution to this problem is to make a process graph 
that simulates the multiplication process. 
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The process graph for the multiplication routine 
is shown in Figure 9. The symbols used have the following 
meaning: 




table value of product of the 
digits X Q times 



addition 



The following examples is used to show the flow 
of the hex-digits throughout the routine. 



Input: X = 1 2 3 ^ Hex 

Y=5678 Hex 

Calculation by hand: 

5678 
1 2 3 'I 

0 8 
07 0 E 

06 OC 12 

05 0A OF 1 ^ 

0 1 2 2 py 

CY CY CY ox 



0 6 2 ^ 



Figure 10 shows the complete flow of the numbers throughout 
the routine. The value obtained agrees with that computed 
by hand and can be used as a check of the process graph. 
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PROCESS GRAPH FOR MULTIPLICATION ROUTINE 
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SAMPLE MULTIPLY USING PROCESS GRAPH 









d. Constructing the Program 

The process graph was used to assign and keep 
track of the Index Registers being used (Fig. 11). 

The subroutine first stores all the values in 
the Index Registers into RAM, R2 and R3 are used to address 
values in the RAM. RO and R1 are used to fetch values from 
the multiplication table. R*l RF are the working registers 
and are assigned as shown in Figure 11. 

e. Expansion of Multiplication Routine 

The multiplication routine was written to give 
a truncated product of two positive numbers accurate to 
four significant digits. The multiplication routine which 
handles both positive and negative numbers is shown in the 
process graph. Figure 12. 

Each number inputted into the multiplication is 
tested for its sign. The sign bit is then shifted out to 
normalize the number to four significant digits. The pre- 
viously described multiplication routine is then executed. 
The resulting number is shifted back to its proper form 
and the required sign bit is set. 

f. Summary 

The multiplication routine developed required 
the memory space of three ROMs and computed its result in 
five milliseconds. A complete listing of the multiplication 
routine is found in Appendix C. 
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MULTIPLY INDEX REGISTER ASSIGNMENT 
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2. The Cosine Routine 



The development of a Cosine routine for the micro- 
computer which could compute sufficiently fast was one of 
the major programming tasks. Cosine Routines written for 
general purpose computers are usually written as series 
approximations in order to save memory space. The object 
of the routine written for this program was to increase 
the speed of calculation. 

Two Cosine routines previously programmed on the 
MCS-4 Microcomputer were investigated. The first was a 
Chebyshev approximation routine which required 750 milli- 
seconds to compute the Cosine. The second routine investi- 
gated was a Cordi c approximation which required 350 milli- 
seconds. Both methods were too time consuming for this 
project . 

The procedure developed in this project was a table 
look-up, linear interpolation routine. The Newton Divided- 
Difference Interpolating Polynomial was used because of its 
simplicity. The size of the table required and the accuracy 
of the results are both functions of the degree of the 
Polynomial used. For simplicity and speed, a first-order 
divided-difference table was used which resulted in a linear 
interpolation of the form: 

F(X) = F(X q ) + (X - X Q ) F[X 1 ,X q ] 
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where 



X = 6 



F(X) = COS 0 



F[Xi,X c ] - 



F(X 1 ) - F(X q ) 
X 1 " X 0 



The table consisted of all values of F(X i ) and 
F[X.,, X.] for F(X) = COS 0, 0 <_ 0 <_ 1.88 radians In hexa- 
decimal. The increment used for each value in the table 
was .08 radians hexadecimal. The table was constructed 
from a Fortran program which used a decimal increment of 
q , equivalent to .08 hexadecimal, and outputted the 
desire table values in hexadecimal, 
a. Table Length 

The size of the table loaded into the program 
was a function of the required accuracy of the Cosine rou- 
tine. The data supplied to the program from the navigation 
devices was accurate to three significant figures. 

Different size tables were constructed and 
tested for accuracy. Since the interpolation was linear, 
the largest error occurred at the mid-point between each 
table value. 
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Maximum Error 



Linear Approx. 

The table size was adjusted until the maximum 
error was within three units in the fourth significant 
figure, thus guaranteeing three significant figure accuracy. 

The table values were loaded sequentially into 
ROM number three. The first entry was the value of F(X n ) 



The remaining part of the ROM was used for the interpolating 
routine . 



from outside of the ROM. The input registers RC, RD, RE, 
RF with RF the least significant figure. 




b. Description of the Program 



The Cosine routine was designed to be called 



Examples : 



RC RD RE RF 



6 = 1.520 
0 = 0.7^8 



would be loaded 0 2 5 1 

would be loaded 8^70 
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The description of the programs function is best shown by 



an example. Let the 


input be 


CD 

II 

O 

XT 


Therefore 


RD 


RD 


RE 


RF 


0000 


0100 


0111 


0000 


0 




7 


0 



The first object of the program was to find the 
required values (F(X^) ,F(X^ + ^ } X^) ) in the table. The ROM 
can be thought of as a 16 by 16 matrix of 8 bit words. Each 
byte therefore will hold two significant figures of each 
table value. Since each table value has four significant 
figures and two table values were required for each compu- 
tation . only four bytes needed to be retrieved from the 
table. If the four required bytes for each computation are 
thought of as a unit, the table could be thought of as rows 
of units with four units in each row. The increment in 
radian value therefore between each row is .20 radians hexa- 
decimal. The required row in the table is easily found by 
dividing the inputted angle by two. 

9JA9. = 03 AO (hexadecimal) 

The required row in this example is row three. To find the 
required unit in row three, the remainder of the inputted 
angle is divided by four. 
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2 °^. = 0028 



The required unit in this example is unit two. The final 
remainder of the inputted angle, 0008, is equal to the 
difference (X - X^) . 

In summation, the location of F[X i+1 ,X jL ] in 
the table was found by simple divisions by two and four 
performed on the input and the difference. Divisions by 
powers of two can be obtained by shifting operations. 

After the desired table values and difference 
were found, the multiplication routine was called which 
multiplied (X - X ± ) times FCX^^] and stored the result 
in index registers RC to RF. The final computation tooK 
the value in RC to RF and subtracted it from F(X^) and stored 
the result in index register RC , RD, RE, and RF to be re- 
turned to the executive program. The step by step procedure 
of the Cosine routine is shown in the Cosine process graph. 
Figure 13. 

c. Expansion of Routine 

The table used in the Cosine Routine was also 
the table required by a Sine routine. It was noted that 
this routine could also be used to find the sin 0 for 
0 < 6 < 90 by subtracting the input G from 90 (1.88 rads-0 
rads) and using the same routine since cos(90 - 0) = sin(O). 
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d . Summary 

The memory space required by the Cosine routine 
was one ROM plus the space taken up by the subroutines called 
by the Cosine routine. The main part of the Cosine routine 
contains only 46 instructions. The time required to execute 
the Cosine routine was basically the time required to execute 
the multiplication. The Cosine routine required only a 
total of 5.17 milliseconds. This computational speed repre- 
sents a 70 fold decrease in the computational time to com- 
pute the Cosine by previously available routine. By table 
look-up schemes, it was proven that the computational speed 
of the microcomputer could be competitive with that of a 
general nurnose computer. The feasibility of using the 
MCS-4 as a navigation computer was also proven by demonstra- 
ting its ability to find the Cosine and Sine at competitive 
speeds and within the memory limits and accuracy limits 
required. 

3 . Common Routines 

The common subroutines were written to do the basic 
housekeeping operations such as storing data, simple arith- 
metic, shift operations, and transfer of data between RAM 
and IR. These subroutines were called by the executive 
routine and the multiply and Cosines routines to aid in the 
data handling. The functions handled by the common subrou- 
tines were those best suited for the MCS-4 and therefore 
could be written in a straight forward way requiring little 



speed or memory space. The functions of the common routines 
were broken into three groups: Arithmetic, Shifting, and 

Data Handling. 

a. Arithmetic 

The arithmetic routines handled the simple addi- 
tions and subtractions required in the navigation program. 
These routines were handled well in the MCS-4 by the 4-bit 
ripple-through carry type adder incorporated in the 4004 
CPU. This allowed direct addition or subtraction of two 
hexdigits in either the accumulator or RAM. There were 
two addition routines, two subtraction routines, and two 
special purpose arithmetic routines written. 

The two addition routines written were ADDRAM 
and ADDRAMIR. ADDRAM adds the contents of the RAM addressed 
by IRO and IR1 to the contents of the RAM addressed by IR2 
and IR3 and stores the results in the RAM addressed by IR2 
and IR3. The contents of the RAM addressed by IRO and IR1 
remains unchanged. 



RAM (IRO, IR1 ) > 

RAM (IR2, IR3) ■> 



->RAM ( IR2 , IR3) 



ADDRAM occupies 12 words of ROM and takes 288.8 microseconds 
of computation time. 

ADDRAMIR adds the contents of the RAM addressed 
by IRO and IR1 to the contents of IRC thru IRP and stores 
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the results in IRC thru IRF. The contents of the RAM 
addressed by IRQ and IR1 remains unchanged. 



RAM (IRO, IR1 ) > 

IR(C F) j 



■> IR ( C + F) 



ADDRAMIR occupies 22 words of ROM and takes 237-6 micro- 
seconds of computation time. 

The two subtraction routines written for the 
navigation program are SUBIR and SUBRAMIR. 

SUBIR subtracts the contents of IRC thru IRF 
from the contents of IR4 thru IR7 and stores the result in 
IRC thru IRF. The contents of IR^I thru IR7 remains unchanged. 



IR(H -v 7) > 

IR(C + F) > 



-» IR( C -v F) 



SUBIR occupies 17 words of ROM and takes 182.6 
microseconds of computational time. 

SUBRAMIR subtracts the contents of IRC thru IRF 
from the contents of the RAM addressed by IR2 and IR3 and 
stores the result In the RAM addressed by IR2 and IR3- The 
contents of IRC thru IRF remains unchanged. 



RAM( IR2 , IR3) •> 

IR(C - F) v 



->RAM( IR2 , IR3) 
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SUBRAMIR occupies 32 bytes of ROM and takes 



3*15.6 microseconds of computational time. 

The two special arithmetic routines written for 
the navigation program are COMPLEMENT and COMANGLE. 

COMPLEMENT takes a four character hex number 
addressed by IR2 and IR3 in RAM, takes the ones complement 
of that number, and stores the result back into the location 
in RAM addressed by IR2 and IR3. 



RAM(IR2 , 



IR3) 



» vyA 



*RAM(IR2, IR3) 



COMPLEMENT occupies 28 bytes of ROM and takes 302.4 micro- 
seconds of computational time. 

COMANGLE takes the angle loaded in IRC thru 
IRF, subtracts it from 90 in hex radians, and stores the 
result back into IRC thru IRF. 



IRC ( C 



F) 




*IR(C -> F) 



COHANGLE occupies 18 bytes of ROM and takes 194.4 micro- 
seconds of computational time, 
b. Shifting 

Multiplication and division of hex-digits by a 
multiple of two was accomplished by shifting the variable 
either left or right the required number of bits. To take 



advantage of this capability, two subroutines were written 
for the navigation program, DIV2IR and DIV2 . 

DIV2IR divides the contents of IRC thru IRP by 
two and stores the result in the RAM addressed by IRO and 
IR1. The division is accomplished by shifting the variable 
right one bit. The original contents of IRC thru IRF are 
unaffected. 




DIV2IR occupies 24 bytes of ROM and takes 259.2 microseconds 
of computational time. 



DIV2 divides the contents of the RAM addressed 
by IR2 and IR3 by two and stores the result in the RAM 
addressed by IR2 and IR3. 



RAM( IR2 , 



IR3) 




IR3) 



DIV2 occupies 33 bytes of ROM and takes 356.4 microseconds 
of computational time. 

c. Data Handling 

The transfer of data within the navigation pro- 
gram was accomplished by the data handling routines. These 
routines saved much memory space in the Executive Routine 
and the Multiply and Cosine Routines by grouping these tasks 
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into separate subroutine calls. There were five data handling 
routines written for the Navigation program. Three routines, 
RAMIRC, RAMIR8, and IRRAMC, were written to transfer data 
between the IRs and RAM. TRANRAM was written to transfer 
data between different locations in RAM. TIME was a special 
routine written to load the proper time interval of one navi- 
gation cycle into the IRs. All five routines were written 
to handle data of four hex-digit size. 

RAMIRC transfers the contents of the RAM addressed 
by IR2 and IR3 into IRC thru IRF. The original contents 
of RAM are unaffected. 



RAM(IR2 , 



IR3) 




*IR(C + F) 



RAMIRC occupies 16 bytes of ROM and takes 172.8 microseconds 
of computational time. 

RAMIR8 transfers the contents of the RAM addressed 
by IRO and IR1 into IR8 thru IRB. The original contents 
of RAM are unaffected. 



RAM (IRO, 



IR1) 




->(IR(8 h. B) 



RAMIR8 occupies 16 bytes of ROM and takes 172.8 microseconds 
of computational time. 
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IRRAMC transfers the contents of the IRC thru 
IRP into the RAM addressed by IR2 and IR3. The original 
contents of IRC thru IRF are unaffected. 



IR(C + F). 




V RAM(IR2 , IR3) 



IRRAMC occupies 16 bytes of ROM and takes 172.8 microseconds 
of computational time. 

TRANRAM transfers the contents of RAM addressed 
by IR2 and IR3 into RAM addressed by IRO and IR1. The 
original contents of RAM addressed by IR2 and IR3 are 
unaffected. 



RAM(IR2 , 



IR3) 




■^RAM ( IRO , IR1 ) 



TRANRAM occupies 8 bytes of ROM and takes 86.4 microseconds 
of computational time. 

TIME is a special subroutine written to load 
the time interval of one navigation cycle into IR8 thru IRB. 
The time interval used in this program was 200 milliseconds. 
To input a different time interval into the navigation pro- 
gram, a new TIME routine is substituted into the program. 

>IR(8 - B) 
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TIME occupies 9 bytes of ROM and takes 97.2 microseconds 
of computational time, 
c. Summary 

Thirteen Common Subroutines were written for 
the Navigation program. Each routine was designed to 
handle the particular needs of the Navigation program. 

All thirteen routines fit in one ROM. A complete listing 
of the Common Subroutines is in Appendix A. 

4 . Summary of Subroutines 

Each operation defined in the operational process 
graph. Figure 6, was successfully programmed within the 
memory and time constraint of the program analysis. Each 
subroutine was written in a modular form allowing for easy 
addition or subtraction of new code. The memory size, 
computational time, and capability of each subroutine is 
listed in Table VII. 

E. EXECUTIVE ROUTINE 

The Executive Routine was written to call up the sub- 
routines in the order described by the operational process 
graph. Figure 6. The Executive Routine established the 
priorities of each function and was designed to make all 
the decisions in the execution of the Navigation program. 
The variables used by the Executive Routine were all stored 
in RAM. The designated location of each variable in RAM 
is shown in Table VIII. 
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Length 

(8-bit Time 



Subroutine 


v/ords) 


(yse 


MULT 


743 


5000 


COS 


255 


5170 


ADDRAM 


12 


289 


ADDRAMIR 


22 


238 


SUBIR 


17 


183 


SUBRAMIR 


32 


3^6 


COMPLEMENT 


28 


302 


COMANGLE 


18 


194 


DIV2IR 


24 


259 


DIV2 


33 


356 


RAMIRC 


16 


173 


RAMI R 8 


16 


173 


IRRAMC 


16 


173 


TRANRAM 


8 


86 


TIME 


9 


97 



Purpose 

[IR(8+B) ] • [IR(C-FF) >IR(C+F) 
COS[IR(C-*F) ] + IR(C+F) 

RAM( 0 , 1 ) + RAM(2 ,3)^RAM(2 ,3) 
RAM( 0,1) + IR(C->-F)-»-IR(C->-F) 
IR(4-7)-IR(C+F)->-IR(C+F) 
RAM(2 > 3)-IR(C-^F)^RAM(2,3) 

~ [ RAM (2 j 3 ) ] +RAM (2,3) 

(90° ) -IR(C-F)->TR(C-tF) 
[IR(C+F) ]/2+RAM(0,l) 

[ RAM (2,3)] /2->-RAM (2,3) 

RAM (2,3) + IR(C-*F) 

RAM (0,1) IR(8-*-B) 

IR(C+F) + RAM (2,3) 

RAM (2,3) + RAM(0 , 1) 

(Nav Cycle Time )->-IR( 8+B) 



Table VII NAVIGATION SUBROUTINES 
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The first operation of the Executive Routine was to 
compute the proper Cosine and Sine of the current true 
heading. Before calling the COSINE subroutine, the Execu- 
tive Routine determined which quadrant the true heading was 
in and set up the input angle in the proper form so that the 
Sine or Cosine could be computed by the COSINE subroutine. 
After calling the COSINE subroutine, the Executive Routine 
stored the output with the proper sign in its location in 
RAM. A graphical display of the decisions made by the Execu- 
tive Routine in computing the Cosine and Sine of the true 
heading is illustrated in Figure 14. 

Next the Executive Routine smoothed the current TAS 
indication with the previous TAS computed. The smoothed 
TAS was then used to compute VAX and VAY as described in 
the Navigation equations and illustrated by the operational 
process graph. 

After computing VAX and VAY, the Executive Routine tested 
to see if the INS was good. If the INS was good, the Execu- 
tive Routine computed the wind acting on the vehicle and 
the distance increments traveled as measured by the INS. 

If the INS was down, the Executive Routine computed the dis- 
tance increments traveled as measured by the AIR-MASS system 
using the last computed wind components. Next the Doppler 
System was tested to see if it was up or down. If the 
Doppler was up, the distance Increments as measured by 
the Doppler were computed. 
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TRUE HEADING PROCESS GRAPH 



The distance Increments outputted from the Navigation 
Program were selected by the Executive Routine in the 
following priority according to the following conditions: 



Priority 



Conditions 



1. 


INS and Doppler 
Smoothed 


INS 


up and Doppler Up 


2. 


INS 


INS 


up and Doppler down 


3. 


Doppler 


INS 


down and Doppler up 


4. 


Air-Mass 


INS 


down and Doppler down 



The Executive Routine was loaded into two ROMs with 
space left for addition of new code. A complete listing 
of the Exeuctive Routine is in Appendix B. 

F. ERROR ANALYSIS 

Microcomputers have a limited arithmetic capability. 

The basic arithmetic operation of the MCS-4 is a four bit 
addition. More complex arithmetic processes had to be 
reduced to this basic operation. It w as therefore very 
important to avoid unnecessary precision throughout the 
calculations. Since the inputs into the Navigation program 
came from instruments whose precision is limited to three 
hexadecimal digits, the choice of four hexadecimal arith- 
metic was considered to be sufficiently accurate. The 
following error bound analysis was performed to show that 
the input errors' dominate the total error. 

The starting point for the error analysis is the opera- 
tional process graph. Figure 6. It was apparent from the 
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operational process graph that the outputs DX and DY are 
symmetric, therefore an analysis of only the computations 
for DX were made. A process graph that involved only the 
operations which have an influence on DX was constructed, 
as shown in Figure 15, from the operational process graph. 
Figure 6. 

•The errors corresponding to each operation in Figure 
were designated e^, egj e^j •••» e g* The errors of 

the inputted data were expressed as e(TH), e(VGXI), e(SDD), 
and e(SHD). Due to the smallness of the errors, the products 
of errors were considered negligible when compared to the 
linear terms. The propagation error was derived as discussed 
in Chapter 2 of Reference 28 in the Bibliography. 

There were two means by which each operation contributed 
to the error propagation. 

1. ) Transmitting the errors which were inputted into 

the operation. 

2. ) Adding an error of its own, which is due to the 

rounding or truncating process which limits the 
number of digits carried to the next operation. 

The transmitted errors were calculated by calculating 
the differentials of the expression. 



d(x-i-y) = dx + dy 



d(xy) = y dx + x dy 




y 



dx 



- x dy 



y 
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Figure 15. PROCESS GRAPH OF DX 



d(sin x) = cos x dx 



d(cos x) = - sin x dx 

The rounding or truncating errors were simply added to 
the transmitted error and thereafter propagated through the 
remainder of the calculation. 

The error after the first operation in Figure 15 is 
given by: 



e (cos(TH) ) = - (sin(TH)e(TH) + e • 

This error is further transmitted by operation three. 



e ( cos (TH ) SDD) = SDD(-sin(TH)e (TH) + e^ + cos(TH)e(SDD) 



Similarly, the error after the second and fourth operation 
is given by: 



e ( sin (TH) SHD) = SHD(cos(TH)e (TH) + e ? ) + sin(TH)e(SHD) 

+ e H 



This gives the total error after operation six as: 
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e(cos(TH)SDD) + e (sin(TH)SHD) + e 6 

In a similar fashion, the error after operation five is: 

e(T.VGXI) = VGXI-e(T) + T-e(VGXI) + e,- 

b 

Combining the above, the total error is expressed as: 

e (DX) = |{e(T-VGXI) + [e ( cos (TH) SDD)+e ( sin (TH) )SHD + e^+e^+eg 

The following is obtained by expanding the expression for 

e ( Dx ) : 

e(Dx) = i{VGXI-e(T)+e(VGXI) •T+e c -+SDD(-sin(TH)e(TH)+e 1 ) 

2 pi 

+ cos (TH ) .e(SDD)+e 3 +SHD(cos(TH)e(TH)+e 2 ) 

+ sin(TK)e(SHD)+e lj +e 6 +e 7 }eg 

The error bound for e(Dx) is computed by using the triangular 
inequalities . 

| e (Dx) | < |{| VGXI| • |e(T) | + | T| • | e(VGXI) | + | SDD | [ | sin (TH )| [e (TH)| 
+ |e 1 |]+ | cos(TlI) | |e(SDD) |+ | SHD | [ | cos (TH ) | | e (TH) | 

+ |e 2 |] + | sin(TH) | | e(SHD) | + | e^ | + | e ^ | + |e,,| 

+ k 6 l + |e ? | } + | eg | 
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The time increment T is controlled by the programmer 
and hence may be chosen to be a power of sixteen, therefore 
the inherent error e(T) = 0. The addition process carried 
out in operations six and seven is carried out exactly in 
the range of values permitted, therefore eg = 0 and e^ = 0. 
This reduces the error bound expression for the total 
propagated error to: 



e(Dx) < |{|T| • |e(VGXI) | + | SDD | [ | sin (TH) | | e (TH) | + lej] 

+ | cos(TH) | | e (SDD) | + | SHD | [ | cos (TH ) | | e (TH) | + |e 2 |] 
+ | sin(TH) | | e (SHD) | + |e 5 | + |e 3 | + |e^|) + |eg| 

The error bound for each operation in Figure 15 due to 
roundoff and truncation was found to be: 



Operation 



Error Bound (Decimal) 



Cosine 

Sine 

Multiplication 

Division 

Subtraction 

Addition 



3 x 10 
3 x 10 
3 x 10 
8 x 10 
8 x 10 
0 



-5 

-5 

-5 

-6 

-6 



Corres. Error 








The error bound for the inputs was obtained from published 
sources. Since an actual maximum error for each system could 
not be found, the (3o) value, 99 - 7 % CEP, was used. The 
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systems onboard the P3C naval aircraft were used as repre- 
senting the current "state-of-the-art" systems In operational 
use today. 



System 


Error (30) 


Corresponding Error 


Inertial Navigation 
System 


5.4 NM/HR 


| e (VGXI) | 


Inertial True Heading 


.5 degrees 


| e (TH) | 


Doppler Along Heading 


1.8 


| e (SHD) 1 


Doppler Across Heading 


3.6 


| e (SDD) | 



The error bound for the Doppler is a function of velocity. 
In order to compute the error bound for the system, an exam- 
ple was used. The error bound analysis was calculated for 
an aircraft with a velocity of >100 knots along the true 
heading and a velocity of 50 knots across the true heading. 

The navigation cycle time used was 150 milliseconds. The 
true heading of the vehicle was 045 degrees. The inputs 
and error bounds inputted into the system were calculated 
as follows: 

TH = 045 degrees = 0. 785*1 radians 

e(TH) = .5 degrees = 0.0087266 

VGXI = 3.8.2 knots 

e(VGXI) = 5.*! knots 

TIME = .15 sec = .0000*116 HR 



SHD = ( MOO knots) (. 0000*116 HR) = .0166*1 nm. 



e(SHD) = (.0166*1 nm) (.018) = .0002995 nm. 

SDD = (50 knots) (. 0000*116 HR) = .00208 nm. 
e(SDD) = (.00208 nm)(.036) = .00007*18 nm. 

Although the values calculated above are dimensionally 
correct, they are not in the proper form for inserting into 
the computer. The MCS-*1 microcomputer used in this system 
works with a l6-bit fixed point variable. In order not to 
lose significant digits in each number, the inputs are all 
normalized. The fixed point number used in the calculation 

is in the form of (. _) , therefore the value of VGXI 

used in the computations is .3182. The time variable is 
normalized to . *1 1 6 0 which results in SHD equalling .166*1. 

Now that the inputs are dimensionally correct and in the 
proper format, they can be substituted into the error bound 
equation as follows : 

|e(Dx) | < i{(.*ll6)(. 005 * 1 ) + (. 0208) [ (SIN( . 785*1 ))(. 0087266) 

+ . 00003 ] + (C0S( . 785 j O )(. 0007 * 18 ) + (.166*1) 

[(C 0 S( . 785 * 1 ) )(. 0087266) + .00003] + (SIM( . 785*1)) 

(.002995) + .00003 + .00003 + .00003) + .000008 
| e ( Dx ) | < .0030798 
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This is the maximum error bound for the total naviga- 
tion system. The maximum error bound equation was broken 
into two parts, one part indicated the maximum propagated 
errors due to the inputs and the other indicated the maximum 
error developed by the navigation program. 

The error bound for the inputs alone is • calculated as 
follows : 

l e(Dx) linput - ! { 1 T I l e(VGXI) I + I SDD 1 [ | SIN (TH) | | e(TH) | ] 

+ |COS (TH) | |e(SDD) | + | SHD | [ | COS (TH) | | e(TH) | ] 
+ j SIN (TH) | | c ( SHD ) | } 

|e(E)x)|input . **l6) ( . 005-O + ( . 208) ( . 70711) ( . OO87266) 

+ 70711) ( • 0007^8) + (. 166*0 (. 70711) (. 0087266) 

+ (. 70711) ( .002996) } 



l G(Dx) knput = ‘ 0030211 



Next the maximum error due to the calculations in the 
navigation program alone were computed. 
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I e (Dx) I < i{(.0208)(. 00003) + (.1664)(. 00003) + .00003 
1 'prog — 2 

+ .00003 + .00003) + .000008 
|e(Dx) l prog £ .0000557 

The results of the error bound analysis indicated that 
the maximum error created in the navigation microcomputers 
computations was only 1.8 per cent of the total maximum 
error. It was concluded from this analysis that the accu- 
racy of the MCS-4 navigation program using a 16-bit fixed 
word data length was well within the limits required for 
the navigation problem. 

The best way to check the results of this error bound 
analysis would be to fly the system in an actual aircraft. 
Since an aircraft was not available for this purpose, a 
detailed FORTRAN simulation program was written to test 
the functions of the navigation program. Section VI dis- 
cusses the FORTRAN simulation program and the results 
obtained from examples tested on this program. 

G. SUMMARY OF PROGRAM 

The navigation program written in this thesis for the 
MCS-Jl microcomputer takes the outputs of the Inertial navi- 
gation system, Doppler navigation system, and Air-Data 



] 27 



system and computes the vehicle’s distance traveled in the 
east/west and north/south directions. The total program 
consists of 1768 instruction words on seven 4001 ROM chips 
broken into an executive routine and fifteen subroutines. 

Two 4002 RAM chips are required to store the data and varia- 
bles used in this program. The total computational time 
required for one navigational cycle is between 36 and 80 
milliseconds depending on the navigational mode used. The 
computational error developed by the navigation program from 
error bound analysis represents only .1 per cent of the 
total error. The program uses a sixteen bit fixed point 
variable which allows it to accept inputs up to 2047-99 
knots w5 th an accuracy of ± .0625 kts. The total cost of 
the one CPU chip, seven ROMs, and two RAMs used by the system 
is $95.00. 

A complete listing of the navigation program written 
in the MCS-4 assembler language is found in Appendices A, 

B, and C. 
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VI . FORTRAN SIMULATION PROGRAM 



The development of the MCS-lJ Microcomputer navigation 
program was a very time-consuming process that involved the 
construction and testing of many subroutines. The testing 
of the total program in its final form was limited by the 
capability of modeling the inputs to the navigation program 
in a machine level language. Hand calculated examples were 
tested and used to debug the completed program, it was de- 
cided to write a detailed bit-by-bit digital simulation of 
the navigation program in FORTRAN on the IBM 360/67 compu- 
ter. The FORTRAN simulation was written so that major 
changes in word length and scaling could be tested without 
reQuiring major rewriting of the microcomputer program. 

Parameter changes were made to the Inertial, Doppler, 
and Air-Mass inputs without incurring the added cost and 
time of developing an actual system. The FORTRAN simulation 
was a comprehensive and flexible means of testing the micro- 
computer program through a series of realistic and unrealis- 
tic tasks. A complete listing of the FORTRAN simulation 
program is found in Appendix D. 

A. DEVELOPMENT 

The FORTRAN simulation program was developed as an exact 
simulation of the navigation program developed for the MCS-'I 
microcomputer. The program was written by utilising the 
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process graphs developed in Section V. The program includes 
a parallel solution of the navigation equations utilizing 
the FORTRAN routines available on the IBM 360/67 computer. 

The results of the solution by the MCS-4 simulation and the 
FORTRAN solution provided a comparison of- the 16-bit fixed 
point solution of the MCS-4 microcomputer program and the 
24-bit floating point routines on the IBM general purpose 
computer. 

The FORTRAN Simulation provided the programmer with the 
ability to input many variables. The outputs of the Iner- 
tial, Doppler, and Air-Mass systems were modeled to include 
system failures and input errors. The outputs of the simu- 
lation indicated how well the programmed microcomputer could 
handle changing variables as compared to a programmed 
general purpose computer. 

The main body of the simulation program was written as 
a subroutine in order that systems programmed or simulated 
in FORTRAN could easily be tested with the navigation routine. 
The main body of the simulation program declares the index 
registers, random access memory locations, and read only 
memory locations as dimensioned arrays. The numbering and 
size of each array corresponds exactly with the addressing 
and size of these areas within the MCS-4. The main body of 
the simulation program prepares the inputs and initial con- 
ditions for the execution of a simulated navigation cycle. 

The hexadecimal Cosine table stored in the MCS-4 ROM is 
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computed and stored in the ROM array by the main body of 
the simulation program. The simulation begins when the 
subroutine NAV is called with the initial inputs. 

Subroutine NAV is a FORTRAN simulation of the MCS-iJ 
navigation executive program. The inputted variables to 
NAV are 32-bit floating point decimal numbers. Each varia- 
ble is converted to a 16 bit fixed point hexadecimal number 
by the subroutine CONVRT, Appendix D. After each input is 
in the proper form, subroutine NAV stores each variable in 
the RAM array at the same location corresponding to the RAM 
addresses in the MCS-iJ, Table VIII. After the inputs have 
been stored in the proper locations and in the proper format, 
subroutine NAV calls a FORTRAN simulation of each subroutine 
called by the MCS-i] navigation executive routine. Each 
subroutine is executed in the same order as in the MCS-^ 
executive program. The 16-bit hexadecimal results are 
printed to provide a check with the results of the actual 
MCS-*J microcomputer navigation program. The NAV subroutine 
lastly calls subroutine CONVD, Appendix D, which converts 
the 16 bit hexadecimal results into floating point decimal 
results which are then compared with the solution of the 
navigation equations by floating point FORTRAN calculations. 

A complete listing of subroutine NAV is found in Appendix D. 

Each common subroutine programmed on the MC3- , I micro- 
computer was duplicated in the FORTRAN simulation. A com- 
plete description of each subroutine and its graphical 
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representation can be found in Section V.D.3* A complete 
listing of the FORTRAN simulation of each of the common 
subroutines is found in Appendix D. 

The multiplication and Cosine subroutines used in the 
FORTRAN simulation were developed from the multiplication 
process graph, Figure 9 , and the Cosine Process graph. 

Figure 13. The Cosine routine uses the Cosine table stored 
by the main body of the simulation program. The multiplica- 
tion table used by the simulated multiply routine was com- 
puted individually for each pair of hexadecimal digits 
called by the multiply routine. 

The last subroutine written for the simulation program, 
subroutine PRINT, was written to incorporate some of the 
features of the MCS-4 Assembler and Interpreter. Subroutine 
PRINT, Appendix D, was developed to allow the programmer to 
dump the contents of the index registers and RAM locations 
at any point during the execution of the simulation. 

A complete listing of the FORTRAN simulation program 
is in Appendix D. 

B. APPLICATION 

The FORTRAN simulation was first used to optimize the 
navigation program. The first navigation program written 
on the MCS-JJ had been designed to handle only positive 
numbers within the true heading limitations of zero to ninety 
degrees. Before expanding the program to handle nl] Quadrant 



and both positive and negative inputs, different ways to 
implement these changes were investigated. The Fortran 
simulation allowed numerous techniques to be attempted with 
the advantages and limitations of each readily apparent 
in the output. The Fortran simulation program was easy to 
change, whereas the MCS-lJ program would have required major 
program changes to implement each technique . After the 
Fortran program had been modified to do the required tasks, 
the MCS-4 navigation program was modified to incorporate 
these required changes. The true heading was used in four 
separate quadrants of ninety degrees each with the Execu- 
tive routine testing the true heading input and computing 
the Cosine and Sine as required. The sign of each variable 
was incorporated into the program by using a ones complement 
scheme with the first bit of the l6-bit word becoming the 
sign bit. The required changes to each subroutine was first 
tested in the simulation program before being programmed 
into the navigation routine. 

The second use of the simulation program was to debug 
the MCS-JJ program in its final form. The simulation program 
was run through a series of examples. The simulated MCS— ^ 
program results were compared v/ith the lortran computed 
results to insure the correctness of the simulation program. 
After the simulation program had been completely debugged, 
the actual MCS-^1 program was run with the same inputs as the 
simulation program. By dumping the contents ol the regi^tei ... 
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and RAM memory at key points in the actual program and 
comparing these values with those computed and printed out 
by the simulation program, the sections of the actual .pro- 
gram in error were easily identified. After changing a 
section in error, the actual program and simulation pro- 
gram were run to see if the correction was correct and also 
see if it affected any of the other sections. This method 
of checking and rechecking proved to be an excellent tool 
for debugging microcomputer programs. The only problem 
encountered in this method was the required time for outputs 
to be printed from the IBM 360/67. 

The third use of the simulation program was to check 
the results of the error bound analysis. The same inputs 
used in the error analysis were inputted into the simulation. 
Simulation runs were made using correct inputs in order to 
compare the results of the MCS- J I calculation with those of 
the Fortran program in order to determine the errors due to 
the navigation program. Simulation runs were made for each 
navigation mode of operation and for each true heading quad- 
rant. Next the same runs were made with the inputs at 
maximum errors. The results of these runs and the comparison 
of the errors developed are summarised in Table IX. 

The results, summarized in Table IX, confirmed the re- 
sults of the error bound analysis. It was noted that the 
greatest computational error occurred when the vehicle tra- 
veled a direct path with constant Inputs. 



This was due to 



Computed X-Loc. Computed Y-Loc. Maximum Error Due to 
Type of Path May. Mode Type Input FORTRAN MCS-4 FORTRAN ICS-4 NAY Program Input 
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TABLE IX. RESULTS OF FORTRAN SIMULATION 



the linear addition of the truncation error when the inputs 
remain constant. It was also noted that the computational 
error was zero when the vehicle returned to the departing 
position indicating that the truncation error cancelled in 
the opposite direction. 
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VII. CONCLUSIONS 



The results of this design study indicated that a micro- 
computer is both fast enough and powerful enough to handle 
the complex task of navigation. The total computational 
time of the navigational cycle, 80 milliseconds, is well 
below the 200 millisecond cycle time used in current naviga- 
tion computers. Table look-up routines for such complex 
tasks as multiply. Cosine, and Sine proved an effective 
means of trading inexpensive memory for computational speed. 
The total hardware cost of the ten LSI cips, excluding a 
board to hold the DIP packages, was $ 95 - 00 . 

The software aids investigated in this study were very 
effective in decreasing programming and debugging time. 

Graph Theory, in the form of process graphs, was an excellent 
means of visualizing the actual flow of the data throughout 
the computations. The initial problem was successively 
broken into smaller and smaller discrete parts until each 
operation could be programmed easily in the assembly language 
of the MCS-4. The operational process graph was then used 
to combine the programmed discrete operations into the re- 
quired navigation computations. The actual writing of the 
program was greatly simplified by this method. 

The error bound analysis of the microprogram was greatly 
simplified by the use of process graphs. The maximum error 
in each of the discrete operations was initially computed. 



137 



The summing of these errors as they propagated throughout 
the program was easily identified with the aid of the opera- 
tional process graph. The testing of actual erros in the 
computation, however, proved very difficult. Sample problems 
were hand calculated with the aid of the process graphs and 
the results were compared with the results computed by the 
actual program. When these results did not agree, both the 
hand calculated procedure and the computer program had to 
be debugged. This procedure was very time consuming with 
only a limited number of tests being made. The FORTRAN 
simulation program was written to aid in testing and debugging 
the navigation program. 

The FORTRAN simulation program proved to be an effective 
aid in testing the navigation program. The set up of the 
simulation program allowed the results of the navigation 
program to be compared at different points in the computa- 
tions with a FORTRAN solution of the same required equations. 
With the aid of the IBM 3 6 0/6 7 computer, many tests were run 
on the FORTRAN simulation program. The results of these 
runs identified those areas in the navigation program that 
needed rewriting. Before expanding the navigation program 
to include negative inputs from different true heading quad- 
rants, the FORTRAN simulation program was changed to incor- 
porate the required changes. After all the changes had been 
tested and optimized on the simulation program, the MCS-^I 
navigation program was rewritten to incorporate those changes. 
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The simulation program proved very effective in optimizing, 
expanding, debugging, and testing the navigation program. 

It is concluded from this design study that many of 
the dedicated computational tasks being done by large 
general purpose computers on board naval ships and aircraft 
can be done by microcomputers. It is recommended that 
research and development of a multi-microcomputer system be 
begun to replace one or more of the costly general purpose 
computer systems currently being used. The development of 
the MCS-808O microcomputer with its increased computational 
speed and power over the MCS-iJ will make this task easier. 

It is also recommended that before programming a micro- 
computer, the desired program should be written in a higher 
level language simulation program before investing in the 
software costs of a microprogram. The development of PL/M, 
derived from PL/I, a high level language for the MCS-8 and 
MCS-8080 greatly increases the ease of programming and 
testing programs for the microcomputer. 

The rapid improvements being made in microcomputers, 
and the current availability of inexpensive microcomputers, 
make it imperative that the Navy begin now in the development 
of systems utilizing this technology. 
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APPENDIX A 

❖ <5 # # ***** £ «««#** £ £ $ * * # ajs # * £ # & * £«**:<: sjs £*£* ^ # $ $$$$$ :* «** # 

THIS IS THE MCS-4 MICROCOMPUTER NAVIGATION PROGRAM 
USEO TO COMPUTE THE PRESENT POSITION OF AN AIRCRAFT 
FROM DISTANCE, VELOCITY, AND AIRSPEED INPUTS FROM 
INERTIAL, DOPPLER, AND/OR AIR MASS SYSTEMS; 

£ £ £ * £ ****** * ****** * * ** * 



***** *********** * * * ** * *** **** * * ****** * £ $ * ** *** ***** *** 
THIS SECTION CONSISTS OF THE COMMON SUBROUTINES 
USED BY THE NAVIGATION PROGRAM AND THE MULTIPLY 
AND COSINE SUBROUTINES; 



# sjc afc 



$$$$$ £ $$ $ : 



~ <v* «v* -r* •v 



THIS SUBROUTINE ADDS THE CONTENTS OF THE RAM MEMORY 
ADDRESSED BY INDEX REGISTERS 0 AND 1 TO THE CONTENTS 
OF INDEX REGISTERS "C" THRU "F" AND STORES THE RESULT 
INTO THE INDEX REGISTERS "C" THRU "F"; 



sjc 



# %V ^ vU «. 
" T V V T ' 



/ %V \U >JU ^ X O# %U ^ X vU X X x X x X X X X %•# 

% ^ ^ «y ^ <Y“ ‘V 'i’* -v 'i' t 5, ^ ^ ^ -p -*p -*p v 



AD DRAM IR: 
CLC; 
SRC PO 
RDM; 
ADD RC 
XCH RC 
INC R1 
SRC RO 
RDM; 
ADD RD 
XCH RD 
INC Ri 
SRC RO 
RDM; 
ADD RE 
XCH RE 
INC Rl 
SRC RO 
RDM; 
ADD RF 
XCH RF 
CLC; 
BBL 0; 



.', -j- -j.. 






£ y,t -Jx y; sj: j(j ?' ;’c ;'; *A «] 

T HI S S UBF OUT I NE TAKES A FOUR CHAR HEX NUMBER ADDRESSED 
BY THE CONTENTS OF INDEX REGISTERS 0 AND 1, AND ADDS 
IT TO A FOIJP CHAR HEX NUMBER ADDRESSED BY INDEX RFGS 
2 AND ?. , AND LOADS THF RESULT INTO RAM MEMORY ADDRESS 
BY INDEX REGISTERS 2 AND 3; 

* i £***** * * * * * * ** * * * * if. * * if * * * * * ** *** ****** * X* * ****** if if 



AD DP AM : 

CLC; 

LDM 12; 

XCH R5; 

AD: 

SRC RO; 
RDM; 

SRC R2 ; 
ADM; 

WRM: 

INC F’.l ; 

INC R3 ; 

ISZ R5 AD; 
BBL O; 



l'lO 



THIS SUBROUTINE TAKES A FOUR CHARACTER HEX NUMBER 
LOCATED IN INDEX REGISTERS "C" THRU "F", DIVIDES 
THIS NUMBER BY TWO AND STORES THE RESULT IN RAM ; 

*##$#****£** a***# *****#*£#***£#*##*****£#*****££***«£#* 



DI V2IR: 

CLC; 

SRC RO; 
LD RF; 
PAR; 
WR.M; 

LDM 14; 
XCH Ri; 
SRC RO; 
LD RE; 
RAR r 
WRM; 

LDM- 13; 
XCH Ri; 
SRC RO; 
LD rd; 
RAR; 
WRM; 

LDM 12; 
XCH p 1 ; 
SRC RO; 
LD RC; 

par; 

WRM; 

BBL o; 






J/ 4- J-« 4# ^ ^ ^ -J 

>|< #1* ^ *]*> r|i* >p #p rf. -I- 

THIS SUBROUTINE DIVIDES THE CONTENTS OF RAM MEMORY 
ADDRESSED BY INDEX REGISTERS 2 AND 3 BY TWO: 



J» V* ■>*» »•* j/ <ju «JU **' »u V 1 .* j< »i» «>/ «v ' 

*0T •T- *t* t* *r* “V* <v •v* *v 4' “v -V -i' Or' 'r» *v ** 



y- %«* J- »<• -J* 

*>' 4* »i* *r '<* t ¥ * 1 * *v 



* sj: # % s}c # ijc # # £ sjc 5jc £ # # # ## # # 



D 1 V2 : 

CLC; 

SRC R2 ; 

RDM; 

RAL; 

JCN 2CY PDI V ; 
CLC; 

RAR; 

stc; 

jun'ndiv; 

PDIV: 

RAR; 

NO IV: 

RAP; 

WRM; 

LDM 14; 

XCH R3; 

SRC P.2; 

RDM; 

RAR; 

WRM; 

LDM 13; 

XCH R3; 

SRC R2 ; 

RDM; 

FAR; 

WF:M; 

LDM 12 ; 

XCH P3; 

SRC R2 ; 

RDM; 

RAR; 

WRM; 

BBL o; 



l'l] 



# #*%%%*&**#**%******##* *£&£ £#£#:><* 



THIS SUBROUTINE LOADS THE DESIRED TIME CYCLE FOR THE 
NAVIGATION ROUTINE INTO INDEX REGISTERS "8" THRU "B"; 



TIME: 



LDM 


0; 


XCH 


R8 ; 


LDM 


0; 


XCH 


R9 ; 


LDM 


0; 


XCH 


P.A ; 


LDM 


4; 


XCH 


RB ; 


BBL 


0; 



^-. sjc^: :J:s(: sjcsjt # * * jf: # sfrj^sf; * $ s}e s^sjcsjc :>; $ # a!:# 3j: sje 



THIS SUBROUTINE LOADS A FOUR CHAR HEX NUMBER FROM 
INDEX REGISTERS "C" THRU "F" INTO RAM MEMORY ADDRESS 
BY INDEX REGISTERS 2 AMD 3; 



■ V'c^'r tv^v vr'i'r*'T r ' i'v'pt'T Tv'vrvv ttt'i'Tt Vr 'i ,, f'i'T r v «v- «*r“ +r ^ *r -t* *v r v 'i' <T“r 



IRRAMC : 

SRC R2; 
LD RC; 
WRM; 

INC R3 ; 
SRC R2 ; 
LD RD; 

wrm; 

INC R3 ; 

r r» r r* ^) . 

Cr\U r\ t. $ 

LD RE; 
WRM; 

INC R3 ; 
SRC R2 ; 
LO RF; 
WRM; 

BBL 0; 



ij; •%. £ £ # # * 5^ # ^ ^ ^ ^ s* # ?Jt $ * # :!s ;Jc * si sjc # sjt # s'c * 3j<* ^ s,V sfc sisjc sf: # 4c :\c ^ # :£ # sjr 



THIS SUBROUTINES LOADS A 
INDEX REGISTERS "C" THRU 



FOUR CHAR HEX NUMBER INTO 
"F» FROM RAM MEMORY LOCATION 



ADDRESSED BY THE CONTENTS OF INDEX REGISTERS 2 AND 3: 

3* £ j}: * * # <•. £ £ * # * * £ ** *4: & * >}e * $ # « £ * $ # « £ *# * *# ^ 



RAM IRC : 

SRC R2 ; 
RDM; 

XCH RC ; 
INC R3 ; 
SRC R2 ; 
RDM; 

XCH RD; 
INC R3 ; 
SRC R 2 ; 
RDM; 

XCH RE ; 
INC R3 ; 
SRC R2 ; 
F DM ; 

XCH RF; 
BBL 0; 



Hi 2 



# * £ ** * * * * * * * *# * £ # * **$ * * * * ** ** * * **** 

THIS SUBROUTINE LOADS A FOUR CHAR HEX NUMBER INTO 
INDEX REGISTERS 8 THRU "B" FROM A RAM MEMORY LOCATION 
ADDRESSED BY INDEX REGISTERS 0 AND 1; 



RAM I R8 : 

SRC RO; 
RDM; 

XCH R8 ; 
INC Rl; 
SRC RO; 
RDM; 

XCH R9 ; 
INC Rl ; 
SRC RO; 
RDM; 

XCH RA; 
INC Rl ; 
SRC RO; 
RDM; 

XCH R6; 
BBL o; 



sC %V V/ \!r ^ «X> > V %U V/ %*L* -s* 1 -* X ■**/■ %** V # »** **- «**> 'Jf J# *.C v*# ^ %U> ^ *U J- ^ vV 'V s*<* ^ vU %*» V- \V 

^ ^ v t* ^ ^ / w r V 'r t *** # r ^r »r v 'i* r # i* v *v '»* r i' n' *i % v "r V ^ ^ V -r 



THIS SUBROUTINE SUBTRACTS THE CONTENTS OF INDEX 
REGISTERS "C" THRU "F" FROM THE CONTENTS OF INDEX 
REGISTERS 4 THRU 7 AND STORES THE RESULTS IN INDEX 
REGISTERS "C" TH»U "F"; 



«<* J> %V <X> %V *?* »Jl» »W- *'< j« <JUr *U J/ «JU »V 'V ^ *•« »*» %•* »'< ^ *** j» »V »U »(» «Ar *V **» >•*— Jf «<•■' >4 »•» »*• 4 ^V *V «A* »V J. ^ vV *'< *Xr <JLr mJLr 

v *r -t' v «v* v *v *,- v 'A *>' 'i* 'r • ' *•* ^ ' *%' -A v •,' 'A *i» -r- y r v •>» -i' v v 'c- v w r -r -V v n* v v*r v v v n- v •'i' *»- -i- -v* 



SUB1R: 

LD R4; 
SUB RC ; 
XCH RC ; 
CMC; 

LD R5; 
SUB RD; 
XCH RD; 
CMC; 

LD R6 ; 
SUB RE; 
XCH RE ; 
CMC; 

LD R7; 
SUB RF ; 
XCH RF ; 
CMC; 

BBL 0; 



1'I3 



# if # £ £ # £ ^ * # # * £ A # A # & * £ £ £ # A -r £ * ;*c :* **£ $ * * * * * s$s # * £ * # £ * ^c * £ * £ 

THIS SUBROUTINE SUBTRACTS THE CONTENTS OF INDEX REGS 
"C" THRU "F" FROM THE CONTENTS OF RAM MEMORY ADDRESSED 
BY INDEX REGS 2 AND 3 WITH THE RESULT LOADED INTO RAM 
MEMORY ADDRESSED BY INDEX REGISTERS 0 AND 1; 

% * # =£ * =? £ =S * =£ # # $ # * £ £ ❖ * $ # t- # *f # $ $ # £ « « =£ *■ « * A £ £ £ £ * # * $ * * * # Jfc # $ * * £ 

SUBRAM IR: 

CLC; 

SRC R2 ; 

RDM; 

SUB P.C ; 

SRC RO; 

WRM; 

CMC; 

INC R3 ; 

INC Rl; 

SRC R2 ; 

RDM; 

sub p.d; 

SRC ro; 

WRM; 

CMC; 

INC P3; 

INC Rl; 

SRC R2 ; 

RDM; 

SUB RE ; 

SRC ro; 

WRM; 

CMC; 

INC R3 ; 

INC Rl ; 

SRC R2; 

RDM; 

SUB RF ; 

SRC RO; 

WRM; 

CMC; 

BBL o; 



£ * # # i'. # # jjc A if * ifi, c # # # £ # £ #:)5 A if 3}c if if ^ if ;Jc # sic ;Js # * £ if if if :£ sj< sjc # # if # # 

THIS SUBROUTINE TRANSFERS THE CONTENTS OF THE RAM 
MEMORY ADDRESSED BY INDEX REGISTERS 2 AND 3 INTO 
RAM MEMORY ADDRESSED BY INDEX REGISTERS 0 AND 1; 

if if it A sjt if if itif if s*c if * # £ & if # ^ 5}: :Jc i< if # y.t s{s # j{: ^ if sj: # A ii # ^ £ j|c if jjc # # $ if # £ £ 



TR ANRAI-1 : 

SRC R2 ; 

RDM; 

SRC RO ; 

WRM; 

INC Pa ; 

ISZ R 3 TRANRAM; 

seL o; 



THIS SUBROUTINE TAKES A FOUR CHARACTER HEX NUMBER THAT 
IS LOADED IN RAM MEMORY, TAKES THE COMPLEMENT CF THAT 
NUMBER AND RETURNS IT BACK TO RAM MEMORY; 

£ 3(r 5 }: !*c 4! $ * * * # £ sj« ■if. 5* # -if. :•< * # £ # # # ^ 5* * * * 5* :,'« * * # $ £ £ # * # # =£ ❖ ❖ # * # £ £ * ❖ 

COMPLEMENT : 

CLB; 

SRC R2; 

RDM; 

C M A ; 

I AC; 

WRM; 

INC R3 ; 

SRC R2 ; 

LDM 0; 

XCH ro; 

RDM; 

CM A; 

ADD RO; 

WRM; 

INC R3 ; 

SRC R2 ; 

RDM; 

CMA; 

ADD RO; 

WRM; 

INC R3 ; 

SRC R2 ; 

RDM; 

CMA; 

ADD RO; 

WRM; 

CLB; 

B8L o; 



vC «// *x# %v ^ -X- v'/ ^ vC* %*/ >V x# %V \C %V ^ 

•t' ^ t v n** *r* 'S *<+ ^ 'r 'i ' t *r V # i* t t* ^r* ^ n 



> \U * f / «v 



v<v %V ^ <V <JU ^ %C- X X J f 

#,» r|% ^ ^j •> <',*> 



THIS SUBROUTINE 
LOADED IN INDEX 

%*# %U JU %U %♦/ sU «JU %V X/ 4/ 4/ 

T 3 ? ^'T V ^ -IT 'r T V 'r 'i v T V 



FINDS THE COMPLEMENT OF AN ANGLE 
REGISTERS "C" THRU "F"; 

JU 4/ 4/ 4/ *U 4/ 4« 4/ 4^ 4/ 4# ^ 4/ 4/ 4^ 4 f 4^ 4^ 4^ 4/ 4# %V 4/ 4^ 4# %V vV 4# 4/ 4* 4# JU %C 

v *i' ^ t v 'i* v y v ^ -i'* v v v *r- ^ ^ v -v* -r^ V 



COMANGLE: 

CLC; 

LDM 2; 

SUB RC; 

XCH RC ; 
CMC; 

LDM 2; 

sub rd; 

XCH RD; 
CMC; 

LDM 9; 

SUB RE; 

XCH RE; 
CMC; 

LDM I ; 

SUB RF ; 

XCH RF ; 
CLB; 

BBL 0; 

ENUSU3: NOP; 



APPENDIX B 

*##* 5 !: ££##*;( 

THIS SECTIGN IS THE EXECUTIVE PART OF THE 
NAVIGATION PROGRAM; 

# # £ if =4= * # * <S $ # * * * * * # # 5^ £ $ * # =* * £ * * sfc $ * « * # # $ aj: * £ * xjc >jc sjc * # $ * * # # * # 



PROG: REM 

CLC; REM 



LDM 


2; 




REM 


XCH 


RC 






LOM 


2; 






XCH 


RD 






LDM 


9; 






XCH 


RE 


5 




LDM 


I; 






XCH 


RF 


• 

t 




FIM 


R2 


92; 


R E M 


FIM 


RO 


72; 


REM 


JMS 


TRANRAM; 


REM 


FIM 


RO 


56; 


• REM 


FIM 


R2 


72; 


REM 


JMS 


SU 


BRAMIR; 


RFM 


JCI1 


NZCY QO; 


REM 


FIM 


RO 


72; 


REM 


FIM 


R2 


56; 




JMS 


SU 


BRAMIR; 


REM 


JCN 


NZCY Ql; 


REM 


FIM 


RO 


56; 


REM 


FIM 


R2 


72; 




JMS 


SUBRAMIR; 


REM 


JCN 


NZCY Q2 ; 


REM 








REM 


F 1 1-1 


R2 


56 ; 


REM 


JMS 


RA 


MI RC: 


REM 


JMS 


CO 


SINE; 


REM 


FIM 


RO 


79 : 


REM 


JMS 


D I 


V2 I R ; 


REM 


FIM 


R2 


56 ; 




JMS 


RA 


MIRC; 


REM 


JMS 


CO MANGLE; 


REM 


JMS 


CO 


sine; 


REM 


FIM 


RO 


31 ; 


REM 


JMS 


D I 


V2 1 p ; 


REM 


JUN 


THSET ; 


REM 


: 






REM 


CLC; 


1 






FIM 


R2 


56; 




JMS 


CO 


MPLEMFNT ; 


REM 


FIM 


R2 


56 ; 




JMS 


RAMI RC; 




JMS 


CO 


SINE; 


REM 


FIM 


RO 


79 ; 




JMS 


Di 


V2I f; : 


. REM 


F I M 


P.2 


72 ; 




JMS 


RA 


MI RC; 




JMS 


CC 


S I N E ; 


REM 


F I M 


RO 


3 1 ; 




JMS 


01 


V2IR; 


P E M 


FIM 


R. 2 


2 8 ; 




JMS 


CO 


MPLEMFNT ; 


REM 


JUN 


TH 


SET ; 


REM 


l 






REM 


CLC; 


r 






FIM 


R2 


72 ; 


P EM 


JMS 


CO 


y.PL PMr.NT ; 


F I M 


R 2 


72; 




JMS 


RA 


MI RC; 





BEGIN NAV CYCLE: 
CLEAR CARRY; 

LOAD INDEX REGISTERS 
90 DEG AND TEST THE 
TH QUADRANT; 



LOCATION OF TRUE HEAD; 
LOCATION OF TH WORK SPC; 
TRANSFER TH TO WORK SPC; 
LOCATION OF SECOND 
WORK SPACE; 

TH INPUT; 

TH-90DEG LOADED IN RAM; 
JU^P IF TH IN QUAD 0; 
TEST IF TH IN QUAD 1; 

TH-90DEG ; 

JUMP IF QUAD U 
NOW TEST FOR QUAD 2; 

SUB ANOTHER 90DEG; 

JUMP IF QUAD 2; 

TH IN QUAD 3; 

tu-ovahcq • 

LOAD' I R t S’ WITH TH-270; 
COS (TH-270) ; 

MEMORY SPACE FOR SINE; 
SHI FT FOR SIGN BIT ; 

LOAD IR WITH TH-270; 
COMPLEMENT TH; 

COS OF TH IN QUAD 3; 
MEMORY LOCATION FOR COS; 
SHI FT FOR SIGN SIT; 

QUAD 3 TH SET 
JUMP TO MAIN PROGRAM; 
FIND COS AMD SINE FOR TH 
IN QUAD 2; 



SHIFT TH INTO 0-90DEG ; 

FIND S I N ( T H ) ; 

SHIFT FOR SIGN BIT; 



COS( TH) ; 

SHIFT FOR SIGN BIT; 

COS NEGATIVE IN QUAD 2; 
JUMP 10 MAIN PROGRAM; 
COS AND SIN IN QUAD I; 



] 1|6 



SHIFT TH INTO 0-90DEG ; 



JMS COSINE; 

FIM RO 31 ; 

JMS DIV2IR; 

FIM R2 28; 

JMS COMPLEMENT; 
FIM R2 56; 

JMS RAMI RC ; 

JMS COSINE; 

FIM RO 79; 

JMS DI V2 I R ; 

FIM R2 76; 

JMS COMPLEMENT; 
JUN THSET ; 

QO : 

CLC; 

FIM R2 56; 

JMS COMPLEMENT; 
FIM R2 56; 

JMS RAMI RC ; 

JMS COSINE; 

FIM RO 79; 

JMS DIV2IR; 

FIM R2 76; 

JMS COMPLEMENT; 
FIM R2 72; 

JMS RAMI RC ; 

JMS COSINE; 

FIM RO 31; 

JMS OIV2IR; 



REM FIND COS(TH) ; 

REM SHIFT FOR SIGN BIT; 

REM COS (TH) NEG IN QUAD 1 ; 

REM FIND SIN(TH) ; 

REM SHIFT FOR SIGN BIT; 

REM SIN(TH) NEG IN QUAD 1; 
REM JUMP TO MAIN PROGRAM ; 
REM COS AND SIN IN QUAD 0; 

REM SHIFT TH INTO 0-90DEG ; 

REM SIN(TH); 

REM SHIFT FOR SIGN; 

REM SIN(TH) NEG IN QUAD 0; 

REM COS(TH) ; 



THSET: 

CLB; 



FIM 


RO 32; 


FIM 


R2 12; 


JMS 


ADDRAM; 


FIM 


R2 15; 


JMS 


DI V2 ; 


FIM 


RO 12; 


JMS 


RAMI RC; 


FIM 


R2 28; 


JMS 


RAMI RC; 


JMS 


MULT; 


DY I NS : 


FIM 


R2 o; 


SRC 


R2 ; 


RDM; 


JCN 


ZAC DNIY; 


IN SUP: 


FIM 


R2 60; 


JMS 


I KRAMC; 


FIM 


P 2 60; 


JMS 


COMPLEMENT ; 


F i M 


RO 63; 


FIM 


R2 60; 


JMS 


ADDRAM; 


FIM 


R2 60; 


F I M 


RO 52; 


JMS 


ADDRAM; 


FIM 


R2 oO; 


f I M 


PO 52; 


JMS 


ADDRAM; 


FIM 


R 2 60; 


F IM 


P.0 52 ; 


JMS 


ADDRAM; 


r IM 


R2 63 ; 


JM S 


Div? ; 


F I M 


P 2 63? 


JMS 


D I V2 ; 


F ) M 


PO 52; 


FIM 


R2 60; 


JMS 


TRANRAM; 



REM START MAIN PROGRAM; 

REM CLEAR AC AND CY ; 

REM MEM LOC OF INPUT TAS; 
P.EM MEM LOC OF SMOOTH TAS; 
REM INPUT TAS ADDED TO 

ci'onTu T A C • 

REM HfVlOC'of’ SUM; 

REM DIV SUM BY 2; 

REM MEM LOC NEW SMOOTH TAS; 
REM LOAD TAS INTO IR8-IRB; 
REM MEM LOC OF COS(TH); 

REM LOAD C-OS(TH) ; 

REM MULT TAS*C0S (TH)=VAY; 
REM TEST IF INS GOOD; 

REM MEM LOC INS FLAG; 

REM LOAD FLAG INTO AC; 

REM JUMP IF INS DOWN; 

REM INS UP COMPUTE WIND; 

REM MEM LOC VAY; 

REM LOAD VAY INTO RAM MEM; 
REM MEM LOC CF VAY; 

REM NFG CF VAY; 

REM MEM LOC OF VGYI; 

REM SIT WIND SMOOTHING ROUT 
REM (VGYI-VAY) =VWYR; 

REM MEM LCC VWYR; 

REM MEM LOC VWY ; 

REM Vi-.YR+VWY; 

REM MEM LOC (VWYR+VWY); 

REM MEM LOC VWY; 

REM ( VWYR+VWY) + VWY ; 

REM MEM ( VWYF +VWY+VWY ) ; 

REM ME M VWY; 

REM ( V ri YR + VWY + VWY ) + VWY : 

RLM MEM (VWYR+3VWY) 

REM (VWYR +3 VWY J / 2 ; 

REM (VWYF 3 VWY ) / 4 ; 

REM MEN LCC VWY; 

REM MEM LOO ( V.% Y R+ 3 VWY ) / A ; 
REM VWY= (VwYR + .iVWY ) M ; 



DY I : 




REM 


FIN 


R2 68; 


REM 


JMS 


RAM IRC; 


REM 


JMS 


time; 


REM 


JMS 


MULT; 


REM 


FIM 


R2 24; 


REM 


JMS 


I rpamc; 


REM 


JUN 


up iy; 


REM 


ONI Y: 




REM 


FIM 


RO 52; 


REM 


JMS 


addramir; 


REM 


JMS 


time; 


REM 


JMS 


MULT ; 


REM 


FIM 


R2 24; 


REM 


JMS 


1 RRANC; 


REM 


UP I Y: 




REM 


FIM 


R2 7 6 ; 


REM 


JMS 


RAM IRC; 


REM 


FIM 


p.0 12; 


REM 


JMS 


RAMI R8 ; 


REM 


JMS 


MUL T ; 


REM 


DX INS: 




REM 


F I M 


R2 0; 


REM 


SRC. 


R2; 




RDM; 




REM 


JCN 


ZAC DM IX; 


REM 


FIM 


ft 2 44; 


REM 


JMS 


I RPAMC; 


REM 


FIM 


R 2 44; 


REM 


JMS 


COMPLEMENT; 


REM 


FIM 


R2 44; 


REM 


FIM 


RO 64; 


REM 


JMS 


A DDR AM; 


REM 


FIM 


P.2 44; 


REM 


FIM 


<n / r \ » 

r u *to j 


P. E M 


JMS 


ADDPAM; 


REM 


FIM 


R 2 44; 


REM 


FIM 


RO 46; 


REM 


JMS 


AD DR AM; 


REM 


FIM 


R2 44; 


REM 


FIM 


RO 48; 


KEM 


JMS 


ADDPAM; 


REM 


F I M 


R 2 47; 


REM 


JMS 


D I V2 ; 


REM 


FIM 


P 2 47; 


REM 


JMS 


PI V2; 


REM 


F I M 


RO 48; 


REM 


FIM 


R2 44; 


REM 


JMS 


1 RANRAM; 


REM 


DXI : 




REM 


FIM 


R2 64; 


REM 


JMS 


RAMI PC; 


REM 


JMS 


TIME; 


REM 


JMS 


MULT ; 


REM 


F I M 


R2 6; 


REM 


JMS 


IRRAMC; 


REM 


JUN 


UP IX; 


PEM 


ON IX: 




F EM 


FIM 


RO 48; 


REM 


JMS 


A D D° A M I P. ; 


REM 


JMS 


TIME; 


REM 


JMS 


MULT ; 


REM 


FIM 


F 2 6; 


REM 


JMS 


IRRAMC ; 


REM 


UP IX: 




REM 


FIM 


R2 I ; 


REM 


SRC 


R2 ; 




POM; 




REM 


JCN 


/.AC POWND; 


RDM, 


FIM 


R2 76; 


REM 


JMS 


RAMI RC; 


REM 



COMPUTE OYI ; 

MEM LOC VGY I ; 

VGY 1 LOADED IPC-IRF; 
TIME INC IR8-IRB; 

(VGYI )*(TIME)=DYI ; 

MEM LCC DY; 

DY=DYI ; 

JUMP TO COMPUTE 
VWX AND DXI ; 

COMPUTE DYW; 

MEM LOC VWY; 

( VAY+VWY)=VGYW; 

LOAD TIME INCR; 
(VGYW)*(TIME)=DYW; 

MEM LOC DY; 

DY = DYW ; 

COMPUTE VWX AND DXI; 

MEM LOC SIN(TH) 5 
LOAD S I N ( TH ) ; 

MEM LOC TAS; 

LOAD TAS; 

(TAS)*(SIN(TH) )=VAX; 
TEST IF INERTIAL GOOD; 
INERTIAL FLAG; 

READ FLAG INTO AC; 

JUMP IF INERTIAL DOWN; 
MEM LOC VAX; 

STORE VAX; 

MEM LOC VAX; 

NEGATE VAX; 

(-VAX ) ; 

( VGXI ) ; 

( VGXI )+(-VAX ) =VWXR; 

(vwxi ; 

/ W 1„; V \ • 

( VWXP l-t-(vwx) ; 

( v wxr-4-vwx I ; 

(VWX) ; 

(vwxR+vuxi+tvwxi; 

( VWXR+2VWX) ; 

(vwx); 

( VWXF+2VWX ) + ( VWX) ; 

( VWXR+3VWX) ; 

(VWXR+3VWX) /2; 

( ( VWXR+3VWX I / 2 ) ; 

( ( VWXR+3VWX ) / 2 ) / 2 ; 

(vwx ) ; 

( ( VWXR+3VWX ) /4) ; 

V W X = ( (VWXR + 3VWX )/4) ; 
COMPUTE DXI; 

(VGXI); 

LOAD VGXI; 

LOAD TIME INCR; 

(VGXI l*(TIME)=DXI; 

DX; 

D X = D X I ' 

COMPUTE DYD AND DXD; 
COMPJTE DXW IF INS DOWN 
(VWX) ; 

(VWX ) + (VAX) ^VGXW; 

LOAD TIME INCR; 

( VGXW)*(TIMF:)=DXW; 

(DX) ; 

dx=pxw; 

COMPUTE DYD AND DXD; 
(DOPPLER FLAG) ; 

LOAD DGPP FLAG INTO AC; 
JUMP IF OOP P DOWN; 

( SIN(TU) ) ; 

LOAD SIN( TH) ; 
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FIM RO 84; 

JMS RAMI Re; 

JMS MULT; 

FIM R2 96; 

JMS IRRAMC; 

FIM R2 96; 

JMS COMPLEMENT; 
FIM P 2 20; 

JMS RAMIRC; 

FIM RO 80; 

JMS RAMI R 8 ; 

JMS MULT; 

FIM RO 96; 

JMS ADDRAMIR; 

FIM R2 140; 

JMS IRRAMC; 

FIM R2 28; 

JMS RAMIRC; 

FIM RO 84; 

JMS RAMIR8; 

JMS MULT; 

FIM R2 100; 

JMS IRRAMC; 

FIM R2 76; 

JMS RAMIRC; 

FIM RO 80; 

JMS RAMIR8; 

JMS MULT; 

FIM RO 100: 

JMS ADDRAMIR; 

FIM R 2 0; 

SRC R2 ; 

RDM; 

JCN ZAC DOWN I ; 
FIM RO 8; 

JMS ADDRAMIR; 
FIM RO 127; 

JMS D I V2 I F ; 

FIM R2 124; 

FIM RO 8; 

jms tranram; 

FIM RO 24: 

FIM R2 140; 

JMS ADDRAM; 

JMS DIV2J 
FIM RO 24; 

JMS TRANRA”; 

jun do wnd ; 

DOWN 1 : 

FIM R2 8; 

JMS IRRAMC; 

FIM R2 140; 

FIM P.0 24 ; 

JMS TRAN PAM; 

DOWND: 

JUN PRCG; 



REM 

REM 

REM 

REM 

REM 

PEM 

REM 

REM 

REM 

RBM 

REM 

REM 

PEM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

PEM 

REM 

REM 

REM 

REM 

R E M 

REM 
REM 
REM 
REM 
REM 
REM 
P.EM 
REM 
PEM 
REM 
REM 
REM 
REM 
REM 
PEM 
R E M 
REM 

PEM 

REM 

PEM 

REM 

REM 

REM 

REM 



(SDD) ; 

LOAD SDD; 

(SOD)*(SIN(TH) ) ,* 

( SDD*SIN(TH ) ) ; 

STORE SDD* SIN(TH); 

( SOD*SIN(th ) ) ; 

NEG SDD*S I N ( TH ) ; 
(COS(TH) ) ; 

LOAD COS< TH) ; 

(SHD) ; 

LOAD SHD; 

( S HD ) * (COS ( TH ) ) ; 

-( SDD*SIN(TH) ) ; 
SHD*COS(TH)-SDD*SIN(TH) 



= DYD ; 

(DYD) ; 

STORE DYD; 

(COS (TH) ) ; 

LOAD COS(TH) ; 

(SDD) ; 

LOAD SDD; 

( SDD ) *CQS ( TH ) ; 

WORK AREA; 

STORE SDD*COS(TH); 
(SIN(TH) ) ; 

LOAD SIN(TH) ; 

(SHD) ; 

LOAD SHD! 

(SHD)*(SIN(TH) ) ; 
(SDD*COS(TH) ) ; 
SHD*SIN(TH)+SDD*COS(TH) 

= dxd; 

TEST INERTIAL FLAG; 

LOAD INS FLAG INTO AC; 
JUMP IE INS DOWN; 

(DX) ; 

( DX ) + ( DXD ) ; 

(DX WORK AREA); 

LOAD (DX+DXD)/2; 

( D X + DXD ) /2 J 

(DX) ; 

DX= ( DX+DXD ) / 2 ; 

(DY) ; 

(DYD) ; 

( DY ) + ( DYD ) ; 

( DY + DYD ) /2 ; 

(DY) ; 

DY = ( DY+DYD ) / 2 ; 

JUMP TO END OF CYCLE; 

IF INS DOWN USE DXD 
AND DYD; 

(DX) ; 

ox=dxd; 

(DYD) ; 

(DY) ; 

DY = DYD ; 

END OF NAY CYCLE; 

RETURN AND START 
NEXT NAV CYCLE; 
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APPENDIX C 






THIS SECTION CONTAINS THE NAVIGATION SUBROUTINES 

* £ * * $ £ a*** 5* Sis * * * * *$$* * ***** $ ***** * * # * * * £*** * #« $ * 



af::^ ajcfcajc :}c 

THIS IS THE TABLE OF VALUES USED BY THE COSINE ROUTINE 
THE FIRST COLUMN CONSISTS OF THE VALUE OF THE COSINE 
FROM 0 TO 90 DEGREES IN INCREMENTS OF .08 HEX RADIANS; 
THE SECOND COLUMN CONSISTS OF VALUES USED TO 
INTERPOLATE BETWEEN VALUES IN THE FIRST COLUMN; 



CON "FFFF" 
CON "OEFF" 
CON "08 FF" 
CON "OEEF" 
CON "IOEF" 
CON "2 ECF" 
CON "38BF" 
CON "6E9F" 
CGN "B 0 8F 11 
CON "1F5F" 
CON "A93F" 
CON "60 IF" 
CON "63 EE" 
CON "A2 BE" 
CON "3E7E" 
CON "364F" 
CON "9A0E" 
CON "7BCD" 
CON "F88D" 
CON "F24D" 
CON "B9 FC " 
CON "3 DAC " 
CON "9D3C" 
CON "CAOC" 
CON "05BB" 
CON "5C5B" 
CON "COOB" 
CON "82AA" 
CON "814A" 
CON "0ED9" 
CON "1879" 
CON "BF 09" 
CON "1 5 A8" 
CON "5838" 
CON "89C7" 
CGN "B 857" 
CON "26 E6" 
CON "Cl 76" 
CGN "DBFS" 
CON "6485" 
COM "9605" 
CON "8194" 
CON "4614" 
CUN "0A93" 
CGN " CC 13" 
CCN "OF 92" 
CON "6022" 
CCN "4 I, A 1 " 
CCN "C 1 2 1 " 
CCN "F 1 AO" 



CON "0040" 
CON "FF BO" 
CON "BF31" 
CON "2F31" 
CON "1E32" 
CON "8CB2" 
CON "4A33 " 
CON "37B3" 
CCN "3334" 
CON "3EA4" 
CON " F 72 5 " 
CON "70A5" 
CON "9716" 
CON "2D 8 6" 
CON "1107" 
CON "4377" 
r n v mo3E7” 

CON "FI 58" 
CGN "3EB6" 
CON "4829" 
CON "1099" 
CON "85F9" 
CON " 68 5 A " 
CON "C8B A" 
CON "661 B" 
CON "417B" 
CON "49CB" 
CON "5E1C" 
CON "607C" 
CON "5FEC" 
CON "1B00" 
CON "835D" 
CON "B69D" 
CON "7ADD" 
COM "CG1E" 
COM " 83 5 E " 
CON "BA 8 F." 
COM "AEBE" 
CON "2EEE" 
CON "4A1F" 
CON "924 F" 
COM "276F" 
CON "D78F" 
CON "A4AF" 
CON "90 BF" 
CON "82 OF" 
COM "83 E F " 
CCN "90FF" 
CON "A9FF" 
CCN "BEFF" 
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***#*£$££**££*$£*££££#*❖$#$£:£££#****£❖**$#£ $❖ £*$$#*#** 
THIS SUBROUTINE COMPUTES THE COSINE OR SINE OF A VALUE 
INPUTTED IN HEX RADIANS USING A SECOND ORDER NEWTON 
DIVIDED-DIFFERENCE INTERPOLATING ROUTINE WITH TABLE 
LOOKUP; 

S)C « 3(t 3? * 5^ * ❖ ❖ ^ ^ * 5^ * 4= * * =|: ^ ^ ^ * 5.= ❖ « ^ ^ ^ ^ * ^ * * ❖ ❖ ❖ =Jt ❖ ❖ ❖ * 



COS INE :CLB; 

FIN RO 
SPC 
LDM 

wrm; 

LDM 
XCH 

rar; 

XCH 
RAR; 
XCH 
JCN 
LDM 
XCH 
NCR: 

LD RD; 

ral; 

JCN 
JUN 
CR : 

LDM 
XCH 
JUN 
ACR: 

CLC; 
RAR; 
XCH 
LDM 
ADD 
XCH 
OUT: 

FIN 
INC 
FIN 
INC 
FIN 
I NX 
FIN 
LD RD; 
XCH RE 
LDM 
XCH 
XCH 
JMS 
JMS 
FIM 
SRC 
LDM 
WRM ; 
BBL 



RO ; 

i; 

0; 

RF; 

RE; 

RO ; 

NZCY 

0; 

Rl ; 



NZCY 

OUT; 

8 ; 

Rl ; 
NCR; 



RD 

A; 

Rl 

Rl 

R4 
Rl 
R6 
Rl 
R8 
Rl 
R ' 



0; 

RC 
RD. 
MULT 
SUBI R; 
RO 2; 
RO ; 

0; 



CR; 



ACR; 
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THIS ROM CONTAINS A 16*16 HEX MULTIPLICATION TABLE. 
EACH VALUE REPRESENTS A MULTIPLICATION OF TWO SINGLE 
HEX-DIGIT NUMBERS. THE ZERO ROW IS USED FOR TABLE 
INSTRUCTIONS; 

*****#*****************************************❖£#*££* 



TABLE: 
ROW 1 : 
ROW2 : 
R0W3 : 
ROW 4 : 
ROW 5: 
ROW 6 : 
ROW 7: 
ROWS : 
ROW 9 : 
ROW A : 
ROWS: 
ROWC : 
ROWO : 
ROWE : 
ROW F : 



FIN 


P.O ; 


BBL 


0 ; 


nop; 




NOP; 




NCP ; 




NOP ; 




nop ; 




nop; 




NOP; 




NOP ; 




mop; 




nop; 




MOP; 




NOP; 




nop; 




nop; 




CCN 


" 00 " 


CON 


" 01 " 


CCN 


" 02 " 


CGN 


" 03" 


CON 


"04" 


CON 


"05" 


CON 


"06" 


CCN 


"0 7" 


CON 


"03" 


CON 


" 05" 


CON 


"OA" 


COM 


"OB" 


CCN 


"OC" 


CON 


» 00 " 


CON 


"OE" 


CON 


"OF" 


CON 


" 00 " 


CON 


" 02 " 


CON 


"04" 


CON 


"06" 


CON 


"08" 


CON 


"OA" 


CON 


"OC" 


CON 


11 OE" 


CON 


" 10 " 


CON 


" 12 " 


CON 


"14" 


CON 


"16" 


CON 


"18" 


CON 


"1A» 


CON 


" 1 C" 


CON 


"IE" 


CON 


" 00 " 


CON 


"03" 


CON 


"06" 


CON 


"09" 


CON 


"OC" 


CON 


"OF" 


CON 


" 12 " 


CON 


"15" 


CON 


"18" 


CON 


"IB" 


CON 


"IE" 


CON 


" 21 " 


CON 


"24" 


CON 


"27" 


CON 


"2 A" 


CGN 


" 20 " 


CON 


" 00 " 


CON 


"04" 


CON 


"08" 


CCN 


"OC" 


CON 


" 10 " 


CON 


"14" 


CON 


"3.8" 


CON 


"1C" 


CON 


" 20 " 


CON 


"24" 


CON 


"28" 


CON 


»2C" 


CON 

CON 


"30" 


CGN 


ii 34 H 


CCN 


"38" 


CON 


» 3 C " 


" 00 " 


CON 


11 05" 


CON 


"OA" 


CON 


"OF" 


CON 


"14" 


CON 


"19" 


CON 


"IE" 


CON 


"2 3" 


CON 


"28" 


CON 


" 20 " 


CON 


"32" 


CCN 


"37" 


CON 


«3C" 


CON 


"41" 


CON 


"46" 


CON 


" c, b » 


CON 


" 00 " 


CCN 


"06" 


CON 


"OC" 


COM 


" 12 " 


CON 


"18" 


CON 


"IE" 


CON 


"24" 


CON 


"2 A" 


CON 


"30" 


CON 


"36" 


CON 


II 1. Q II 


CON 


"42" 


r n m 
v^. j • ' 


"43" 


CON 


11 (, t:n 


C n N 


II 94. II 


CON 


USA " 


CON 


" 00 " 


CON 


"07" 


CON 


"OE" 


CON 


"15" 


CON 


" 1 C" 


CON 


"23" 


CON 


" 2 A " 


CON 


"31" 


CCN 


"38" 


CON 


II -3 p II 


CCN 


"46" 


CON 


"40" 


CCN 


"54" 


CON 


" 5B" 


CON 


"62" 


CON 


"69" 


CON 


" 00 " 


CON 


"08" 


COM 


" 10 " 


CCN 


"18" 


CON 


" 20 " 


CON 


"28" 


CON 


"30" 


CON 


"38" 


CON 


"40" 


CON 


"48" 


COM 


"50" 


CON 


"5 8 " 


CON 


"60" 


CCN 


"63" 


CON 


"70" 


CON 


"78" 


CON 


" 00 " 


CON 


"09" 


CON 


" 12 " 


CON 


"IB" 


CON 


"24" 


CON 


" 20 " 


CON 


"36" 


CON 


"3F" 


CON 


"48" 


CON 


"51" 


CCN 


"5 A " 


CON 


"6 3" 


CON 


" 6 C " 


CON 


"7 5" 


CON 


"7E" 


CON 


"87" 


CON 


" 00 " 


CON 


"OA" 


CON 


"14" 


CON 


"IE" 


CON 


"28" 


CON 


"32" 


CON 


"3C" 


CON 


"46" 


CON 


"80" 


CON 


"5 A" 


CON 


"64 " 


CCN 


" 6 E" 


CON 


"73" 


CON 


"82" 


CON 


" 8 C " 


CON 


"96" 


CCN 


" 00 " 


CON 


"OB" 


CON 


"16" 


CON 


" 21 " 


CON 


" 2 C" 


CUN 


"37" 


CON 


"42" 


COM 


"4 0" 


CON 


"58" 


CON 


" 63" 


CON 


" 6 E " 


CON 


11 7 q it 


CCN 


"84" 


COM 


" 8 F" 


CON 


"9A" 


CON 


" A 5 " 


CON 


"GO" 


CON 


"OC" 


CON 


"18" 


CON 


"2 4" 


CON 


"30" 


CON 


11 30 ," 


CON 


"48" 


CON 


" 54 11 


CON 


" 60 " 


CON 


" 6 C" 


CON 


"76" 


CON 


"84" 


CCN 


"90" 


CON 


"9C" 


COM 


" A 8 " 


CON 


"04 " 


CON 


" 00 " 


CON 


" 00 " 


CON 


" 3 A" 


CON 


11 p 7 n 


COM 


"34" 


CON 


"41" 


COM 


"4 L " 


CON 


"5 8 " 


CON 


" 68 " 


CON 


11 75 " 


CON 


"62" 


CON 


"OF" 


CON 


hoc" 


COM 


11 /\ 9 11 


CON 


"P.6" 


C n M 


"03" 


COM 


" 00 " 


COM 


"Of:" 


CON 


" 1 C. " 


CON 


"2 A" 


CCN 


"38" 


COM 


"46" 


CON 


II C. 11 


COM 


"62" 


CON 


"70" 


CON 


"YE" 


CON 


" 6 C " 


CON 


" 9 A " 


CON 


"AO" 


CON 


" p 5 " 


CON 


Iir / t 11 


COM 


H 02 " 


CON 


" 00 " 


COM 


» or" 


CON 


j'c" 


CON 


"2 0 " 


CON 


ti ■> H 


CON 


" 9 p " 


ton 


"5 A" 


CON 


"6 9" 


CON 


"70" 


COM 


"87" 


r 0 


"94" 


CCN 


" A 5 " 


CON 


"84 " 


CON 


"C3" 


CON 


" 02 " 


CON 


"El" 
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$ # s* 3!-. if * if af: $ $ £ # « a* # sit Sp 5?: X; # fc £ a* * a* * $ ❖*$*#**# * 

THIS SUBROUTINE MULTIPLIES A FOUR HEX-DIGIT NUMBER 
BY A HEX-DIGIT NUMBER RESULTING IN A FOUR HEX-DIGIT 
CHOPPED NUMBER ; 

**#**«*«^J}: 5X*«5{:*3{t^^j{t2j:*«5!5**«^5lc^* £ * ** afrafraj:### # :£=)£:£# 



MULT: 



CLC; 
FIM R2 
SRC R 2 
LDM o; 

wrm; 

LDM 4; 



160; 



XCH 

INC 

INC 

INC 

INC 

INC 

INC 

INC 

INC 

INC 

INC 

INC 



R3 
R 3 
R3 
R3 
R3 
R3 
R3 
R3 
R3 
P3 
R3 
R3 



SRC 


R2 


LD 


R4 


WRM; 


SRC 


R2 


LD 


R5 


WRM: 


SRC 


R2 


LD 


R6 


WRM; 


SRC 


R2 


LD 


R7 


wrm; 


SRC 


R2 


LD 


P 8 


wrm; 


SP.C 


R2 


LO 


R9 


WRM; 


SRC 


R2 


LD 


RA 


WRM; 


SRC 


P2 


LO 


RB 


WRM; 


SRC 


R2 


LD 


RC 


WRM; 


SRC 


R2 


LD 


RO 


WRM; 


SRC 


R2 


LD 


RE 


wrm; 


SRC 


R2 


LO 


RF 


wrm; 



REM ZERO SIGN COUNTER 



REM STORE CONTENTS OF INDEX 
REGISTERS INTO MEMORY; 



REM THIS SECTION OF THE MULTIPLY SUBROUTINE 

TESTS THE TYPE AND SIGN OF THE INPUT VALUES, 



FIM RO 2; 

SRC RO; 

RDM; 

JCN NZ AC COS; 
LDM 15; 

XCH R3 ; 

SRC R2 ; 

r-, r>. t * • 

r\un i 

RAL ; 

JCN ZCY PM; 

LDM 12 ; 

XCH R3 ; 

JMS COMPLEMENT; 
LDM 0; 

XCH R3 ; 

SRC R2; 

LDM 8; 

WRM; 

PM: 

LDM 11 ; 

XCH R3 ; 

SRC R2 ; 

RDM; 

RAL; 

JCN ZCY NM; 

LDM 8; 

XCH R3 ; 

JMS COMPLEMENT; 
LDM 0; 

XCH R3 ; 

SRC R2 ; 

LDM 8; 

ADM; 

wrm; 

NM : 

LDM 8; 

XCH R3 ; 

CLC; 

MULT2: 

SRC R2 ; 

RDM; 

RAL; 

WRM; 

! SZ R3 MULT 2 ; 



REM LOAD COS FLAG; 



REM JUMP IF COS; 

REM TEST SIGN OF FIRST INPUT; 



REM JUMP IF POSITIVE; 

REM COMPLEMENT NEGATIVE; 



REM SET SIGN COUNTER; 



REM TEST SIGN OF SECOND INPUT; 



REM JUMP 1 F POSI TI VE; 

REM COMPLEMENT NEG INPUT; 



REM SET SIGN COUNTER; 



REM NORMALIZE INPUTS; 
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REM THIS BEGINS THE MAIN PORTION OF MULTIPLY SUBROUTINE. 

INDEX REGISTEFS 0 AND 1 ARE USED TO FETCH VALUES FROM 
MULTIPLICATION TABLE. THE TABLE VALUES ARE ADDED AS IN 
LONG-HAND MULTIPLICATION; 



POS1 

LOCI 



POS2 

LOC2 



COS: 




LDM 


15; XCH R3 
ZAC P0S1 ; 


JCN 


XCH 


ro ; 


LDM 


8; XCH R3 


XCH 


Rl ; 

TABLE; 


JMS 


XCH 


RO ; 


XCH 


R4 ; 


JUN 


LOCI ; 


: XCH 


R4 ; 


: LDM 


14; XCH R3 


JCN 


ZAC P0S2 ; 


XCH 


RO; 


LDM 


9; XCH R3 


XCH 


Ri ; 


JMS 


TABLE; 


XCH 


RO; 


XCH 


R5 ; 


JUN 


L0C2; 


: XCH 


R5 ; 


: LD R4; 


ADD 


R5; 


XCH 

TCC; 


R4 ; 



XCH 

LDM 

JCN 

XCH 

LDM 

XCH 

JMS 

XCH 

XCH 

XCH 

XCH 

JUN 

POS3: XCH 
LDM 
XCH 

L0C3: LD I 



R7; 

15; XCH R.3; 
ZAC PQS3; 

R 9; XCH R3; 
Ri ; 

TABLE ; 

Ri ; 

R5 ; 

ro ; 

R6 ; 

LCC3 ; 

R5 ; 

0 ; 

R6 ; 

A • 



ADD R5; 

XCH R4 ; 

TCC; 

XCH R5 ; 

LDM 14; XCH R3; 
JCN ZAC P0S4; 
XCH RO; 

LDM 10; XCH R3 ; 
XCH RI ; 

JMS TABLE; 

XCH RI ; 

XCH R 8 ; 

XCH RO; 

XCH R9; 

JUN L0C4; 

PO S4 : XCH RG; 

LDM 0; 

XCH R9 ; 

LD R6: 

ADD P 7 ; 

ADD R9 ; 

XCH R6; 

TCC; 

XCH R7 ; 



SRC R2 ; RDM; 
SRC R2; RDM; 

SRC R2 *, RDM; 
SRC R2 ; RDM; 



SRC R2; RDM; 
SRC R2 ; RDM; 



SRC R2; RDM; 
SP.C R2 ; RDM; 



l r ^^4 



L0C4 : 





LDM 


15; XCH R3 ; 


SRC 


R2 ; 


RDM 




JCN 


ZAC P0S5; 








XCH 


RO; 










LDM 


10; XCH R3 ; 


SRC 


R2 ; 


RDM 




XCH 


Ri ; 










JMS 


TABLE; 










XCH 


Ri ; 










XCH 


R9 ; 










XCH 


RO; 










XCH 


RA; 










JUN 


L0C5; 








P0S5: 


XCH 


R9 ; 










LDM 


0; 










XCH 


RA ; 








L0C5 : 


LD P 


5; 









lu r< j , 

ADD R6 ; 
XCH R5 ; 
TCC; 



XCH R6 ; 
LO R5; 
ADO R9 ; 
XCH R5; 
TCC; 



ADO R6 ; 
XCH R6 ; 



LDM 

JCN 

XCH 


14 ; 
ZAC 
RO ; 


XCH R3; 
PGS6 ; 


SRC 


R2; 


RDM; 


LDM 


li ; 


XCH R3 ; 


SRC 


R2 ; 


RDM; 



XCH R1 ; 

JMS TABLE; 

XCH PI ; 

XCH RB; 

XCH RO; 

XCH R9 ; 

JUN L0C6; 

P0S6: XCH RB; 

LDM 0; 

XCH P9 ; 

L 0C6 : LD R7; 

ADD R9 ; 

ADD RA; 

XCH R7 ; 

TCC5 
XCH R9 ; 

LD.’-' 15; XCH R3; SRC R2 ; RDM; 
JCN ZAC R0S7 ; 

XCH RO; 

LDM 11 ; XCH R3; SRC R2 ; RDM; 
XCH P.l ; 

JMS TABLE; 

XCH R 1 ; 

XCH PA; 

XCH RO; 

XCH RC ; 

JUN L0C7; 

P0S7 : XCH RA; 

LDM 0; 

XCH RC ; 

L0C7 : LD R9; 

ADD RC ; 

XCH R 9 ; 

LD R6 ; 

ADD R 7 ; 

XCH P6 ; 

TCC; 

XCH R7 ; 

LD P.6; 

ADD RA ; 

XCH f<6 ; 

TCC; 

ADD R7 ; 
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ADO R9 ; 

XCH RF; 

LDM 13; XCH R3; 
JCN 7.AC POS8 ; 
XCH RO; 

LDM 1.0; XCH R3; 
XCH Ri; 

JMS TABLE; 

XCH RO; 

XCH R7 ; 

JUN LCC8 ; 

P0S8: XCH R7; 

L0C8 : LD R4; 

ADD R7 ; 

XCH R4 ; 

TCC; 



XCH 


R7 ; 






LD P.4 ; 






ADD 


R8 ; 






XCH 


R4 ; 






tcc; 








XCH 


R 8 ; 






LDM 


13; 


XCH 


R 


JCN 


ZAC 


P0S9 


t 


XCH 


RO ; 






LDM 


ll ; 


XCH 


R 


XCH 


pi ; 






JMS 


TABLE; 




XCH 


Rl ; 






XCH 


R 9 ; 






XCH 


RO; 






XCH 


RA ; 






JUN 


L0C9 







P0S9: XCH R9 ; 

LDM 0; 

\ / r* I • r\ A « 

ALH r\ 9 

L0C9: LD RA; 

ADD R7? 

ADD R5 ; 

XCH R5 ; 
TCC; 

ADD R6 ; 

XCH R65 
TCC.; 

XCH RA ; 

LDM 12; 

JCN ZAC 
XCH RO; 

LDM ll; 

XCH Rl ; 

JMS TABLE; 
XCH RO; 

XCH P.7: 

JUN loCa; 
POSA: XCH R7; 
LOCA: LD R7; 

ADD R9 : 

XCH R7 ; 
TCC; 

ADD P.5 ; 

XCH R5 ; 
TCC; 

XCH PE; 

LD P.5; 

ADD P.8 ; 

XCH R 5 ; 
ICC; 

ADD RE ; 

ADD R6 ; 

XCH RE; 

TCC ; 

ADD RF ; 



XCH R3; 
POSA; 

XCH R3 ; 



SRC R2 ; RDM; 
SRC R2; RDM; 



SRC R2; RDM; 
SRC R2 ; RDM; 



SRC R2 ; RDM; 
SRC R2; RDM; 



i r ;6 



xch rf; 

LD R4: 
ADD R? ; 
XCH RC; 
TCC; 

ADD R8 ; 
ADD R5 ; 
XCH RD; 
TCC; 

ADD RE; 
XCH RE ; 
TCC; 

ADD ra; 
ADD RF; 
XCH RF; 



CLB ; 

FIM RO 2; 

SRC RO; 

RDM; 

JCN NZ AC CSIN; 
FIM RO 111; 

JMS DIV2IR; 

FIM R2 108; 

FIM RO 168; 

JMS TRAN RAM; 
CLB; 

FIM R2 160; 

SRC R2; 

RDM; 

ral; 

JCN ZCY PPRCD; 

F I H R 2 168; 

JMS COMPLEMENT; 

PPRCD: 

FIM R2 168; 

JMS RAMI RC; 

JUN ENDMULT; 

CSIN: 



LDM 4 ; XCH R3 ; 
XCH R4 ; 

INC R3; 

SRC R2‘, 

RDM; 

XCH R5 ; 

INC R3 ; 

SRC R2 ; 

POM; 

XCH R6; 

INC R3 ; 

SRC R2 ; 

RDM; 

XCH R 7 ; 

ENDMULT: BBL 0; 

END 



REM END TABLE LOOKUP MULT, PRODUCT 
LOCATED IN IR*S "C" THRU "F"; 



REM TEST COS FLAG; 



REM JUMP IF COS; 

REM SHIFT FOR SIGN BIT; 



REM TEST SIGN OF PRODUCT; 

PEM JUMP IF POSITIVE; 

r n un« rnrMT r. n n ru i r T « 

r\ c im l/Unr l era in i r mJiJuc i t 

REM LOAD PRODUCT INTO IRC TO IRF; 



REM LOAD INDEX REGISTERS 4 THRU 7 
WITH INITIAL VALUE BEFORF 
MULT ROUTINE CALLED; 

SRC R2 ; RDM; 
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APPENDIX D 



* SS******** #*** *#*##**#**«#**##«* 

THIS PROGRAM SIMULATES THE MICROCOMPUTER 
IN THE NAVIGATIONAL PROGRAM TO COMPUTE 
THE PRESENT POSITION OF THE AIRCRAFT 

**<-.***:*=***$«* ***sS*5M‘$fc£*S*£#***S**^**** *££*********❖** 



IMPLICIT INTEGER (A-B,D-V) 

DIMENSION R(Ifc),RAM(I6,16), ROM( 16 , 32 ) , Y ( 53 ) , Z ( 52 ) , T ( 4 ) 
DATA XT A S/O. 400/ , XVGVI /-. 2774/ , XVGXI /O. 3236/, XVWY/00.0 
1/, IF/O/ ,DF/0/, XTASR/0.400/, XTH/O. 7654/, XSDD/.01810/, XS 
2HD/. 1418/, XT/O. 347/ 

DO 1 1=1,16 
R ( I ) = 0 
DO 1 J=l, 16 
1 RAM{ I , J )=0 
1=2 

Y( 1)=. 99999 
X=. 03125 
30 Y ( I I = COS ( X ) 

X=X+. 03125 
1 = 1+1 

IFIX.GT. 1.63251 C-0 TO 40 

GG TO 30 
40 DO 60 J=1 , 52 
60 Z( J)=(Y(J)-Y( J+l ))/. 03125 
L = 1 

DO 63 J = l, 13 
DO 63 K = 1 , 7,2 
FIX1 = Y( L 1*65536. 0*16. 0 
FIX2 = Z ( 0*65536 . 0*16.0 
P.EM- F I X 1 — { FIX1/16)*16 
F I XI = F I XI / 16+RCM/ 8 
R EM=F I X 2 — ( FI X2/ 16) *16 
F 1X2= FI X2/ 16+REM/8 
DO 64 1=1,4 

R C M ( J , 4 * ( K - 1 ) + I ) = F I X 1 - ( F I X 1 / 1 6 ) * 1 6 
F I X1 = F I XI / 16 

ROM( J, 4*K+ I ) =F I X 2 - ( FIX2/16) *16 
64 F I X2= F I X2/ 16 
L =L + 1 

63 CONTINUE 

DO 1001 I F = 1 ,2 . , 

DO 1002 DF=1,2 
XTH= 3.1415/4.0 
XTH= XT H+ 0.0067266 
XVGXI =-XVGXI 
XVGYI =- XVGY I 
XLCCF = 0.0 
YLGCF = 0.0 
XLC'CM = 0.0 
YLOCM = 0.0 
DC 100 M= ] ,10 

CALL IJAV(R ,PAM,RGM, T, TE , XTAS, XVGY I , XVC-X I ,XVWY,XVWX, IF 
1 XTASP , XTH,XSDD,XSHD,XT , XL OCF , YLOC F, XLOCM , YLOCM) 

100 CONTINUE 

XTH=3 .1415 + XTH 
XVGYI = - XVGY I 
X VGX! = -XVGXI 
DO 110 N = 1 , 1 0 

CALL NA V ( P , R A M , RDM , T , TE , XT AS , X VGY I , X VG X I , X VWY , XV WX , I F 
1 XTASR . XTH , XSDD , X SHU , XT , XL GCF , Y L C C I , XLOCM , Y L OCM ) 

no continue 

1002 CONTI ML) t 
1001 CONTINUE 
RETURN 
END 
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# Jfc# ## # # 1^3* £ £ ** # :*; * $ ^ # £ ajc £ £ £ £ # # # * # *: £ * ^ # a{c ❖ £ # £ £ ❖ $ $ 5$C $ # * 

THIS SUBROUTINE SIMULATES THE FUNCTION 
OF THE MICROCOMPUTER EXECUTIVE ROUTINE 



SUBROUTINE NAV( R, RAM, ROM, T , TE , XT AS, XVGYI , 

IXVGXI i XVWY, XVWXt I F , DF , XTASR , 

2XTH,XSDD, XSHD,XT , XLOCF , YLOC F, XLOCM , YLOCM ) 

IMPLICIT INTEGER (A-R,D-V) 

DIMENSION R ( 1 6 ) , RAM 16,16), RC'Mt 16,32), T(4), TE<4) 
CALL CONVRT ( 1 , 13, XTASP ,PAM, R,0 ) 

CALL CONVRK 1 , 5 , XT , RAM , R, 0 ) 

CALL CONVRT ( 3 1 1 , XT AS , R AM , R , 0 ) 

CALL CONVRT ( 3, 5 , X T H , R A M , R , 1) 

CALL CONVRT ( 4, 1, XVWX, PAM, R , 0) 

CALL CONVD ( 4 , 1 , P. AM , P , XVA, 0 ) 

WRITE (6,103) XV A 
103 FORMAT ( 10X, ' VKX= • , F15.6) 

CALL CC NVRT ( 4 , 5 , X V W Y , P AM, R, 0) 

CALL CONV RT( 5, 1,XVGXI ,RAM,R,0) 

CALL CONVRT( 5, 5, XVGYI , RAM, R,0) 

CALL CONV RT (6»1,XSHD,RAM,R, 0) 

CALL CONVRT (6, 5, XSDD,RAM,R,0) 

XSIN- -COS (3. 1415/2 .O-XTH) 

XCOS=CCS( XTH) 

WR ITE ( 6 , 2 00 ) XS I N , XCOS 

200 FORMAT! 10X, 1 XS I N= ' , F10 . 6, 2X , « XCOS= ' , F10 . 6 ) 

CALL CO SP L ( P , RDM , RAM ) 

SIN AND COS HAVE JUST BEEN CALCULATED 



CALL CONVD ( 5 » 13 , RAM , R, XVAL1,0) 

CALL CONVD ( 2 , 1 3 , R AM. , R , XVAL2 , 0 ) 

WPITE(6,200) XVAL1 , XVAL2 
X t AS= (XiASR+XTASi/2.0 
XVAY=XTAS*XCCS 
XVAX=XTAS*XSIN 
WRITE(6,201) XT AS , XVAY , XVAX 

201 FORMAT ( 10X, ' XTAS = ' ,F10.3,2X, 'XVAY=* , FI 0.3, • XVAX=' , F10, 
13) 

DO THE ABOVE CALCULATIONS WITH MCS-4 
R ( 3 ) = 3 

P ( 4 ) = 1 

CALL DI V2 ( R, RAM ) 

R ( 3 ) = 1 
R (4)=13 

CALL DI V2 ( R, RAM ) 

P ( 1 ) = 3 
R ( 2 ) = 1 
R < 3) = 1 
R ( 4 )= 13 

CALL ADDRAM(R, RAM) 

R ( 3 ) = 1 
R ( 4 ) = 1 3 

CALL RAMI RC(R, PAM) 

R ( 1 ) = 1 
R { 2 ) = 1 3 
R ( 3 ) = 2 
P. ( 4 ) = 1 3 

CALL MULT PN( R, RAM) 

R ( 3 ) = 4 
R ( 4 ) = 9 

CALL IRP.AMCtR, RAM) 

R (3) =5 
R ( 4 ) = 1 3 

CALL MULT PN ( P , 0 AM ) 

R ( 3 ) = 5 
R ( 4 ) = 9 

CALL I RPAMC( F , RAM ) 

CALL C 0 N V D ( 4 , 9 , R A M , R , X V A L 1 , 0 ) 

CALL CO NVL> ( 6 , 9 , R A M , R , XVAL2 , 0) 
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c 



202 

C 



203 

C 

C 

600 

601 

C 



WRITE (6,201) XT A S , XVAL1 » XVAL2 

WE HAVE JUST COMPARED OUTPUTS VAY AND VAX 
IF(DF.NE.O) GO TO 600 
XDXD=XSDD*XCO$+XSHD*XS IN 
XDYD=XSHD*XCOS-XSDD*XSIN 
WRITE (6,202) XDXD , XDYD 

FORMAT ( 10X, 'XDXD=» ,F10.4,2X, ’ XDYD=' , F 1 0 « 4 ) 

DO THE ABOVE STATEMENTS ON THE • MI CROCOMPUT ER 

R ( 1 ) = 6 
R ( 2 ) = 5 
R ( 3 ) = 2 
R ( 4 ) = 1 3 

CALL MULTPNIR, RAM) 

R (3 ) = 7 
R ( 4 ) = 5 

CALL IRRAMC ( R, RAM ) 

R ( 1 ) = 6 
R ( 2 ) = 1 
R (3)=2 
R ( 4 ) = 1 3 

CALL MUL T PN ( R, RAM ) 

R (3) = 7 
R (4)=1 

CALL I RRAMC ( R, RAM ) 

R (3)=5 
R ( 4 ) = 1 3 

CALL MULTPN(R, P AM) 

R ( 1 ) = 7 
R ( 2 ) = 5 



CALL ADRMTR(R, RAM) 
R ( 3 ) = 6 
P ( 4 ) =9 

C * • I t n o ^ nr / r. DAM' 
ALL 1 .'M-MI 'L VP ) KAf-i | 

R ( i ) = 6 



R ( 2 ) = 5 
R ( 3 ) = 5 
R ( 4 ) = 1 3 



CALL MULTPN(R, RAM) 

R ( 1 ) =7 
R ( 2 ) = 1 

CALL COMPLC ( Rt RAM ) 

CALL ADRM IRC R, RAM) 

R (3 )=3 
R ( 4 ) =9 

CALL I RRAMC (R, RAM) 

CALL CO NV D ( o , 9 , R A M , R , X V A L 1,0) 

CALL CQNVU ( 3 , 9 , R A M , R , X V AL 2 , 0) 

WRITE (6 ,203) X V A L 1 , XV A L 2 

FORMAT ( 10X, ’ DXD= ’ , F10.4,2X, * DVD = ’ ,F10.A ) 

WE HAVE JUST COMPLETED THE CALCULATIONS WHICH 
MAKE USE OF THE DOPPLER RADAR 

IF ( IF. ME. 0) GO TO 300 
XVW Y= ( 3 .0*XVwY+( XVAY-XVGYI ) ) /4. 0 
XVWX= ( 3.0-XVWX+ ( XVAX-XVGXI ) )/4.0 
WRITE(6,60l) xvwy, xvwx 

FORMAT \ iOX, * X VI. Y - • , F10 . 4, 2X, • XVWX" ’ ,F10. 4) 

WE ARE GOING TO PROCEED WITH MCS 4 



R. ( 3 ) - 5 
R ( 4 ) = 6 

CALL RAMI P.C(P, RAM) 
CALL COMP LC ( P , PAM ) 
R ( J ) = 4 
P (2)=9 

CALL ADRM IR(P, PAM) 
R ( 1 ) a 4 
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R ( 2) = 5 

CALL ADRM IR( R, RAM ) 

CALL ADRMIP(R, RAM) 

CALL ADRMIRC R f RAM) 

R ( 3 ) = 4 
R ( 4 ) = 5 

CALL IRRAMCtR, RAM) 

CALL DI V2 ( R, RAM) 

CALL D I V2 ( R» RAM ) 

R ( 3 ) = 5 
R ( 4 ) = i 

CALL PAFIRCtR, RAM) 

CALL COMPLC(R, RAM) 

R ( 1 ) = 5 
R ( 2) =9 

CALL AORMIR(R, RAM) 

R ( 1 ) = 4 
R < 2 ) = 1 

CALL A'DRM1R(R, RAM) 

CALL ADRMIRtR, RAM) 

CALL A DRM I R ( F; , RAM ) 

R (3)=4 
P <4) = 1 

CALL IRRAMC(R, RAM) 

CALL D I V2 ( R» RAM ) 

CALL DI V2 ( R» RAM ) 

CALL CCNVD(4, 5, RAM, R, XVAL 1 , 0) 

CALL C0NVD(4,1 , RAM , R, XVAL2 , 0 ) 

303 FORMAT ( 10X, 'VWY=' , F10.4,2X, ' VWX= ' ,F10.4) 



C 



30i 

C 



304 

C 

300 



C 



WE HAVE COMPLETED CALCULATING 

XDY=XT*XVGYI 

XDX=XT*XVGXI 

WRITE (6,301) XDX,XDY 

FORMAT ( iOX, 5 XDX- ! , F10.5,2X, ' XDY = 



VWY, VWX 



, F10, 



REPEAT THE ABOVE CALCULATIONS ON MCS-4 

R ( 1 ) = 1 
P 1 2 ) = 5 
R { 3 ) = 5 
R { 4 ) = 5 

CALL MULT PN( R, RAM) 

R ( 1 ) = 2 
R ( 2 ) - 9 

CALL DI V2 IR( R, RAM) 

R ( 1 ) = 1 
R ( 2 ) = 5 
R ( 3 ) - 5 
R ( 4 ) = 1 

CALL MULTPN(R, RAM) 

R ( 1 1 = 1 
R < 2 ) = 9 

CALL r>IV2 IR(R, RAM) 

C ALL COMVD ( 1 , 9 , R AM , R , XVAL 1,0) 

CALL COMVD ( 2 , 9 , RAM , R , XVAL2, 0) 

WRITE (6,304) XVAL 1, XVAL 2 
FORMAT ( IOX , ' DX = ' , (10.4 , 2X, ' DY= ' , F10.4) 

WE HAVE COMPLETED CALCULATING DY, DX 

I F( .MOT . ( ( IF .LG.O) . AND. (DF. EQ.O ) ) ) GO TO 400 
XDY=I XDY+XDYD) /? .0 
XGX=( XDX+XDXO) / 2 . 0 
W R I T E ( 6 , 3 0 1 ) XDX , XDY 

REPEAT THE ABOVE SEQUENCE WITH MGS -A 



R ( l ) - 3 
R ( 2 ) = 9 
R ( 3 ) = 2 
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c 

400 

C 



c 

500 

C 



R (4) =9 

CALL ADDRAMtR, RAM) 

CALL DI V2 ( R, RAM ) 

R ( 1 ) =6 
R ( 2 ) =9 
R ( 3 ) = 1 
R ( 4 ) = 9 

CALL AODR AM ( R , RAM ) 

CALL DI V2 ( Rt RAM ) 

CALL C0NVD(2,9,RAM,R,XVAL2,0) 

CALL CONVD ( 1 , 9 , RAM , R , X VAL 1 , 0 ) 

WRITE (6, 304 } XVAL1,XVAL2 

DY AND DX HAVE JUST BEEN CALCULATED 

GO TO 540 

I F ( . N OT « < (IF.EQ.l). AND. ( DF.EO .0) ) ) GO TO 500 

XDY=XDYD 

XOX=XDX D 

WRITE (6,301) XDX,XDY 

REPEAT THE ABOVE SEQUENCE WITH MCS-4 

R ( 3 ) =3 
R ( 4 ) = 9 
R ( .1 ) = 2 
R ( 2) =9 

CALL TRANRMtR, RAM) 

R ( 3 ) =6 
R ( 4 ) =9 
R < 1 ) = 1 
R ( 2 ) = 9 

CALL TRANRMtR, RAM) 

CALL CGNVD(2,9,RAM,R,XVAL2, 0) 

CALL CONVD ( 1 , 9 , RAM , R , X VAL 1 , 0 ) 

! ,p * T TC f ^ *"> r> | v/ W a i ■» \/ W A i n 

W f\ i. ’ C L Uf JU-'H 

DY AND DX HAVE JUST BEEN CALCULATED 
GO TO 540 

IFtUF.NE.D.CR. (DF.NE.l) ) GO TO 540 
XDY=XT* ( XVWY+XVAY ) 

XOX=XT* (XVWX+XVAX) 

WRITE (6,301) XDX , XDY 

REPEAT THE ABOVE CALCULATION WITH MCS-4 

R 1 1 ) -4 

R ( 2 ) == 5 
R (3 )=4 
R ( 4 ) =9 

CALL ADDRAMtR, RAM) 

R ( 1 ) = i 
R ( 2 ) = 5 

CALL MULTPNtR, RAM) 

R ( 3 ) = 2 
R ( 4 ) =9 

CALL IRRAMCtRt RAM) 

R ( 1 ) - 4 
R < 2 ) = 1 
R ( 3 ) = 5 
R ( 4 ) - 9 

CALL ADDRAMtR, RAM) 

R( 1 ) = 1 
R ( 2 ) = 5 

CALL MULT PN(R, RAM) 

R ( 3 ) = 1 
R(4)*9 

CALL I R P A KG ( R, RAM) 

call cn :;v o 1 2 , 9 , u t m , a , x v a l 2 , 0 ) 

CALL CONVD ( 1 ,<>, c AM,R , X VAL 1,0) 

WRITE (O, 304) XV AL 1 , XV AL 2 
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c calculations are now complete for dx and dy 

540 XLOCF = XLOCF + XDX 

XLOCM = XLOCM + XVALL 
YLCCF = YLOCF + XOY 
YLOCM = YLOCM + XVAL2 
WRITE (6t550) XLOCF, YLOCF 
WRITEJ.6, 560) XLOCM. , YLOCM 

550 FORMAT (3X, 'XLOCF = ' , F 10 . 5 , 4X , • YLOCF = ',F10.5) 

560 FORMAT (8X, 'XLOCM = • , F 10 . 5 , 4X , ' YLOCM = • ,F10.5,//) 
RETURN 
END 



I 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE COSINE 



SUBROUTINE COS I NE ( I R , RAM ) 

IMPLICIT INTEGER (A-Z) 

DIMENSION 1 R ( 16 ) , RAMI 16 , 32 ) , T ( 4 ) , TE ( 4 ) 

I = I R ( 16)*S+IP( 15 ) / 2 

J=(IR( 15J-I IRI 15)/2)*2)*16+IR< 14) /8*8 
DO 1 K=5» 12 

1 IR(K)=RAM.I H-l, J + K-4) 

I R( 15 i = 1 R I 1 4 ) - I R I 1 4 ) / 8 * S 
I R ( 14 ) =1 R 1 1 3 ) 

I P ( L3 > = 0 
I R ( 16 ) = 0 
DO 2 K = 1 , 4 
T ( K ) = I R ( 1 2+K ) 

2 T E I K ) = I R( 8+K ) 

CALL MULTI T,TE, IR) 

ACC= I I R I 8 I - 1 R ( 16) )*16**3+( IP.I 7 ) -I Rt 15) ) *16**2+! I R I 6 ) - 1 
1RI14) )*16+( IR(5)-IR( 13) ) 

ACC=ACC/2 
DO 3 K= 1 , 4 

I R I 12 + K)=4CC-( ACC/16)*16 

3 ACC = * CC / lo 
RETURN 
END 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-A SUBROUTINE ADDRAM 

^ # £ *3* ** 5jG * 3* # # £ * * * *r 5jf * * 



jjc afc if & # j;< if. jjc £ # :J; # if. :}t 



if # :Jt rf; if # # :Jc # £ sjc 3*.: ajs j{c 



SUBROUTINE ADOS AM ( IR.RftM) 

IMPLICIT INTEGER (A-V) 

DIMENSION IF ( 16) , RAM I 16, 16) 

S UM = 0 

DO 1 J - 1 , 4 

1 S UM= SUM <-( RAMI I R I 1 ) , IR I 2H-J-1 H RAM( IP. (3 ) , IR I 4) + J-1 ) K16 
1**1 J-l ) 

DO 2 J = 1 , 4 

RAMI I R( 3) • I c* I 4 ) + J-l ) = SUM- ( SUM/1 4) *16 

2 SUM-SUM/ 16 
RETURN 
END 



1(3 3 






THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE SUQIR 

# If * J)C sjc ## Jf! *5)C# * Jjc# « $ # a;:# sjt # sjs j* # # s}c « >}c a*## £ #4:# $ * 



subroutine subriri ir,ram) 

IMPLICIT INTEGER < A- V I 
DIMENSION I R ( 16 > ,RAM< 16,16) 

SUM=0 
M = I R ( 3 ) 

N = IR( 4 ) 

DO 1 J = l,4 

T ERM- ( RAMI M,N + J-1 )-IR( 12+J ) >*16**( J-l) 

1 SUM=SUM+TERM 
DO 2 J=l,4 

RAMI IRI I) , IR< 2 ) +J-l)=SUM-( SUM/ 16 >*16 

2 SUM=SUM / 1 6 
RETURN 
END 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE OIV2IR 

^ ^ ^ ^ > r vV \V ^ *Ar %*# sA# •»*# »>V 4# 4# V* V# U# 4# O* •. V 4# O/ -Jf ^ ^ ^ %t> 

-f(Ai #j\ #1* *p «,k #p v ^p ^p ^p ^p #p ^p ^p -p ^p /p >p ^p #p ■'p ^p #p #p •‘j" ^p -V* "T V V ^ •T* ^p “T ^ -S’* ^ ^f* ^ >p ^p -p V 



SUBROUTINE D I V2 I R I IR, R AM ) 
IMPLICIT INTEGER (A-V) 
DIMENSION IRI16) , RAMI 16,16) 
SUM=0 
F 1=0 

IF! IR( 161 . LT.8) GO TO 3 
F 1= 1 

CALL COMP L C I I R , R A M } 

3 00 1 J = 1 , 4 

TERM= I IRI 12+J) )*16**< J-l) 

1 SDM=SUM+TERM 
SUM=$ ! JM/2 

00 2 J=l,4 

IRI J+12 )=SUM-( SUM/ 16) *16 

2 SUM=SUM/16 
Ir(Fi.EQ.O) GO TO 5 
CALL COMPLC ( I R » R A M ) 

5 DO 6 J=l,4 

6 RAMI IRI 1) , IRI 2 ) + J- 1 ) = I R I J + l 2 ) 
RETURN 
END 



# £ s;s •/*£:& >!o;t # :)f :*« jfcxc sj; ### #£##### sjs # :fc 



;Jc ^ ^ # £ ajc >Jc # # # jJ; ^ 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE AD DRAM I R 

V #5^ *#*>[! £***## Iflflflfififif ##***#* * 






S UBROUT INF A DR M I R ( I R , R A M ) 

IMPLICIT INTEGER (A-V) 

DIMENSION IRI 16) ,RAM( 16,16) 

SUM=0 
K = I R ( 1 ) 

L = I P I 2 ) 

DO 1 J = 1 , 4 

TERM=( R AMI K , L + J - 1 ) + 1 R I 1 2+ J ) ) * 1 6** I J- 1 ) 

1 SUM- SUM+T E RM 
DO 2 J = 1, 4 

IRI 12 + J ) - SUM- ( SUM/ 16) *16 

2 SUM=SUM/16 
RETURN 
END 



❖ £ 2jC # # * # # # # * * * * # # * *5^ # 2{C £ jjt # * # ##** * £ # # $ £ # # #*$ 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE RAMIR8 

jjt v *t;fc*cs}:#y,s#:* # ##££## A#* sj: * $ ## # ## 



SUBROUTINE R AM I R3 ( I R, R AM ) 
IMPLICIT INTEGER (A-V) 
DIMENSION IR( 16 ) ,RAMl 16, 16) 
DO 1 J = 1 , 4 

1 I R ( 8+ J ) =RAM ( I R ( 1 ) , I R ( 2 ) +J-1 ) 
RETURN 
END 



A S{c ^ £$>}:#;$; ;Jc 2^ £ # 3{c £ 3|s i 






THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MC 5-4 SUBROUTINE IRRAMC 



SUBROUTINE IRRAMC ( I R, RAM ) 
IMPLICIT INTEGER (A-V) 
DIMENSION I R ( 16 ) , RAMI 16,16) 

DO 1 J = 1, 4 

1 R AM( I R ( 3 ) , IR(4)+J— 1)=IR(12+J) 
RETURN 
END 



5jc ^ # ijcjjs ;Jc 2^# # £ ?Jt $ # jfc # # 3$C Jjc^js S}: # # 5js >J< sfc # # 2 {c :£ # 2 $: # 3jt :Jc 3}c a}c 

THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE RAM IRC 

y# j# \V yu v* ^ y* wu sV >v y# >v y • y# J/ y# %v ^ ^ x- x x x X x x x x x x x x x x x x x x x x %*+ x 

*,** ^ -**,v #,k >,% ^ v A|\ p ^,k ^ ^ v v # r # r* *r v v p i % *r* t t 



SUBROUTINE RAM I RC ( IR, RAM) 
IMPLICIT INTEGER (A-V) 
DIMENSION IRC 16) , RANK 16,16) 

DO 1 J = 1 , 4 

1 I R( 12*J ) =RAM( I R ( 3) , IR(4 )+J-l ) 
RETURN 
END 



*-:*«*«* *# * *** *#**:Sr* ****** ****** *4'## **««* 

THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE TRANRAM 

=Jn $ * if}): # % % £ if $ f- # $ if. 1-. * Sf :y: $ # <: if # $ <s # # ;*t * & sfr $ $ $ # # # aje $ >jt if if # *,t # a)c * >Jc # 



S U3R0UT I N E T R ANR M (IP., RAM ) 

IMPLICIT INTEGER (A-V) 

D I MENS I ON I R ( 1 6 ) , R AM ( 1 6 , 1 6 ) 

DO 1 J = 1 , 4 

1 R AM ( I R { 1 ) , I M 2 ) *• J - 1 )=RAM( IR(3 ) , IP ( 4)+J - 1 ) 
RETURN 
END 
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THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE 01 V2 

# ijt ^ ^ yjt ^ # $: * * # # # 



SUBROUTINE DIV2(IR,RAM) 
IMPLICIT INTEGER (A-Z) 

0 1 MENS I ON IR ( 16 ) ,RAM( 16,161 

SUM=0 

K= IR ( 3 ) 

L=IR(4) 

F 1 = 0 

DO I J = 1 , 4 

1 IR(J+12)=RAM(K,L+J-1) 
IF(IR(16) .LT.8) GO TO 3 
F 1= I 

CALL COMPLC ( I R, RAM) 

3 DO 4 J = 1 , 4 

4 SUM=SUM+IR( J+12l*16**( J-l) 
SUM=SUM/2 

DO 2 J= I , 4 

I R ( J + 12 ) = SUM-( SUM/ 16) *16 

2 SUM=SUM/16 
IF(FI.EQ.O) GO TO 5 
CALL COMPLC ( IR, RAM ) 

5 DO 6 J = 1 , 4 

6 RAM { K , L+J- 1 ) = I R ( J -*- 1 2 ) 

RETURN 

END 






* ^ vU ^ V# %*# ^ ^ %V 'V ^ vV ^ 'V ^ ^ %'f 'V *V %4U \<# JU X X X X 

* ^ /,v ^ ^ ^,S "|\ #p ^ ^ ^ #,<» 



THIS SUBROUT I NE SI MULATES THE FUNCTIONS 
OF the KCS-4 MULTIPLICATION TA3LE 



at at at at a 






SUBROUTINE TABLE ( A, 8, C » D) 

IMPLICIT INTEGER (A-Z) 

D=(A*8) / 1 6 

C=A*B-D*16 

RETURN 

END 



# y,c ^ ^ 5lc 5r =r ^ :£ ^;>|s # # :*c # ^ ?;? * # # ijssjs £ :‘,c 2 ^ 3* ( c sjc # * # >}c # >1' # ❖ sjc^c* 

THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE AOORAM 



# 3$; :£ £ 3 ;: # # : 



: # s$c £ # lie 5|c if. £ 2(c sjc # # jfc ^ # # sjc £ s{c :J? # # sjc £ # # '4s # A' & # ajs ^ # £ $ 



SUBROUTINE ADD { A , B , C , 0 ) 

IMPLICIT INTEGER (A-Z) 

D = ( A + B ) / 1 6 

C-(A+B)-D*16 

RETURN 

END 
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THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE MULT 

3}t £ $ £ Jjt Jjc if 5f: # & £ & £ ■$ $ # & # % ^ ^ # £ # * £ + 5.V * :jc :£ $ * £ # j.'c # ;Jc s|; £ * * jjc £ sjc * 4: ^c jjc 



1 



12 

145 



21 



28 



38 

2 



SUBROUTINE MULT!X,Y,R) 

IMPLICIT INTEGER (A-Z) 

D I MENS I ON X(4),Y(4),R( 16) , TEMP ( 16 ) 
DO 1 1=1,16 
TEMPI I) =R( I ) 

CALL TABLE ( X ( 1 ) , Y ( 4 ) , BL , R ( 5 ) ) 

CALL TABLE! X!2),Y(3)»BL,R(6)) 

CALL ADD( R ( 5 ) , R( 6 ) ,R! 5 ) ,R( 8 ) ) 

CALL TABLE(X(2) ,Y!4),R!6) ,R(7) ) 
CALL ADD(R( 6) ,R( 5 ) ,R( 5 ) ,R(6) ) 

CALL TA3LE(X(3) , Y ( 3 ) , R ( 9 ) ,R(10) ) 
CALL A D D ( R ( 7 ) , K ( 8 ) ,R(8),BL) 

CALL ADD(R(8),R( 10) ,R(7),R(8) ) 

CALL TABLE(X(3),Y(4),R!10),k!ll)) 
CALL ADD( R ( 6) , R ( 7 ) , R ( 6 ) , R < 7 ) ) 

CALL ADD(R(10),R(6),R(6),CY) 

CALL A D D ( C Y , R ( 7 ) ,R(7), BL) 

CALL TABLE (X( 4) , Y( 3), R ( 12) , R( 10 ) ) 
CALL ADD(R(8) ,R( 10) ,R(8), BL ) 

CALL ADD(R(8) , BL , R l 8 ) , BL) 

CALL A D 0 ( R { 8 ) , R ( 1 1 ) , R ( 8 ) , R ( 10 ) ) 
CALL TABLE ( X ( 4 ) , Y ( 4 ) , P ( 1 1 ) , R C 1 3 ) ) 
CALL ADD(R(13),R( 10)»P(10) ,BL) 

CALL ADD! R < 7 J , R ( 8 ) , R ! 7 5 , R ( G ) } 

CALL ADDlR(ll) , R ( 7 ) , R ! 7 ) , CY ) 

CALL ADD(CY,R(8) , A C C » B L ) 

CALL ADD( ACC, R( 10 ) , R( 16 ) , BL ) 

CALL TABLE! X(3 ) , Y( 2 ) , BL,R(8 ) ) 

CALL ADD(R(8),R(5),R(5) , R ( 8 ) ) 

CALL ADD! R(9 ) ,R( 5) ,R(5 ) ,R(9) ) 

CALL TABLE (X(4) , Y( 2) , R( 10) ,R(1 1 ) ) 
CALL ADD(Rdl) ,R!8) ,ACC,BL) 

CALL ADD(ACC,R(6) ,R(6 ),CY) 

CALL ADD! CY,R (7) , R( 7) , R( 1 1) ) 

CALL TABLE! X(4 ) ,Y( 1),BL,R(8)) 

CALL ADO!R( 10) ,R( 8) ,R(8 ) , CY ) 

CALL ADD! CY , R ! 6 ) , R(6),P(15) ) 

CALL ADD! R(6) ,R! 12 ), R !6 ) , CY ) 

CALL ADD(CY,R{ 15) ,ACC, BL) 

CALL ADD! ACC , P. ( 7 ) , R ( 1 5 ) , CY ) 

CALL ADD(CY,R(16) , R ( 1 6 ) , B L ) 

CALL A 0 D ( R ! 3 ) , R ! 5) ,R( 13) ,CY) 

CALL AOD(CY,R (9 ) , R (9 ) , BL ) 

CALL A DO ( R ( 6 ) , P. ( 9 ) , R ( 1 4 ) , C Y ) 

CALL ADD! R( 15 ) , CY, R! 15) ,CY) 

CALL ADD! CY, R( 11 ) , ACC , CY) 

CALL ADD(ACC+CY,R!16),R(16) ,BL) 

DO 2 1=1,12 
R ! I ) = T E MP ! I ) 

RETURN 

END 
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## ###### £#£:$;#£ £ # ££#:}; ##:*** £* ### 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 ASSEMBLER DUMP AND TRACE 






4 £* #*❖#:{:** 



S USROUT I N E PR I NT ( N , IR , R AM , ROM , K , L ) 

IMPLICIT INTEGER (A-Z) 

DIMENSION I R ( 1 6 ) , RAMI 16,16) ,ROM( 16,32) 

IF(N.NE.l) GO TO 5 
WRITE (6, 101) IR 

101 FORMAT! 5X, ' IR' ,2X, ' 1 2 3 4 5 6 7 8 9 1 0 1 1 12 13 14 1 

19X,9 12,713,/) 

RETURN 

5 IFIN.NE.2) GO TO 10 

WRITE (6, 102) K , L , I R AM ( K , L + 1 - 1 ) , 1 = 1,4) 

102 FORMAT I 5X, 'RAM LOC ' , 12 , 1 3 , 4 1 3 , / ) 

RETURN 

10 I F ( N . N E ♦ 3 ) RETURN 

WRITE (6, 103) K,L, (R0M(K,L+I-1) ,1=1,8) 

103 FORMAT ( 5X , 'ROM LOC ' , I 2 , I 3 , 8 I 3, / ) 

RETURN 

END 



£ ## # # i*x 3&C 3^' :£ $■ 3^ # # ^ # £ $ Xfi # ^ X' ❖ ❖ ££ # 



4^ 4^ 4^ 



^ vj^ 4> 4/ 4j jV jV jJ# \U 4^ 4» ^ 



^U4^ J- X 



Jf 4/ \ 

V -s* - 



THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 SUBROUTINE COMPLEMENT 



# % ajc* ;$c 4c # sjojt tp. 5Jcj£ i 



c £ £ :$r £ s|c :£ sjcsjt sjisjc^ 



SUBROUTINE COMPLC ( R, RAM) 

IMPLICIT INTEGER ( A-Z ) 

DIMENSION R I 16 ), RAMI 16, 16) 

I FI (R 1 16) .LT.O) .OR. (R I 15) . LT.O) .OR. (R( 14) .LT.O) .OR. (R( 
1G0 TO 5 
L = 1 

DO 1 1=1,4 

R I 12+ I ) = 1 5— R (12 + 1 )+L 
I F ( R( 1 2 + ! ) .ME. 16 ) GO TO 3 
R ( 12* I ) =0 
L = 1 

GO TO 1 
3 L = 0 

1 CONTINUE 
RETURN 

5 DO 2 I = 1 , 4 

2 R ( 12 + 1 ) =-R ( 12 + 1 ) 

RETURN 

END 
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THIS SUBROUTINE SIMULATES THE FUNCTIONS 
OF THE MCS-4 EXPANDED MULTIPLY 

❖ # ❖ £ ❖ # — # * « sjs £ # # # * j|: * rt * * sjs * ,-|5 £ * £ * * £ * * * $ $ * * $ $ # jjc * $ * * £ # £ £ jfc £ # $ 



SUBROUTINE MULT PN ( R.R AM) 

IMPLICIT INTEGER! A-Z) 

DIMENSION R(16) ,RAM(16,16) ,T<4) ,TE(4) 

F 1 = 0 

F2=0 

DO 1 1=1,4 

1 R ( 12+ 1 ) = R A Vi ( R ( 1),R(2)+I— 1) 

IF(P( 16) .LT.8) GO TO 2 
F 1 = 1 

. CALL COMPLC ( R, RAM) 

2 SUM=R (16) 

DO 3 1=1,3 

3 SUM=SUM*16+R(i6-I ) 

SUM=SUM*2 

• DO 4 1=1,4 
T ( I )=SUM-( SUM/ 16 )*16 

4 SUM=SUM/16 
DO 11 1=1,4 

11 R( 12+1 ) =RAM(R(3) , R(4)+I-l ) 

IF(R( 16) .LT. 8) GO TO 12 
F2=l 

CALL COMPLC (R, RAM) 

12 SUM=R (16) 

DO 13 1=1,3 

13 SUM=SUM* 1 6+R ( 16-1 ) 

SUM=SUM*2 

DO 14 1=1,4 

TE ( I J --SUM- (SUM/ 16) *16 

14 SUM=SUM/16 

CALL MULT { T,TE»R ‘ 

SUM=R { 16) 

DO 23 1=1,3 

23 SUM=SUM*16+R( 16-1 ) 

SUM=SUM/2 

SUM=SUM/2 

DO 24 1=1,4 
R ( 1+ 1 2 ) =S U M- ( SUM / 1 6 ) * 1 6 

24 SUM=SUM/16 

IF( (F1 + F2) . N E . 1 ) GO TO 30 
CALL COMPLC! R, RAM ) 

30 RETURN 
END 



^Jc 5^ # * * * # * ❖ sj: ?' r - * sjc £ 3^ £ * * 3* # * £ 

THIS SUBROU T I KE SIMULATES THE FUNCTIONS 
OF THE. TH QUADRANT TEST IN THE 
MC5-4 EXECUTIVE ROUTINE 

# 3^ # T* S?# 3(4 £ 4= * 3^ 3*t 3j= # # # $ # # ### Jjc # s(c # # 3(t 3^3^ # £ £ 3,’: 3 (c S}e * 3J; 



SUBROUTINE COSPL ( R , PGM , RAM ) 

IMPLICIT INTEGER ( A-Z ) 

DIME NS I ON R ( 16) , ROM! 16 , 3 2 ) , R AM ( 1 6 , 16 ) 

C PLACE PI/2 INTO REGISTERS 

R ( 13) =2 
R ( 14) =2 
R ( 15) =9 
R ( 16) =1 

C LOCATION OF TH AND INTERMEDIATE STORAGE 

P. ( 3 ) = 3 
R ( 4 ) = 5 
R( 1)=10 

R ( 2 ) = 1 

CALL SUBR IR(?.,PAM) 

I f ! (RAM ! 10 , l ) .LT .0) .OR . (RAM! 10, 2 ) .LT .0) .OP . (RAM! 1 0, 3) . 
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1 <RAM( 10,4) .LT. 0) ) GO TO 1 
R ( 3 ) = 10 
R (4)=1 
R ( 2)=5 

CALL SU8RIR( R,RAM) 

I F ( (RAM { 10,5) .LT .0) . OR . ( RAM ( 1 0 , 6 ) . LT . 0 ) .OR. (RAM (10 t 7). 
1 ( R AM ( 10,3) .LT.O) ) GO TO 2 
R (4) =5 
R ( 2) = 1 

CALL SUBR IR( R , RAM ) 

I F ( (RAM (10,1) .LT.OJ.QR. (RAM (10, 2) .LT.O) .OR. (RAM(10,3). 
1 (RAM( 10,4) .LT.O) ) GO TO 3 
C TI-PI/2 IS NOW IN R AM ( 10 , 1 ) 

R ( 4) = 5 

C TRANSFER PI/2 TO RAM(10,5) 

CALL IRRAMC(Rt RAM) 

R ( 4 ) = 1 

CALL RAM I RC ( R , RAM ) 

CALL CCS I NE( R , ROM ) 

R ( 3 ) = 5 
R (4)=13 

CALL IRRAMCtR, RAM) 

C TRANSFER TH- PRIME TO R(13) 

R ( 3 ) = 1 0 
R ( 4) = 1 

CALL RAMI RCtR, RAM) 

R ( 4)=5 
R ( 1 ) = 1 0 
R ( 2 ) = 1 

CALL SUBRIRt R,RAM ) 

C COMPLEMENTARY ANGLE IS NOW IN RAM(ltlO) 

R (4 ) = 1 

CALL RAMI RC(R, RAM) 

CALL COSINE (R , ROM ) 

R ( 3) =2 
R ( 4 ) = 1 3 

CALL I R R AMC ( P , R AM ) 

C COS(TH) HAS NOW BEEN STORED 

RETURN 
3 R ( 3 ) = 1 0 
R ( 4 ) = 1 

C COMPLEMENT RAM (10,1) 

CALL RAMIRC(P, RAM) 

CALL CGMPLC ( R , RAM ) 

CALL COSINE(P., ROM) 

R ( 3 ) = 5 
R ( 4 ) = 1 3 

CALL 1 RRAMCtR, RAM) 

C SIN(TH) HAS BEEN FOUND 

P ( 3 ) = 10 
R<4)=5 

CALL RAMI RC(R, RAM) 

CALL C0SI NE ( R , ROM ) 

CALL CCMPLCtR, RAM) 

R ( 3 ) = 2 
R(4)= 13 

CALL I RRAMCtR, RAM) 

C COS(TH) HAS BEEN FOUND 

RETURN 
2 R ( 3 ) = 1 0 

R ( 4 ) = 5 

CALL RAMI RC(R, RAM) 

CALL COMPLC ( R , R A“. ) 

CALL C OS I NE ( R , ROM ) 

CALL CHMPLC(ftt RAM) 

P. 13 )- 2 

R ( 4) = 1 3 

C/LL 1 RR AMC ( P , P AM ) 

C CCS IN E HAS BEEN FOUND 

R <31 = 10 
P (4)=1 

CALL RAMI P.C( R, RAM ) 
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CALL COSINE! R» ROM) 

CALL COMPLC! R, RAM ) 

R ( 3 ) = 5 
R ( 4) = 13 

CALL I RRAMC (R » RAM ) 

C SIN(TH) HAS BEEN FOUND 

RETURN 
1 R ( 3 )= 10 

R(4)=l 

CALL RAMIRClRt RAM) 

CALL CGMPLC(R,RAM) 

CALL CCS I NE ( R» ROM ) 

CALL COMPLC(R, RAM) 

R ( 3 ) = 5 
R ( 4 ) = 1 3 

CALL I RRAMC (Rt RAM) 

C SIN(TH) HAS BEEN FOUND 

R (3 1=3 
R (4)=5 

CALL RAMI P.C( R» RAM ) 

CALL COSINE ( R » ROM ) 

R( 3)=2 
R ( 4 ) = 1 3 

CALL ! RRAMC (R. RAM) 

C COS(TH) HAS BEEN FOUND 

RETURN 
END 
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SUBROUTINE CGNVD ( N, M. RAM, R, XVAL 
IMPLICIT INTEGER ( A-V) 

DIMENSION RAM ( lt>,16),R< 16) 

F LAG=0 
R ( 3 ) =N 
R ( 4 ) = M 

CALL RA MIR.C( R » RAM ) 

IF(R( 16) .LT.8) GO TO 3 
FLAG= 1 

CALL COMPLC ( R» RAM ) 

3 K 1=K- 1 

I F ( K i ) 8,9,10 

8 XVAL= 0. 0 
GO TO 4 

9 XVAL=R( 16) 

GO TO 4 

10 X VAL=R ( 16) 

DO 1 I = 1 , K 1 

1 XVAL=R( 16-1 HXVAL*16.0 

4 K 1 =4— K 

IF(Kl.LE.O) GO TO 6 
DO 2 I = 1 » K 1 

2 XVAL = XVAL + R{ 13+K1-T ) * 1 6 . 0** ( - I ) 
IF (K.EC.G) XVAL = X VA L*2 . 0 

I F ( FLAG - EO • 0 ) RETURN 

XVAL- -X VAL 

RETURN 

END 



v i 

• ' j 
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THIS SUBROUTINE CONVERTS THE DECIMAL INPUT 
INTO HEXADECIMAL FOR THE MCS-4 PROGRAM 

^c a^^c^t^cafc^r ##3,* ^ # £ ;Jc $ •£. sfc sfc ?Jc :£ afc A # ^ :£ # 3^ 5*c # sjc 3fc $ $ # # 



3 

1 

5 



2 

10 



SUBROUTINE CONVRT (N,M, XV AP , RAM, R,K) 
IMPLICIT INTEGER(A-V) 

DIMENSION RAMI 16, 16) ,R( 16) 

FLAG=0 

IFIXVAR .GE.O.O) GO TO 3 
XVAR=— X VAR 
FLAG=1 
VAR=XVA R 

IFIK.LT .1 J GO TO 5 
DO 1 1=1, K 

R AM ( N , M+3-K+ 1 ) = VAR— VAR/ 16* 1 6 
VAR=VAR/16 
K 1=4— K 

IF(Kl.LE.O) GO TO 10 

FR=XVAR 

XFR=XVAR 

DO 2 I = 1 , K 1 

XFR=(XFR-FR)*16 

RAMIN, M+4-K-I )=XFR 

P R = X F R 

IFIFLAG .EQ.O) RETURN 
XVAR=-X VAR 
R I 3 ) = N 
R ( 4 ) = M 

CALL RAMI RC I R » RAM ) 

CALL COMPLCI R, RAM ) 

CALL IRRAMCIR, RAM) 

DCTIIOM 

i « i_ i wvU 

END 



172 



BIBLIOGRAPHY 



1. Adams, J.R. and others, The Preliminary Design of a 
Long Range Navigation System for Tactical Aircraft , 

Group Thesis Project 0910, United States Naval 
Postgraduate School, Monterey, 1970. 

2. AGARD Conference Proceedings No. ^3, Inertial Naviga- 
tion-Systems and Components , May 1968. 

3. AGARD-LS 52, Guidance and Control of Tactical Missiles, 

May 1972. ~ 

AGARD-AG-158, Computers in the Guidance and Control of 
Aerospace Vehicles , February 1972. 

5. Altman, L. , "Single-Chip Microprocessors Open Up A 
New World Of Applications," Electronics, V. *17, 

p. 81-87, 18 April 197^ . 

6. Beck, G.E. , and others. Navigation Systems, Van Nostrand 
Reinhold, 1971. 

7. Busacker, R.G., and Saaty, T.L., Finite Graphs and 
Networks: An Introduction with Applications, McGraw- 
Hill , 1965. 

8. Carnahan, B. , Luther, H.A., and Wilkes, J.O., Applied 
Numerical Methods , Wiley, 1969. 

9. Cushman, R.H., "Microprocessors are Changing Your Future. 
Are You Prepared:," EDN, V. 18, p. 26-32, 5 November 
1973. 

10. Cushman, R.H., "Understanding the Microprocessor is No 
Trivial Task," EDN, V. 18, p. ^2-^9, 20 November 1973- 

11. Cushman, R.H., "Understand the 8-bit P: You'll see a 
lot of it," EDN, V. 19, P. 48-5^, 20 January 197 1 !. 



12. Cushman, R.H., "Don't overlook the J l-bit PP: They're 
here and they're cheap," EDN , V. 19, p. M-50, 20 
February 197^. 

13. Cushman, R.H., "The Intel 8080: First of the Second- 
Generation Microprocessors," EDN, V. 19, p. 30-36, 

5 May 197^. 



173 



14. Cushman, R.H., "Microprocessors are rapidly gaining on 
Minicomputers," EDN , V. 19, p. 16-20, 20 May 1974. 

15. Forsyth, G.E., "Pitfalls In Computation, or Why a Math 
Book is not Enough," American Mathematics Monthly , 

p. 931-956, November 1970. 

16 . Holt, R.M. , "Current Microcomputer Architecture," 
Computer Design , p. 6 5— 7 3 > February 1974. 

17. Intel, MCS-4 Micro Computer Set , 1973- 

18. Intel, MCS-8 Micro Computer Set , 1973- 

19. Intel, PL/M Programming Manual , 1974. 

20. Johnson, E.E. and Johnson, J.R., Graph Theory with 
Engineering Applications , Ronald, 1972. 

21. Kayton, M. , Fried, W.R., and others, Avionics Naviga - 
tion Systems , Wiley, 1969 

22. Keele, R.V. , Microprocessor Trade-Off Study for Project 
2175 , Tentative and Unpublished NELC Technical Note, 

San Diego, California, 6 March 1973. 

23 . Klass, P.J., "New Gyro Nears Operational Use," Aviation 
Week , p. 50-52, 19 June 1972. 

24. Klass, P.J., "New Guidance Technique Being Tested," 
Aviation Week , V. 100, p. 48-51, 25 February 1974. 

25 . Koenig, H.E. , Tokad, Y. , and Kesavan, H.K. , Analysis 
of Discrete Physical Systems , McGraw-Hill , 1967- 

26. Lockheed S-3A Avionics , p. 4-20, Aircraft Engineering, 
January 197-4. 

27. KcCalla, T.R. , Introduction to Numerical Methods and 
FORTRAN Programming , Wiley, 1964 . 

28. McCracken, D.D., and Dorn, W.S., Numerical Methods 
and FORTRAN Programming , Wiley, 1964 . 

29. NASA Technical Report R-329, A New Concept in Strapdown 
Inertial Navigation , by J.E. Bortz, Sr., March 1971. 

30. Naval Air Development Center, Navigation FP- 1 06 , 1-973- 

31. Naval Air Systems Command OI -85 ADF-2-10.1, legrated 
Weapons System Theory for A 6 E Aircraft, 1 971 . 



174 



32 . 



33 . 



34. 



35 . 



36 . 



37 . 



38. 



39 . 

40. 

41 . 

42. 



Ogdin, J.L. , "Survey of Microprocessors Reveals Limit- 
less Variety," EDN, V. 19, p. 38-43, 20 April 1974. 

Patrol Squadron Thirty-One, P3C Electronic Systems 
Operational Training Study Guide , 1970. 

Ralston, A., A First Course in Numerical Analysis, 
McGraw-Hill , 1965. 

Reyling, G. , "Performance and Control of Multiple 
Microprocessor Systems," Computer Design, p. 81-86, 
March, 1974. 

Ringo, R.L., "Cost-of-Ownership Design Philosophy for 
Inertial Navigators," Astronautics and Aeronautics, 
p. 59-62, June 1973. 



Schultz, G.W. and Holt, R.M. , "MOS LSI Minicomputer 
Comes of Age," Proceedings 1972 Fall Joint Computer 
Conference , p . 1069-1080 , October 1972. 



Schultz, G.W. , Holt, R.M. , and McFarland, H.L., "A 
Guide to Using LSI Microprocessors," Computer, 
p. 13-19, June 1973. 



Sperry Rand Corporation, 



10 O n 4 - -1 ion o 1 "P r,* O vn 

4. *“■ vy vpbi ui/xLiiuu. j. ^ * * y 



*1070 

— > 1 ~ • 



The University of Michigan Engineering Summer 
Conferences, Modern Navigation Systems , 

Wiener, H. , "Computers that fit in Your Pocket," 
Computer Decisions , V. 5 , p. 8-13, August 1973. 

Weissberger , A. J. , "M0S/LSI Microprocessor Selection," 
Electronic Design, p. 100-104, 7 June 1974. 



.175 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 2231^ 

2. Library, Code 0212 
Naval Postgraduate School 
Monterey, California 9 39^0 

3. Mr. W. McMahon 
P3C Software 

Naval Air Development Center 
Warminster , Pennsylvania 

4. Assoc. Professor U. R. Kodres, Code 53Kr 
Department of Computer Science 

Naval Postgraduate School 
Monterey, California 939^0 

5. Chairman, Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 939^0 

6. Asst. Professor G. A. Kildall, Code 53Kd 
Department of Computer Science 

Naval Postgraduate School 
Monterey, California 939^0 

7. Lt. W. L. McCracken, USN 
136 Barber's Point Rd. 

NAS Alameda 

Alameda, California 9^501 



Copies 

2 

2 

1 

1 

1 

1 

1 



176 



5 MOV 76 



2 3 5 '■> O- 

U 1 



2 U 9 



1 o 



tv 



Thesis 
KL825 McCracken 
1 r\oQi a 



153415 



c .1 



aa e n 

Design study of an 

avionics navigation 
microcomputer 



5 



: 7 r 

MOV 76 



3 5 9 0. 

U 9 U 1 

3 9 u ti 



, n f*C T 




Thesis 1 C ^-:1S 

M1825 McCracken 

Design study of an 
avi oni cs navigation 
mi crocomputer . 



thesM1825 

Design study of an avionics navigation m 




3 2768 001 88499 2 
DUDLEY KNOX LIBRARY 



'i 



